Chicago Temps

january = read.table("january.txt",header=T,row.names=1)

par(mfrow=c(2,1),mar=c(2,4,3,0))

plot(january$MAXF~jitter(january$DAY),xlab="DAY",ylab="Maximum ˚F",main="Chicago O'Hare January Temperatures, 1958-2014",ylim=c(-20,80))
points(january$MAXF[which(january$YEAR == 2014)],col='red',cex=3,type='b',pch=16)
text(15,65,"Jan-2014-Highs",col="red",cex=2)

plot(january$MINF~jitter(january$DAY),xlab="DAY",ylab="Minimum ˚F",ylim=c(-20,80))
points(january$MINF[which(january$YEAR == 2014)],col='blue',cex=3,type='b',pch=16)
text(15,65,"Jan-2014-Lows",col='blue',cex=2)