tests/testthat/test_pseudo.R

context("Weighted pseudo observations")

testthat("Pseudo values with weights one give same values as original", {

  ckdat <- generate_data(n = 150, "A")

  test1 <- pseudoci.weighted(ckdat$Tout, ckdat$delta, tmax = 26.5, weights = rep(1, nrow(ckdat)))
  test2 <- pseudoci(ckdat$Tout, ckdat$delta, tmax = 26.5)

  expect_equal(mean(abs(test1$pseudo$cause1 - test2$pseudo$cause1)), 0.0, tolerance = .01)
  expect_equal(mean(abs(test1$pseudo$cause2 - test2$pseudo$cause2)), 0.0, tolerance = .01)

})
sachsmc/pseupersims documentation built on June 12, 2019, 1:14 p.m.