Histogram of survival times
# S3 method for modSim
hist(x, ...)
output of modelSim function (must be of type modSim)
supplementary parameters
hist x
library(survMS)
### Survival data simulated from AH model
res_paramLN = get_param_ln(var=170000, mu=2325)
listAHSim_n500_p1000 <- modelSim(model = "AH", matDistr = "unif", matParam = c(-1,1), n = 500,
p = 100, pnonull = 100, betaDistr = 1.5, hazDistr = "log-normal",
hazParams = c(res_paramLN$a*4, res_paramLN$lambda),
Phi = 0, seed = 1, d = 0)
#> Warning: Options "non-linearity with interactions" not available
### Histogram of survival times
hist(listAHSim_n500_p1000)