# NOT RUN {
data(yates.oats)
head(yates.oats)
newyates <- nna(yates.oats, trait="Y")
head(newyates)
plot(newyates$Y, newyates$nnx)
cor(newyates$Y, newyates$nnx)
#### now fit the models and compare #####
m3 <- mmer2(fixed=Y ~ V*N, random = ~ B + B:MP,
data = yates.oats)
yates.oats$res <- residuals(m3)
m4 <- mmer2(fixed=Y ~ V*N + nnx, random = ~ B + B:MP,
data = newyates)
newyates$res <- residuals(m4)
# library(lattice)
# wireframe(res~row*col,yates.oats)
# wireframe(res~row*col,newyates)
# }
Run the code above in your browser using DataLab