## Fit model on john.alpha data from agridat package.
data(john.alpha, package = "agridat")
## Fit simple model with only fixed effects.
LMM1 <- LMMsolve(fixed = yield ~ rep + gen,
data = john.alpha)
## Obtain deviance.
diagnosticsMME(LMM1)
Run the code above in your browser using DataLab