data(lizards)
glm(cbind(grahami, opalinus) ~ height + diameter + light + time,
family = binomial, data=lizards)
brglm(cbind(grahami, opalinus) ~ height + diameter + light + time,
family = binomial, data=lizards)Run the code above in your browser using DataLab