Strain<-rep(c("X","Y"),c(8,8))
Strain<-factor(Strain)
Sex<-rep(c("M","F","M","F"),rep(4,4))
Sex<-factor(Sex)
Status<-rep(c("E","C","E","C","E","C","E","C"),rep(2,8))
Status<-factor(Status)
Tumours<-rep(c("Y","N"),8)
Tumours<-factor(Tumours)
counts<-c(12,4,74,5,12,2,84,3,14,4,80,10,14,1,79,3)

mice.fit1<-glm(counts~Strain+Sex+Status+Tumours)