Here is a bit of my data:
这是我的一些数据:
df <- structure(list(farm = c("F1", "F1",
"F1", "F1"), index = structure(c(1L,
1L, 4L, 4L), .Label = c("HT", "Mid-T", "Outside Mid-T", "Outside South"
), class = "factor"), sensorheight = c(1L, 1L, 1L, 1L), yrmonth = structure(c(1491004800,
1496275200, 1491004800, 1496275200), class = c("POSIXct", "POSIXt"
), tzone = "UTC"), monthindex = structure(c(2L, 2L, 2L, 2L), .Label = c("Spring",
"Winter"), class = "factor"), N = c(2, 2, 1, 1), TempC = c(2.06446759259259,
6.68402777777778, 1.32268518518518, 5.63194444444445), sd = c(1.17081824208967,
0.034373246307681, NA, NA), se = c(0.827893518518518, 0.0243055555555567,
NA, NA), ci = c(10.5193845460483, 0.308831365115372, NA, NA)), .Names = c("farm",
"index", "sensorheight", "yrmonth", "monthindex", "N", "TempC",
"sd", "se", "ci"), row.names = c(7L, 9L, 20L, 22L), class = "data.frame")
df <- s