data(carrots)
str(carrots)
plot(success/total ~ logdose, data = carrots, col = as.integer(block))
## Classical glm
Cfit0 <- glm(cbind(success, total-success) ~ logdose + block,
data=carrots, family=binomial)
summary(Cfit0)
## Robust Fit (see help(robglm)): %%>>>>> FIXME <<<<
Run the code above in your browser using DataLab