Data<-read.csv("NOAA.csv",sep=";")
Year<-Data[,1]
Anomaly<-Data[,2]
YearRel<-Year-Year[1]

lo1<-lm(Anomaly~YearRel)
Fit1<-predict(lo1)

lo2<-lm(Anomaly~poly(YearRel,3,raw=TRUE))
Fit2<-predict(lo2)

lo3<-lm(Anomaly~cos((YearRel)*pi/36))  #35
Fit3<-predict(lo3)

yscale<-c(-0.5,2)
plot(Year,Anomaly,pch=18,col="gray",ylim=yscale,ylab="")
par(new=TRUE)

#lines(Year,Fit1,type="l",col="red",ylim=yscale,ylab="Anomaly",lwd=3)
#title("NOAA linear fit, R2=0.53")

#lines(Year,Fit2,type="l",col="blue",ylim=yscale,ylab="Anomaly",lwd=3)
#title("NOAA cubic fit, R2=0.53")

lines(Year,Fit3,type="l",col="blue",ylim=yscale,ylab="Anomaly",lwd=3)
title("NOAA cosine fit, R2 = 0.51")