tests/testthat/test_breathtest_parameters.R

In breathtestcore: Core Functions to Read and Fit 13c Time Series from Breath Tests

test_that("Bluck-Coward data for testing of parameter functions", {
  cf = data.frame(
    method = rep(c("nls", "bayesian")),
    group = rep(c("lean", "obese"), each = 2),
    k = c(0.576, 0.606, 0.529, 0.608),
    beta = c(5.24, 5.79, 5.95, 7.54),
    t12 = c(3.67, 3.63, 4.23, 3.99),
    t12in = c(2.076, 2.110, 2.422, 2.466),
    tlag = c(2.88, 2.88, 3.34, 3.26),
    tlagin = c(1.632, 1.724, 1.92, 2.101)
  )

  t50_maes_ghoos = t50_maes_ghoos(cf)
  expect_lt(max(abs(t50_maes_ghoos - cf$t12) / cf$t12), 0.014)

  t50_bluck_coward = t50_bluck_coward(cf)
  expect_lt(max(abs(t50_bluck_coward - cf$t12in) / cf$t12in), 0.035)

  tlag_maes_ghoos = tlag_maes_ghoos(cf)
  expect_lt(max(abs(tlag_maes_ghoos - cf$tlag) / cf$tlag), 0.02)

  tlag_bluck_coward = tlag_bluck_coward(cf)
  # Bluck-Cowards estimate not very exact
  expect_lt(max(abs(tlag_bluck_coward - cf$tlagin) / cf$tlagin), 0.09)

  # Names must be removed
  expect_null(names(t50_maes_ghoos))
  expect_null(names(t50_bluck_coward))
  expect_null(names(tlag_maes_ghoos))
  expect_null(names(tlag_bluck_coward))
  expect_null(names(t50_maes_ghoos))
})