data("salamander")
if (FALSE) {
HLfit(cbind(Mate,1-Mate)~TypeF+TypeM+TypeF*TypeM+(1|Female)+(1|Male),
family=binomial(),data=salamander,method="ML")
# equivalent fo using fitme(), but here a bit faster
}
Run the code above in your browser using DataLab