Plots of Peak emergence for Brendens poster:

# plot of cycle length per month:
x11(width=4,height=5)
mat<-(matrix(1:4,2,2))
layout(mat)
layout.show(4)

#x11(width=3,height=3)
par(mar=c(3,2,2,1))
par(mgp=c(3,0.5,0))
par(cex=.9)
par(tcl=-0.3)

boxplot(split(peak$n,peak$mo),ylim=c(0,12),ann=F)
mtext("Length (Days)", side=2,font=2,line=1.2,adj=0.3)
mtext("Month",side=1,font=2,line=1.5)

#Cycle length for all peak days:
hist(peak$n,breaks=seq(from=4,to=12,by=1),main="",axes=F)
axis(2)
axis(1,at=c(4.5,5.5,6.5,7.5,8.5,9.5,10.5,11.5,12.5),labels=c(4,5,6,7,8,9,10,11,12))
mtext("Length (Days)", side=1,font=2,line=1.6)
mtext("Frequency",side=2,font=2,line=1.2,adj=0.4)

#x11(width=3,height=3)
hist(peak$mo,ann=F,ylim=c(0,50),breaks=seq(from=0,to=12,by=1),main="",axes=F)
axis(2)
axis(1,at=c(0.5,1.5,2.5,3.5,4.5,5.5,6.5,7.5,8.5,9.5,10.5,11.5),labels=c(1,2,3,4,5,6,7,8,9,10,11,12))
mtext("Frequency",side=2,font=2,line=1.2,adj=0.4)
mtext("Month",side=1,font=2,line=1.5)



#Cycle length for JA peak days:

hist(peakJA$n,breaks=seq(from=4,to=12,by=1),ylim=c(0,60),main="",axes=F)
axis(2)
axis(1,at=c(4.5,5.5,6.5,7.5,8.5,9.5,10.5,11.5,12.5),labels=c(4,5,6,7,8,9,10,11,12))
mtext("Length (Days)", side=1,font=2,line=1.6)
mtext("Frequency",side=2,font=2,line=1.2,adj=0.4)