sir<-read.csv("CHAP14_Sirindhornia_count.csv")
modelPoiss<-glm(sir$pollinaria ~ sir$flowers*sir$species-1,family="poisson")
summary(modelPoiss)