2012-11-30 27 views
7

Ông có thể xin vui lòng cho tôi biết làm thế nào để có được/tính RMSE giá trị (root có nghĩa là lỗi vuông) trong R khi bạn thực hiện một mô hình hiệu quả hỗn hợpgốc có nghĩa là lỗi vuông trong R - mẫu hiệu ứng pha trộn

Data: na.omit(binh) 
     AIC  BIC logLik 
    888.6144 915.1201 -436.3072 

Random effects: 
Formula: ~1 | Study 
     (Intercept) Residual 
StdDev: 3.304345 1.361858 

Fixed effects: Eeff ~ ADF + CP + DE + ADF2 + DE2 
       Value Std.Error DF t-value p-value 
(Intercept) -0.66390 18.870908 158 -0.035181 0.9720 
ADF   1.16693 0.424561 158 2.748556 0.0067 
CP   0.25723 0.097524 158 2.637575 0.0092 
DE   -36.09593 12.031791 158 -3.000046 0.0031 
ADF2   -0.03708 0.011014 158 -3.366625 0.0010 
DE2   4.77918 1.932924 158 2.472513 0.0145 
Correlation: 
    (Intr) ADF CP  DE  ADF2 
ADF -0.107        
CP -0.032 0.070      
DE 0.978 -0.291 -0.043    
ADF2 0.058 -0.982 -0.045 0.250  
DE2 -0.978 0.308 0.039 -0.997 -0.265 

Standardized Within-Group Residuals: 
     Min   Q1   Med   Q3   Max 
-2.28168116 -0.45260885 0.06528363 0.57071734 2.54144168 

Number of Observations: 209 
Number of Groups: 46 
+1

Điều này phải ở trên crossvalidated.com. –

+7

@ RomanLuštrik: không phải chỉ là về cách tính * RMSE, đúng không? (Trái ngược với "biện pháp giả giả R^2 thích hợp cho các mô hình hỗn hợp", chắc chắn sẽ được dùng để xác định chéo và được thảo luận tại http://glmm.wikidot.com/faq ... –

Trả lời

12

Bạn don 't cho biết chi tiết về những gì chức năng bạn sử dụng để làm mô hình của bạn, nhưng họ có xu hướng lưu trữ dư của họ bằng cách sử dụng cùng tên, mà bạn có thể kiểm tra với str(), và RMSE có thể dễ dàng tính từ dư:

#make a model 
library(nlme) 
r <- lme(conc ~ age, data=IGF) 

#get the RMSE 
r.rmse <- sqrt(mean(r$residuals^2)) 

và trong bình luận dưới đây, Ben Bolker chỉ ra rằng các đối tượng được thực hiện bởi mô hình fitt các hàm ing nên có phương thức còn lại, làm cho nó có thể thực hiện điều này (mặc dù một số kiểu mô hình có thể trả về số dư đã được chuyển đổi):

r.rmse <- sqrt(mean(residuals(r)^2)) 
+3

nói chung, tốt hơn các hàm mô hình hóa được lưu trữ phải thực hiện phương thức 'residuals()', mặc dù bạn phải cẩn thận một chút trong một số trường hợp về loại số dư nào được trả về (thu nhỏ, Pearson, deviance, studentized ...) –

+0

+1 Điểm tốt. Thêm vào câu trả lời – MattBagg

Các vấn đề liên quan