# class 11 examples for correlation
# loader function for Albany, NY Central Park
# station data (monthly mean anomalies)
source("scripts/loadano.R")



# Albany station data: USW00014735
# only January, July anomalies
alb1<-loadano(station="USW00014735",month="Jan",start=1950,end=2012)
alb7<-loadano(station="USW00014735",month="Jul",start=1950,end=2012)


# New York Central Park station data: USW00094728
nyc1<-loadano(station="USW00094728",month="Jan",start=1950,end=2012)
nyc7<-loadano(station="USW00094728",month="Jul",start=1950,end=2012)



# create a 2x2 panel of figures on one page
# to reset the plot window to a 1-panel plot use:
# par(mfrow(c(1,1)))

par(mfrow=c(2,2))
# top left: time series Albany  and NY Central Park January
plot(c(1950,2012),c(-6,6),typ='n',xlab='year',ylab='Tavg [C]',main="January")
lines(alb1$time,alb1$ano,
     typ='b',pch=1,col=4)
lines(nyc1$time,nyc1$ano,
      typ='b',pch=3,col=2)
plot(c(1950,2012),c(-6,6),typ='n',xlab='year',ylab='Tavg [C]',main="July")
lines(alb7$time,alb7$ano,
     typ='b',pch=1,col=4)
lines(nyc7$time,nyc7$ano,
      typ='b',pch=3,col=2)
# bottom left: scatter plot Albany January with NY Central Park 
plot(c(-6,6),c(-6,6),typ='n',xlab='Alb. Tavg [C] ',ylab='NYC Tavg [C]',main="January")
points(alb1$ano,nyc1$ano,pch=3)
# bottom right: scatter plot Albany July with NY Central Park July
plot(c(-6,6),c(-6,6),typ='n',xlab='Alb. Tavg [C] ',ylab='NYC Tavg [C]',main="July")
points(alb7$ano,nyc7$ano,pch=3)