tests/testthat/test1.R

library(rsurv)

n <-  1000
simdata <- data.frame(
  age = rnorm(n),
  sex = sample(c("f", "m"), size = n, replace = TRUE)
) |>
  mutate(
    aft = raftreg(runif(n), ~ age*sex, beta = c(1, 2, -0.5), dist = "lnorm", meanlog = 0, sdlog = 1),
    ph = rphreg(runif(n), ~ age*sex, beta = c(1, 2, -0.5), dist = "lnorm", meanlog = 0, sdlog = 1),
    po = rporeg(runif(n), ~ age*sex, beta = c(1, 2, -0.5), dist = "lnorm", meanlog = 0, sdlog = 1),
    ah = rahreg(runif(n), ~ age*sex, beta = c(1, 2, -0.5), dist = "lnorm", meanlog = 0, sdlog = 1),
    yp = rypreg(runif(n), ~ age*sex, beta = c(1, 2, -0.5), phi = c(-1, -2, 0.5), dist = "lnorm", meanlog = 0, sdlog = 1),
    eh = rehreg(runif(n), ~ age*sex, beta = c(1, 2, -0.5), phi = c(-0.8, -1, 0.2), dist = "lnorm", meanlog = 0, sdlog = 1),
  )

Try the rsurv package in your browser

Any scripts or data that you put into this service are public.

rsurv documentation built on Oct. 30, 2024, 9:08 a.m.