My Converted Project

First I will read in my data and provide a summary of the data.

giss = read.csv("giss_temp.csv")
summary(giss)

      Year          Month          DecDate        TempAnom        
 Min.   :1881   Min.   : 1.00   Min.   :1881   Min.   :-0.860000  
 1st Qu.:1913   1st Qu.: 3.75   1st Qu.:1914   1st Qu.:-0.210000  
 Median :1946   Median : 6.50   Median :1946   Median :-0.050000  
 Mean   :1946   Mean   : 6.50   Mean   :1946   Mean   :-0.009428  
 3rd Qu.:1979   3rd Qu.: 9.25   3rd Qu.:1979   3rd Qu.: 0.150000  
 Max.   :2011   Max.   :12.00   Max.   :2012   Max.   : 0.880000  

Here is the code used to create a graph of temperature anomolies per year with anomolies before 1980 showing up as purple and after 1980 being orange.

allyears = unique(giss$Year)
ann_temp = tapply(giss$TempAnom, giss$Year, mean)
mycols = ifelse(allyears < 1980, "purple", "orange")
plot(allyears, ann_temp, type = 'h', 
     col = mycols, lwd = 3,
     xlab= "Year", ylab = "T anomaly")