# NOT RUN {
 library(minque)
 data(brate)
 head(brate)
 brate$Geno=factor(brate$Geno)
 brate$Pos=factor(brate$Pos)
 brate$Rep=factor(brate$Rep)
 res=lmm(Brate~1|Geno*Pos+Rep,data=brate)
 res$Var
 res$FixedEffect
 res$RandomEffect
 res=lmm.jack(Brate~1|Geno*Pos+Rep,data=brate,JacNum=10,JacRep=1,ALPHA=0.05)
 res$Var
 res$PVar
 res$FixedEffect
 res$RandomEffect
 ## end
# }
Run the code above in your browser using DataLab