powered by
data(swe07)
pop
deaths
sex
age
log.pop
data(swe07) fit <- glm(deaths ~ offset(log.pop) + sex * as.factor(age), family = poisson, data = swe07) drop1(fit, test = "Chisq") ## Proportional hazards?
Run the code above in your browser using DataLab