# NOT RUN {
data(CPdata)
#### look at the data and fit the model
head(CPpheno)
mix1 <- mmer2(Yield~1,
random=~ id
+ Rowf + Colf,
rcov=~units,
data=CPpheno)
summary(mix1)
# compute the variogram
vm0 <- variogram(mix1, xcoor = "Row", ycoor = "Col")
# plot the variogram
plot(vm0$F1)
# }
Run the code above in your browser using DataLab