5 + 6[1] 11
A first example of generating a report using R, RStudio and Quarto. The goal of this report is to:
5 + 6[1] 11
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.009427
3rd Qu.:1979 3rd Qu.: 9.25 3rd Qu.:1979 3rd Qu.: 0.150000
Max. :2011 Max. :12.00 Max. :2012 Max. : 0.880000
this just names the data “giss” and gives you a summary.
##Calculations
Now, let’s calculate the mean temp deviation per year, and make a vector of years:
annualTemp = tapply(giss$TempAnom, giss$Year, mean)
years = unique(giss$Year)And let’s create some plots
First, a histogram of the distribution of all temperature anomaly values.
hist(giss$TempAnom, breaks = 20)And now a line plot showing the annual temperatures
plot(years, annualTemp, type='l', lwd=2, col=2, xlab="Years", ylab="Temp Anom")
abline(h=0, lty=2)To show code but not run it….
plot(years, annualTemp, type='l', lwd=2, col=2, xlab="Years", ylab="Temp Anom")
abline(h=0, lty=2)Other useful options:
You can also stack all of these to use multiple