# NOT RUN {
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