tests/testthat/test-influence.R

In MultiFrailty: Shared Frailty Regression Models with Inverse Gaussian, Generalized Lindley, and Gamma Frailty Distributions

test_that("influence_frailty calculates leverage, Cook's D, DFFITS", {
  set.seed(123)
  dat <- r_frailty(n = 40, baseline = "weibull", bpar = c(2, 1.5), frailty = "gamma", fpar = c(0.8))
  fit <- fit_frailty(time = dat$time, status = dat$status, baseline = "weibull", frailty = "gamma")

  inf <- influence_frailty(fit)
  expect_equal(length(inf$leverage), 40)
  expect_equal(length(inf$cooks_distance), 40)
  expect_equal(nrow(inf$diagnostics_table), 40)
})