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))
  
})