আমি নিম্নলিখিতটি করতে চাই:
1) বিটা ওলস রিগ্রেশন (কোনও দণ্ডের শর্ত নেই) ; হ'ল রিগ্রাস করার জন্য ব্যবহৃত ভেরিয়েবলগুলি। আমি এটি দ্বারা জে
lm.model = lm(y~ 0 + x)
betas = coefficients(lm.model)
২) দন্ডের মেয়াদ সহ লাসো রিগ্রেশন, নির্বাচনের মানদণ্ড বায়েশিয়ান ইনফরমেশন মাপদণ্ড (বিআইসি) হইবে, প্রদত্ত
যেখানে হ'ল ভেরিয়েবল / রিগ্রসর নম্বর, পর্যবেক্ষণের সংখ্যার জন্য এবং পদক্ষেপ 1 এ প্রাপ্ত প্রাথমিক বিটাগুলির জন্য)। আমি এই নির্দিষ্ট মানের জন্য রিগ্রেশন ফলাফল পেতে চাই , যা প্রতিটি প্রতিবেদকের জন্য আলাদা। সুতরাং যদি তিনটি ভেরিয়েবল হয় তবে তিনটি পৃথক মান থাকবে ।টি বি ∗ জ λ জ λ জ
এরপরে ওএলএস-লাসো অপ্টিমাইজেশন সমস্যাটি দেওয়া হয়
লার্জ বা গ্ল্যামনেট প্যাকেজটি দিয়ে আমি কীভাবে এটি আর করতে পারি? আমি ল্যাম্বদা নির্দিষ্ট করার কোনও উপায় খুঁজে পাচ্ছি না এবং আমি যদি চালাই তবে সঠিক ফলাফল পাওয়া যায় কিনা তা আমি 100% নিশ্চিত নই
lars.model <- lars(x,y,type = "lasso", intercept = FALSE)
predict.lars(lars.model, type="coefficients", mode="lambda")
আমি এখানে কোন সাহায্যের প্রশংসা করি।
হালনাগাদ:
আমি এখন নিম্নলিখিত কোড ব্যবহার করেছি:
fits.cv = cv.glmnet(x,y,type="mse",penalty.factor = pnlty)
lmin = as.numeric(fits.cv[9]) #lambda.min
fits = glmnet(x,y, alpha=1, intercept=FALSE, penalty.factor = pnlty)
coef = coef(fits, s = lmin)
লাইন 1 এ আমি আমার নির্দিষ্ট পেনাল্টি ফ্যাক্টর ( λ j = লগ ( টি ) দিয়ে ক্রস বৈধতা ব্যবহার করি), যা প্রতিটি নিবন্ধকের জন্য পৃথক। লাইন 2 ফিটস.সি.ভি. এর "ল্যাম্বদা.মিন" নির্বাচন করে, এটি ল্যাম্বডা যা ন্যূনতম গড় ক্রস-বৈধকরণ ত্রুটি দেয়। লাইন 3ডেটাতেলাসো ফিট ()সম্পাদনকরে। আবার আমি শাস্তি ফ্যাক্টর ব্যবহৃতλ। লাইন 4লাইন 2 এ নির্বাচিত"অনুকূল" - এর সাথে সম্পর্কিতফিটগুলির সহগগুলি বের করে।alpha=1
এখন আমার কাছে রেজিস্ট্রারদের জন্য বিটা সহগ রয়েছে যা মিনিমাইজেশন সমস্যার অনুকূল সমাধান চিত্রিত করে
পেনাল্টি ফ্যাক্টর সহ । সহগের সর্বোত্তম সেটটি সম্ভবত রেজিস্ট্রারগুলির একটি উপসেট যা আমি প্রাথমিকভাবে ব্যবহার করেছি, এটি লাসো পদ্ধতির একটি পরিণতি যা ব্যবহৃত রেজিস্ট্রারগুলির সংখ্যা হ্রাস করে।
আমার বোঝার এবং কোডটি কি সঠিক?
$\alpha$
হয়ে যায় । দয়া করে এটি তৈরি করুন, কারণ এটি আপনার প্রশ্নটি আরও সহজেই মানুষকে বুঝতে সক্ষম করে তুলবে এবং তাই এর উত্তর দিন।