tests/testthat/test-rng.R

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

test_that("r_frailty generates valid data frames across all 8 censoring mechanisms", {
  cens <- c("none", "right", "left", "interval", "type1", "type2", "progressive", "progressive_type1")

  set.seed(123)
  for (c_type in cens) {
    dat <- r_frailty(n = 50, baseline = "weibull", bpar = c(2, 1.5),
                     frailty = "gl1", fpar = c(1.2, 0.5),
                     x = matrix(rnorm(50), ncol = 1), beta = 0.5,
                     cen_type = c_type)

    expect_s3_class(dat, "data.frame")
    expect_equal(nrow(dat), 50)
    expect_true(all(dat$time > 0))
    expect_true(all(dat$status %in% c(0, 1, 2, 3)))
  }
})