DEV/scripts/ppr_makkink.R

In springgbv/Open-Bodem-Index-Calculator: Calculate the Open Bodem Index (OBI) Score

# make makkink table

  # load csv file with makkink per crop per month
  dt <- fread('dev/data/makkink.csv', header = TRUE)
  
  # colnames for the 12 months
  cols <- colnames(dt)[-1]
  
  # replace all missing values with the value 0.36 for braak soil
  dt[,c(cols) := lapply(.SD,function(x) fifelse(is.na(x),0.23,x)),.SDcols = cols]
  
  # write file to data
  crops.makkink <- copy(dt)
  usethis::use_data(crops.makkink,version = 3, overwrite = TRUE, compress = 'xz')