Hi,

I have the following data 

data1 <- data.frame(count = c(0,1,1,2,4,5,13,16,14), weeks = 1:9,
                    treat=c(rep("1mg",3),rep("5mg",3),rep("10mg",3)))
and I am using 

library(splines)

to fit

glm.m <- glm(count~bs(weeks)+as.factor(treat),family=poisson,data=data1)

and I am interested in predicting the count variale for the weeks 10, 11 and 
12 with treat 10mg and 15mg.

Thanks in advance.

Regards,
Samuel


      
	[[alternative HTML version deleted]]

