inst/test-data/no-B.R

In SimSurvNMarker: Simulate Survival Time and Markers

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

  b_ks <- seq(log(1e-1), 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(-0.79, -1.25, -4.46, 0.14)
  B <- NULL
  Psi <- structure(c(2.58, -0.59, 0.87, -0.61, 0.13, 0.52, -0.59, 1, -0.3,
                     0.27, 0.02, 0.2, 0.87, -0.3, 1.12, 0.03, -0.36, 0, -0.61, 0.27,
                     0.03, 2.19, -0.2, 0.06, 0.13, 0.02, -0.36, -0.2, 0.85, 0.32,
                     0.52, 0.2, 0, 0.06, 0.32, 1.32), .Dim = c(6L, 6L))
  sig <- structure(c(0.14, 0, 0, 0, 0.14, 0, 0, 0, 0.02),
                   .Dim = c(3L, 3L))
  delta <- c(-0.05, 0.05, -0.15, -0.16)
  alpha <- c(0.2, -0.25, 0.1)
  gamma <- structure(c(-0.23, -0.15, 0.19, -0.16, -0.22, -0.17, 0.2, -0.21,
                       -0.1, 0.09, 0.13, 0.14), .Dim = 4:3)

  n_obs <- 2000L
})