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))