treatment <- factor(rep(c(1, 2), c(43, 41)), 
                    levels = c(1, 2),
                    labels = c("placebo", "treated"))
improved <- factor(rep(c(1, 2, 3, 1, 2, 3), c(29, 7, 7, 13, 7, 21)),
                   levels = c(1, 2, 3),
                   labels = c("none", "some", "marked"))    
numberofdrugs <- rpois(84, 10)+1 
   
healthvalue <- rpois(84, 5)   
y <- data.frame(healthvalue, numberofdrugs, treatment, improved)
test <- glm(healthvalue ~ numberofdrugs+treatment+improved, y, family=poisson)
summary(test)