tests/testthat/test-obic_farm.R

In OBIC: Calculate the Open Bodem Index (OBI) Score

# make an example set for five fields
dt <- OBIC::binnenveld[ID <=5]
dt[,B_DRAIN := FALSE]
dt[,B_AREA_DROUGHT := FALSE]

test_that("obic_farm works with extra indicators when not using classic obi", {
  expect_no_condition(obic_farm(dt, useClassicOBI = FALSE))
  out <- obic_farm(dt, useClassicOBI = FALSE)
  
  expect_equal(
    names(out),
    expected = c('fields','farm'),
    tolerance = 0.01
  )
  expect_equal(
    length(out$field),
    expected = 39,
    tolerance = 0.01
  )
  expect_equal(
    length(out$farm),
    expected = 2,
    tolerance = 0.01
  ) 
  expect_equal(
    nrow(out$farm$indicators),
    expected = 32,
    tolerance = 0.01
  )
  expect_equal(
    nrow(out$farm$scores),
    expected = 6,
    tolerance = 0.01
  )
  expect_equal(
    ncol(out$farm$indicators),
    expected = 5,
    tolerance = 0.01
  )
  expect_equal(
    ncol(out$farm$scores),
    expected = 5,
    tolerance = 0.01
  )
  expect_equal(
    unique(out$farm$indicators[,S_OBI_NFIELDS_LOW + S_OBI_NFIELDS_MEDIUM + S_OBI_NFIELDS_HIGH]),
    expected = 5,
    tolerance = 0.01
  )
  expect_equal(
    unique(out$farm$scores[,S_OBI_NFIELDS_LOW + S_OBI_NFIELDS_MEDIUM + S_OBI_NFIELDS_HIGH]),
    expected = 5,
    tolerance = 0.01
  )
  expect_equal(
    out$field$S_C_OBI_A,
    expected = c(0.86,0.833,0.825,0.897,0.764),
    tolerance = 0.01
  )
  expect_equal(
    c(sum(out$farm$indicators$S_OBI_NFIELDS_LOW),
      sum(out$farm$indicators$S_OBI_NFIELDS_HIGH)),
    expected = c(41, 89),
    tolerance = 0.01
  )
  expect_equal(
    c(sum(out$farm$scores$S_OBI_NFIELDS_LOW),sum(out$farm$scores$S_OBI_NFIELDS_HIGH)),
    expected = c(9, 10),
    tolerance = 0.01
  )
})

test_that("obic_farm works with classic setting", {
  dtclassic <- copy(dt)
  dtclassic[,B_DRAIN := NULL]
  dtclassic[,B_AREA_DROUGHT := NULL]
  out <- obic_farm(dtclassic, useClassicOBI = TRUE)
  
  expect_equal(
    names(out),
    expected = c('fields','farm'),
    tolerance = 0.01
  )
  not.deprecated.columns <- names(out$fields)[!names(out$fields) %in% c('I_E_NGW', 'I_E_NSW')]
  expect_equal(
    ncol(out$field[,..not.deprecated.columns]),
    expected = 35,
    tolerance = 0.01
  )
  expect_equal(
    length(out$farm),
    expected = 2,
    tolerance = 0.01
  ) 
  expect_equal(
    nrow(out$farm$indicators[indicator %in% not.deprecated.columns]),
    expected = 28,
    tolerance = 0.01
  )
  expect_equal(
    nrow(out$farm$scores),
    expected = 6,
    tolerance = 0.01
  )
  expect_equal(
    ncol(out$farm$indicators),
    expected = 5,
    tolerance = 0.01
  )
  expect_equal(
    ncol(out$farm$scores),
    expected = 5,
    tolerance = 0.01
  )
  expect_equal(
    unique(out$farm$indicators[,S_OBI_NFIELDS_LOW + S_OBI_NFIELDS_MEDIUM + S_OBI_NFIELDS_HIGH]),
    expected = 5,
    tolerance = 0.01
  )
  expect_equal(
    unique(out$farm$scores[,S_OBI_NFIELDS_LOW + S_OBI_NFIELDS_MEDIUM + S_OBI_NFIELDS_HIGH]),
    expected = 5,
    tolerance = 0.01
  )
  expect_equal(
    out$field$S_C_OBI_A,
    expected = c(0.86,0.833,0.825,0.897,0.764),
    tolerance = 0.01
  )
  expect_equal(
    c(sum(out$farm$indicators[indicator %in% not.deprecated.columns, S_OBI_NFIELDS_LOW]),sum(out$farm$indicators[indicator %in% not.deprecated.columns, S_OBI_NFIELDS_HIGH])),
    expected = c(34, 77),
    tolerance = 0.01
  )
  expect_equal(
    c(sum(out$farm$scores$S_OBI_NFIELDS_LOW),sum(out$farm$scores$S_OBI_NFIELDS_HIGH)),
    expected = c(8,14),
    tolerance = 0.01
  )
})

test_that('obic_farm is sensitive to changing B_FERT_NORM when useClassicOBI == FALSE',{
  dt.fnorm1 <- copy(dt)
  dt.fnorm05 <- copy(dt)
  dt.fnorm05[,B_FERT_NORM_FR := 0.5]
  
  out.fnorm1 <- obic_farm(dt.fnorm1, useClassicOBI = FALSE)
  out.fnorm05 <- obic_farm(dt.fnorm05, useClassicOBI = FALSE)
  
  expect_false(all(out.fnorm1$fields$I_E_GW_NLEA == out.fnorm05$fields$I_E_GW_NLEA))
})