Displaying 1 result from an estimated 1 matches for "yieldh".
Did you mean:
yield
2009 Apr 14
0
Fitted values and MSE of individual fits in lmList
...f each value is the "yearloc", I also need to know their "id". Is there a way to ad this info?
I have not figured out how to obtain individual mean square errors from an lmList object. However, I can do it with:
res <- as.data.frame(matrix(NA, ncol = 3, nrow = length(levels(yieldh$yearloc)), dimnames = list(NULL, c("yearloc", "MSE", "rep"))))
for (i in 1:length(levels(yieldh$yearloc))) {
dta <- subset(yieldh,yieldh$yearloc==levels(yieldh$yearloc)[i], drop = T)
x <- anova(lm(yield ~ factor(rep) + id, data = dta, na.actio = na.omit))$&qu...