Boxes Over Time

Here’s some sequential box plots for temperature over each year.


#RStudio errors out on large graphs sizes, write it to file instead
png("rplot.png", width = 4096, height = 4096) 
par(mfrow = c(6,5))
for(e in sort(unique(df$year)))
	{
		temper <- subset(df$temp_c, df$year == e)
		dayer <- subset(df$julian_day, df$year == e)
		boxplot(temper~dayer, ylim=yaxis, xlim=xaxis)
		text(x=20, y=0, col = "red", cex=8, e)
		abline(h=0, v=182.5, col = "red")
		abline(h=24, col = "green")
		lines(df$temp_c[df$year == e],lwd=3,col="red")
	}
dev.off()

The green line is right at 24c/75f. The horizontal red line is at 0c/32f and the vertical red line is at 182.5 'days' - about the middle of the year. I'm not sure if the curved line is meaningful, I believe it would represent a histogram; frequency of the observations, though.

Leave a Reply

Your email address will not be published. Required fields are marked *