## Fit model on oats data.
data(oats.data)
## Fit simple model with only fixed effects.
LMM1 <- LMMsolve(fixed = yield ~ rep + gen,
data = oats.data)
## Obtain table of effective dimensions.
summ1 <- summary(LMM1)
print(summ1)
## Obtain table of variances.
summ2 <- summary(LMM1,
which = "variances")
print(summ2)
Run the code above in your browser using DataLab