library(here)
bbch <- paste0(0:9, "x")
crops <- c(
"Beans (field + vegetable)",
"Peas",
"Summer oilseed rape", "Winter oilseed rape",
"Tomatoes",
"Spring cereals", "Winter cereals")
EFSA_washoff_2017 <- matrix(NA, length(crops), length(bbch),
dimnames = list(Crop = crops, BBCH = bbch))
EFSA_washoff_2017["Beans (field + vegetable)", ] <-
c(NA, 0.6, rep(0.75, 2), rep(0.8, 5), 0.35)
EFSA_washoff_2017["Peas", ] <-
c(NA, 0.4, rep(0.6, 2), rep(0.65, 5), 0.35)
EFSA_washoff_2017["Summer oilseed rape", ] <-
c(NA, 0.4, rep(0.5, 2), rep(0.6, 5), 0.5)
EFSA_washoff_2017["Winter oilseed rape", ] <-
c(NA, 0.1, rep(0.4, 2), rep(0.55, 5), 0.3)
EFSA_washoff_2017["Tomatoes", ] <-
c(NA, 0.55, rep(0.75, 2), rep(0.7, 5), 0.35)
EFSA_washoff_2017["Spring cereals", ] <-
c(NA, 0.4, 0.5, 0.5, rep(0.65, 3), rep(0.65, 2), 0.55)
EFSA_washoff_2017["Winter cereals", ] <-
c(NA, 0.1, 0.4, 0.6, rep(0.55, 3), rep(0.6, 2), 0.4)
save(EFSA_washoff_2017,
file = here("data/EFSA_washoff_2017.RData"))