Example joint longitudinal and time-to-event model
A model for use in the rstanarm examples related to stan_jm.
Calling example("example_jm") will run the model in the
Examples section, below, and the resulting stanmvreg object will then be
available in the global environment. The chains and iter
arguments are specified to make this example be small in size. In practice,
we recommend that they be left unspecified in order to use the default
values or increased if there are convergence problems. The cores
argument is optional and on a multicore system, the user may well want
to set that equal to the number of chains being executed.
# set.seed(123)
if (.Platform$OS.type != "windows" || .Platform$r_arch !="i386")
example_jm <-
stan_jm(formulaLong = logBili ~ year + (1 | id),
dataLong = pbcLong[1:101,],
formulaEvent = survival::Surv(futimeYears, death) ~ sex + trt,
dataEvent = pbcSurv[1:15,],
time_var = "year",
# this next line is only to keep the example small in size!
chains = 1, seed = 12345, iter = 100, refresh = 0)Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.