Sheff_temp <- read.csv("Sheffield_max_temp_data.csv", header = FALSE)
colnames (Sheff_temp) <- c ("Year", "Month", "Max_temp")
Sheff_temp_noNA <- na.omit(Sheff_temp)
par (mfrow = c (1,1))
boxplot (Sheff_temp$Max_temp ~ Sheff_temp$Month,
xlab = "Months", ylab = "Celsius", xaxt = "n",
main = "Average maximum monthly temperature \n in Sheffield, 1883 - 2014",
col = c (rep ( c("dodgerblue", "indianred1", "green3","goldenrod1") , c (2, 3, 3, 3) )),
cex.main = 1.3
)
axis(side = 1, at = 1:12, labels = month.abb,
las = 3, font =3)
abline (v = c (1:12), col = "grey", lty = 2)
abline (h = c (seq(5 , 25, 5)), col = "grey")

plot (Sheff_temp$Year,
Sheff_temp$Max_temp,
type = "l",
col = "blue",
xlab = "year",
ylab = "Average Temperature (°C)",
main = "Average maximum monthly temperature \n in Sheffield, 1883 - 2014")

par (mfrow = c (1,1))