# 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