inst/test-data/no-gamma.R

In SimSurvNMarker: Simulate Survival Time and Markers

args_env <- new.env()
local(envir = args_env, {
  d_m <- 2L
  d_g <- 5L
  d_b <- 5L
  n_y <- 3L
  d_z <- 4L
  d_x <- 0L
  K <- n_y * d_m

  b_ks <- seq(log(1e-2), log(10), length.out = d_b)
  m_ks <- seq(0        , 10     , length.out = d_m)
  g_ks <- seq(0        , 10     , length.out = d_g)

  omega <- c(-1.5, -1.35, -1.9, -4.8, -0.17)
  B <- structure(c(-0.86, -0.28, -0.58, 0.04, 0.18, -0.12, 0.99, -0.23,
                   -0.48, -0.24, 0.33, 0.06, 0.3, -0.57, -0.25),
                 .Dim = c(5L, 3L))
  Psi <- structure(c(2.54, -0.64, 0.95, -0.35, -0.73, 0.79, -0.64, 1.91,
                     -0.66, 0.26, -0.78, 1.08, 0.95, -0.66, 4.13, 1.99, 0.09, -0.94,
                     -0.35, 0.26, 1.99, 3.86, 0.76, -1.34, -0.73, -0.78, 0.09, 0.76,
                     1.99, -1, 0.79, 1.08, -0.94, -1.34, -1, 3.06),
                   .Dim = c(6L, 6L))
  sig <- structure(c(0.02, 0, 0, 0, 0.15, 0, 0, 0, 0.03),
                   .Dim = c(3L, 3L))
  delta <- c(-0.1, -0.08, 0.03, -0.21)
  alpha <- c(0.23, -0.07, -0.15)
  gamma <- NULL

  n_obs <- 2000L
})