Tôi đã ước tính một mô hình tuyến tính mạnh mẽ R
với trọng số MM bằng cách sử dụng rlm()
gói MASS. `R`` không cung cấp giá trị cho mô hình, nhưng tôi muốn có một giá trị nếu đó là một đại lượng có ý nghĩa. Tôi cũng muốn biết liệu có bất kỳ ý nghĩa nào trong việc có giá trị R 2 cân bằng phương sai tổng và dư giống như cách các quan sát được cân trong hồi quy mạnh hay không. Suy nghĩ chung của tôi là, nếu, vì mục đích của hồi quy, về cơ bản, chúng ta với các trọng số cho một số ước tính ít ảnh hưởng hơn bởi vì chúng là các ngoại lệ theo một cách nào đó, thì có lẽ vì mục đích tính r 2, chúng ta cũng nên đưa ra các ước tính đó. Ước tính ít ảnh hưởng?
Tôi đã viết hai hàm đơn giản cho và R 2 có trọng số , chúng ở bên dưới. Tôi cũng bao gồm các kết quả của việc chạy các chức năng này cho mô hình của tôi được gọi là HI9. EDIT: Tôi đã tìm thấy trang web của Adelle Coster của UNSW đưa ra một công thức bao gồm vectơ trọng số trong việc tính toán cả hai và giống như tôi đã làm, và yêu cầu cô ấy tham khảo chính thức hơn: http: //web.maths. unsw.edu.au/~adelle/Garvan/Assays/GoodnessOfFit.html (vẫn đang tìm kiếm sự giúp đỡ từ Hội Chữ thập Validated về cách giải thích có trọng này 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
Cảm ơn bất cứ ai dành thời gian trả lời này. Vui lòng chấp nhận lời xin lỗi của tôi nếu đã có một số tài liệu tham khảo rất tốt về điều này mà tôi đã bỏ lỡ, hoặc nếu mã của tôi ở trên khó đọc (tôi không phải là một người viết mã).