xyplot(Mean_Sapflow~Mean_Temperature,
             data=Summarised_new_mean_sapflow,
             col="red",
             pch=19,
             xlab="Temperature (°C)",
             ylab=expression(paste("Sapflow Litres day"^{-1})),
             type=c('p', 'smooth'),
             col.line='blue')

      ##Produce the date labels for the plot
      Date_labels<-c("June", "July",
                     "August", "September",
                     "October")

      ##Insert the labels into the plot
      text(Mean_Sapflow~Mean_Temperature,
           data=Summarised_new_mean_sapflow,
           labels=Date_labels, cex= 0.9, pos=3)

structure(list(Date = structure(c(3L, 2L, 1L, 5L, 4L), .Label = c("August",
"July", "June", "October", "September"), class
    = "factor"), Mean_Humidity = c(17.6073333333333,
    21.8006451612903, 18.3896774193548, 14.822, 11.3486666666667),
        Mean_Radiation = c(263.673333333333, 270.906451612903, 178.98064516129,
        152.233333333333, 93.6), Mean_Temperature = c(70.5613333333333,
        61.3306451612903, 71.7335483870968, 72.2136666666667, 81.743
        ), Mean_Sapflow = c(16.067, 23.3567741935484, 22.9416129032258,
        19.3093333333333, 6.70066666666667)), class = "data.frame", row.names = c(NA,
    -5L))