data(breslow.dat)
bres.fm <- fit.models(list(Robust = "glmRob", Classical = "glm"),
formula = sumY ~ Age10 + Base4*Trt,
family = poisson(), data = breslow.dat)
bres.sum <- summary(bres.fm)
bres.sum
Run the code above in your browser using DataLab