A first example of generating a report using R, RStudio and knitr. The goal of this report is to:
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
annualTemp = tapply(giss$TempAnom, giss$Year, mean)
years = unique(giss$Year)
hist(giss$TempAnom, breaks=20)
plot(years, annualTemp, type='l', lwd=2, col=2,
xlab="Years", ylab="Temp Anom")
abline(h=0, lty=2)