loo_ebmstate {ebmstate}

Leave-one-out estimation

Description

This function computes leave-one-out estimation of regression coefficients, cumulative hazard functions, and transition probability functions.

Usage

loo_ebmstate(
  mstate_data,
  mstate_data_expanded,
  which_group,
  patient_IDs,
  initial_state,
  tmat,
  time_model,
  backup_file = NULL,
  input_file = NULL,
  coxrfx_args = list(),
  msfit_args = NULL,
  probtrans_args = NULL
)

Arguments

mstate_data

Data in 'long format'.

mstate_data_expanded

Data in 'long format', possibly with 'expanded' covariates (as obtained by running mstate::expand.covs).

which_group

A character vector with the same meaning as the 'groups' argument of the function CoxRFX but named (with the covariate names).

patient_IDs

The IDs of the patients whose cumulative hazards and transition probabilities one wishes to estimate.

initial_state

The initial state for which transition probability estimates should be computed

tmat

Transition matrix for the multi-state model, as obtained by running mstate::transMat

time_model

The model of time-dependency: either 'clockforward' or 'clockreset'.

backup_file

Path to file. Objects generated while the present function is running are stored in this file. This avoids losing all estimates if and when the algorithm breaks down. See argument input_file.

input_file

Path to backup_file (see argument backup_file). If this argument is given, all other arguments should be NULL.

coxrfx_args

Named list with arguments to the CoxRFX function other than Z,surv and groups.

msfit_args

Named list with arguments to the msfit_generic.coxrfx function other than object,newdata and trans.

probtrans_args

Named list with arguments to the probtrans_ebmstate function other than initia_state,cumhaz and model.

Details

In a given bootstrap sample there might not be enough information to generate estimates for all coefficients. If a covariate has little or no variation in a given bootstrap sample, no estimate of its coefficient will be computed. The present function will keep taking bootstrap samples until every coefficient has been estimated at least min_nr_samples times.

Value

A list with: 95% bootstrap intervals for each regression coefficient and for transition probabilities; bootstrap samples of regression coefficients, cumulative hazards and transition probabilities.

Author(s)

Rui Costa