powered by
data(magnesium)
data(magnesium) fit <- glm(cbind(r,m-r)~trial+group,binomial,data=magnesium[1:20,]) anova(fit,test="Chi") summary(fit)
Run the code above in your browser using DataLab