আমি এমএসএস প্যাকেজটিতে R
এমএম ওজন ব্যবহার করে একটি শক্তিশালী রৈখিক মডেল অনুমান করেছি rlm()
। `R`` মডেলের জন্য কোনও মান সরবরাহ করে না , তবে এটি অর্থবোধক পরিমাণের মতো হলে আমি এটি পেতে চাই। আমি আরও জানতে আগ্রহী যে আর -2 মানের যে মোট ও অবশিষ্টাংশগুলি ওজন হিসাবে একইভাবে পর্যবেক্ষণগুলি শক্তিশালী রিগ্রেশনটিতে ওজন করা হয়েছিল তার ওজন রাখার কোনও অর্থ আছে কিনা । আমার সাধারণ চিন্তাভাবনাটি হ'ল, যদি প্রতিরোধের উদ্দেশ্যে, আমরা মূলত ওজনগুলির সাথে কিছু অনুমানকে কম প্রভাব দেয় কারণ তারা কোনও উপায়ে বিদেশী হয়, তবে সম্ভবত আর 2 গণনার উদ্দেশ্যে আমাদেরও সেগুলি দেওয়া উচিত একই অনুমান কম প্রভাব?
আমি এবং ওজনযুক্ত আর 2 এর জন্য দুটি সহজ ফাংশন লিখেছি , সেগুলি নীচে। আমি আমার মডেলটির জন্য এই ফাংশনগুলি চালনার ফলাফলগুলিও অন্তর্ভুক্ত করেছি যার নাম HI9। সম্পাদনা: আমি ইউএনএসডাব্লিউর অ্যাডেল কাস্টারের ওয়েব পৃষ্ঠা পেয়েছি যা এর জন্য একটি সূত্র দেয় যা উভয়ের গণনা গণনা করার ক্ষেত্রে ওজন ভেক্টরকে অন্তর্ভুক্ত করে এবং আমি যেমন করেছিলাম, এবং তাকে আরও একটি আনুষ্ঠানিক রেফারেন্স চেয়েছি : http: //web.maths ma unsw.edu.au/~adelle/Garvan/Assays/GoodnessOfFit.html (এখনও কিভাবে এই ভরযুক্ত ব্যাখ্যা করা উপর ক্রস যাচাই কাছ থেকে সাহায্য খুঁজছেন R 2 ।)R2
SSe
SSt
#I used this function to calculate a basic r-squared from the robust linear model
r2 <- function(x){
+ SSe <- sum((x$resid)^2);
+ observed <- x$resid+x$fitted;
+ SSt <- sum((observed-mean(observed))^2);
+ value <- 1-SSe/SSt;
+ return(value);
+ }
r2(HI9)
[1] 0.2061147
#I used this function to calculate a weighted r-squared from the robust linear model
> r2ww <- function(x){
+ SSe <- sum((x$w*x$resid)^2); #the residual sum of squares is weighted
+ observed <- x$resid+x$fitted;
+ SSt <- sum((x$w*(observed-mean(observed)))^2); #the total sum of squares is weighted
+ value <- 1-SSe/SSt;
+ return(value);
+ }
> r2ww(HI9)
[1] 0.7716264
যে কেউ এই উত্তর দিতে সময় ব্যয় করে ধন্যবাদ। আমি মিস করেছি কিনা সে সম্পর্কে ইতিমধ্যে কিছু খুব ভাল রেফারেন্স থাকলে বা উপরের আমার কোডটি পড়তে খুব কঠিন হয় (আমি কোনও কোড লোক নই) দয়া করে আমার ক্ষমা প্রার্থনা করুন।