library(glmm)
set.seed(1234)
data(salamander)
m<-1000
sal<-glmm(Mate~0+Cross,random=list(~0+Female,~0+Male),varcomps.names=c("F","M"),
data=salamander,family.glmm=bernoulli.glmm,m=m,debug=TRUE,doPQL=FALSE)
coef(sal)
Run the code above in your browser using DataLab