powered by
data(bliss)
data(bliss) attach(bliss) plot(log(dose),r/m,ylim=c(0,1),ylab="Proportion dead") fit <- glm(cbind(r,m-r)~log(dose),binomial) summary(fit)
Run the code above in your browser using DataLab