# Shustack, D.P. and Rodewald, A.D. (2011) Nest predation reduces benefits to early clutch initiation in northern cardinals Cardinalis cardinalis. J. Avian Biol. 42: 204􏰣209. doi: 10.1111/j.1600-048X.2011.05231.x

d<-read.csv("fledge_number.csv")
par(mfrow=c(2,2))
plot(d,xlab="Clutch initiation date",ylab="Number of fledglings",pch=17,col="blue")
freq<-table(d[,2])
#freq
barplot(freq,xlab="Number of fledglings",ylab="Frequency")
mean(d[,2])
var(d[,2])
install.packagres("vcd")
library(vcd)
gf<-goodfit(freq,type= "poisson",method= "MinChisq")
summary(gf)
plot(gf, type = "standing", xlab="Number of fledged chicks",ylab="Frequency", scale = "raw")
plot(gf, type = "hanging", xlab="Number of fledged chicks",ylab="Frequency", scale = "raw")
model_card<-glm(d$fledged ~ d$date,family="poisson")
summary(model_card)
plot(d,xlab="Clutch initiation date",ylab="Number of fledglings",pch=17,col="blue")
abline(model_card,col="orange",lwd=2)