# NOT RUN {
# }
# NOT RUN {
###---Data generation---###
data.sim <-jointSurrSimul(n.obs=400, n.trial = 20,cens.adm=549,
alpha = 1.5, theta = 3.5, gamma = 2.5, zeta = 1,
sigma.s = 0.7, sigma.t = 0.7, cor = 0.8, betas = -1.25,
betat = -1.25, full.data = 0, random.generator = 1,
seed = 0, nb.reject.data = 0)
###---Estimation---###
joint.surrogate <- jointSurroPenal(data = data.sim, nb.mc = 300,
nb.gh = 20, indicator.alpha = 1, n.knots = 6)
print(joint.surrogate)
# or
joint.surrogate
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab