### Biochemical Oxygen Demand data from Marske
BOD <-
  data.frame(Time = c(1, 2, 3, 4, 5, 7),
             demand = c(8.3, 10.3, 19, 16, 15.6, 19.8))
attr(BOD, "reference") <- "A1.4, p. 270"