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