inst/samples/glmEx.R

In proftools: Profile Output Processing Tools for R

Rprof("glmEx.out")
system.time(replicate(100, {
    clotting <- data.frame(
        u = c(5,10,15,20,30,40,60,80,100),
        lot1 = c(118,58,42,35,27,25,21,19,18),
        lot2 = c(69,35,26,21,18,16,13,12,12))
    summary(glm(lot1 ~ log(u), data = clotting, family = Gamma))
    summary(glm(lot2 ~ log(u), data = clotting, family = Gamma))
}))
Rprof(NULL)