Run stan mixture cure joint model

Use pre-compiled C code.

bmcm_joint_stan(
  input_data,
  model_os = "exp",
  model_pfs = "exp",
  tx_name = "IPILIMUMAB",
  iter = 2000,
  warmup = 1000,
  thin = 10,
  chains = 2,
  mean_cf = NA,
  var_cf = NA,
  centre_age = TRUE,
  ...
)

Arguments

input_data	Dataframe with all variables
model_os	Distribution name; exp, weibull, gompertz
model_pfs	Distribution name; exp, weibull, gompertz
tx_name	Treatment name; IPILIMUMAB, NIVOLUMAB, NIVOLUMAB+IPILIMUMAB
iter	int
warmup	Number of iterations burn-in; int
chains	Number of chains; int
mean_cf	Mean cure fraction; default to U0,1
var_cf	Variance of cure fraction
centre_age	Logical to centre regression covariate
...	Additional arguments