sir<-read.csv("CHAP14_Sirindhornia_count.csv")
modelPoiss<-glm(sir$pollinaria ~ sir$flowers*sir$species-1,family="poisson")
modellm<-lm(sir$pollinaria ~ sir$flowers*sir$species-1)
par(mfrow=c(2,4))
plot(modelPoiss)
plot(medellm)
AIC(modelPoiss,modellm)