R/simfitRI_DPBrier.R

In JMH: Joint Model of Heterogeneous Repeated Measures and Survival Data

simfitRI_DPBrier <- function(sim = 10, N = 500, seed = 100, increment = 0.7, beta = c(5, 1.5, 2, 1, 2),
                        tau = c(0.5, 0.5, -0.2, 0.2, 0.05),
                        gamma1 = c(1, 0.5, 0.5),
                        gamma2 = c(-0.5, 0.5, 0.25),
                        alpha1 = 1,
                        alpha2 = -1,
                        vee1 = 0.5,
                        vee2 = -0.5,
                        lambda1 = 0.05,
                        lambda2 = 0.1,
                        CL = 4,
                        CU = 8,
                        covbw = matrix(c(0.5, 0.25, 0.25, 0.5), nrow = 2, ncol = 2),
                        maxiter = 1000, quadpoint = 10, 
                        landmark.time = 2, horizon.time = seq(3, 6, by = 1), n.cv = 3,
                        method = "GH", ncores = 10) {
  
  
  ParaMatrixRaw <- parallel::mclapply(1:sim, bootsfitRI_Brier,
                                      seed = seed, N = N, increment = increment, 
                                      beta = beta, tau = tau, 
                                      gamma1 = gamma1, gamma2 = gamma2,
                                      alpha1 = alpha1, alpha2 = alpha2, 
                                      vee1 = vee1, vee2 = vee2, 
                                      lambda1 = lambda1, lambda2 = lambda2, 
                                      CL = CL, CU = CU, covbw = covbw, 
                                      quadpoint = quadpoint, maxiter = maxiter, 
                                      n.cv = n.cv, landmark.time = landmark.time, 
                                      horizon.time = horizon.time, method = method, 
                                      mc.cores = ncores)
  
  return(ParaMatrixRaw)
  
}