### Cardinal fledging

d<-read.csv("fledge_number.csv")
plot(d,xlab="Clutch initiation date",ylab="Number of fledglings",pch=17,col="blue")
model_card<-glm(d$fledged ~ d$date,family="poisson")
abline(model_card,col="red",lwd=2)


par(mfrow=c(2,2))
model_card2<-glm(d$fledged ~ d$date,family="quasipoisson")
plot(model_card2)
summary(model_card2)

### note that slope and intercept remain unchanged but the p-value has

par(mfrow=c(1,1))
plot(d,xlab="Clutch initiation date",ylab="Number of fledglings",pch=17,col="blue")
abline(model_card,col="red",lwd=2)
abline(model_card2,col="black",lwd=2)