tests/testthat/test-CooksD.R

In MartinSchobben/point: Reading, Processing, and Analysing Raw Ion Count Data

#-------------------------------------------------------------------------------
# minimum functionality (over time with snapshots)
#-------------------------------------------------------------------------------

test_that("residual diagnostics on internal dataset is consistent", {

  # Descriptive an predictive statistics for 13C/12C ratios (note .output
  # argument and remove zero count analysis)
  tb_R <- stat_R(real_IC, "13C", "12C", file.nm, .output = "complete",
                 .zero = TRUE)

  # residual based augmentation of ion count data for isotope ratios
  expect_snapshot(
    CooksD(tb_R, "13C", "12C", file.nm, .X = Xt.pr, .N = N.pr,
           .species = species.nm, .t = t.nm, .output = "flag")
  )
  expect_snapshot(
    Rm(tb_R, "13C", "12C", file.nm, .X = Xt.pr, .N = N.pr,
       .species = species.nm, .t = t.nm, .output = "flag")
  )
  expect_snapshot(
    CV(tb_R, "13C", "12C", file.nm, .X = Xt.pr, .N = N.pr,
       .species = species.nm, .t = t.nm, .output = "flag")
  )
  expect_snapshot(
    QQ(tb_R, "13C", "12C", file.nm, .X = Xt.pr, .N = N.pr,
       .species = species.nm, .t = t.nm, .output = "flag")
  )
  expect_snapshot(
    norm_E(tb_R, "13C", "12C", file.nm, .X = Xt.pr, .N = N.pr,
           .species = species.nm, .t = t.nm, .output = "flag")
  )
  expect_snapshot(
    IR(tb_R, "13C", "12C", file.nm, .X = Xt.pr, .N = N.pr,
       .species = species.nm, .t = t.nm, .output = "flag")
  )
})