VarianceMCEM_NOBOOT: VarianceMCEM_NOBOOT

Description Usage Arguments Details Value Examples

View source: R/VarianceMCEM_NOBOOT.R

Description

The function VarianceMCEM_NOBOOT performs variance estimation for the Multistate Cure Model when the model is fit using a Monte Carlo EM Algorithm. For each imputed dataset, this function estimates parameter variances by fitting the multistate cure model the imputed dataset. Corresponding standard errors come are extracted from the corresponding logistic and proportional hazards model fits. Then, Rubin's rules are used to obtain a single set of parameter estimates and standard errors across imputed datasets. Note: This function has not yet been implemented for PENALTY values other than 'None'. Important!!!: The function ProperDraws_MC should be used to obtain the imputed datasets to be used in this function. This ensures that the imputed datasets are (roughly) proper imputations and that Rubin's rules can then be applied.

Usage

1
2
VarianceMCEM_NOBOOT(fit, datWIDE, CovImp, GImp, YRImp, deltaRImp, ASSUME,
  TransCov, BASELINE, PENALTY = "None")

Arguments

fit

Multistate cure model fit from MultiCure

datWIDE

A data frame with the following columns:

  • Y_R, the recurrence event/censoring time

  • delta_R, the recurrence event/censoring indicator

  • Y_D, the death event/censoring time

  • delta_D, the death event/censoring indicator

  • G, the cure status variable. This takes value 1 for known non-cured, 0 for "known" cured and NA for unknown cure status

CovImp

A list with IMPNUM elements containing the imputations of Cov output from the function ProperDraws_MC

GImp

A matrix with IMPNUM elements containing the imputations of G output from the function ProperDraws_MC

YRImp

A matrix with IMPNUM elements containing the imputations of Y_R output from the function ProperDraws_MC

deltaRImp

A matrix with IMPNUM elements containing the imputations of delta_R output from the function ProperDraws_MC

ASSUME

This variables indicates what equality assumptions we are making regarding the 24 and 14 transitions. The possible options are:

  • 'SameHazard': Lambda_14(t) = Lambda_24(t)

  • 'AllSeparate': No restrictions on Lambda_14(t) and Lambda_24(t)

  • 'ProportionalHazard': Lambda_14(t) = Lambda_24(t) exp(Beta0)

  • 'SameBaseHaz': Lambda^0_14(t) = Lambda^0_24(t), No restrictions on beta_14 and beta_24

TransCov

a list with elements: Trans13, Trans24, Trans14, Trans34, PNonCure. Each list element is a vector containing the names of the variables in Cov to be used in the model for the corresponding transition. 13 is NonCured -> Recurrence, 24 is Cured -> Death, 14 is NonCured -> Death, 34 is Recurrence -> Death. PNonCure contains the names of the covariates for the logistic regression for P(NonCure).

BASELINE

This variable indicates the assumptions about the baseline hazard form. This can take values 'weib' and 'cox'

PENALTY

This variable indicates whether we are using any variable selection in the model fitting. Right now, the options are 'None' (no variable selection), 'Ridge' (ridge regression for all covariates in all models) and 'Lasso' (lasso for all covariates in all models, only implemented for Cox baseline hazards)

bootnum

Number of bootstrap samples used for each imputed dataset

Details

This function provides parameter estimates and estimated variances. The parameter estimates are obtained using Rubin's rules, but an alternative estimate of the multistate cure model parameter can be obtained by averaging the parameter estimates from the last few iterations of the model fitting algorithm. In our experience, we found that the approach that averages across the last few iterations (rather than estimated using Rubin's rules) provides a better estimate of the parameter of interest.

Value

OUT a matrix containing the following:

Examples

1
2
3
4
5
6
7
8
attach(SimulateMultiCure(type = "UnequalCensoring"))
Cov = data.frame(X1,X2)
VARS = names(Cov)
TransCov = list(Trans13 = VARS, Trans24 = VARS, Trans14 = VARS, Trans34 = VARS, PNonCure = VARS)
datWIDE = data.frame( Y_R, Y_D, delta_R , delta_D, G)
fit = MultiCure(iternum = 100, datWIDE, Cov, ASSUME = "SameHazard", TransCov = TransCov, BASELINE = "weib", IMPNUM = 10) ### Note: This will take a moment
Proper = ProperDraws_MC(datWIDE,Cov, CovImp = fit[[9]], GImp = fit[[10]], YRImp = fit[[11]], deltaRImp = fit[[12]], ASSUME = "SameHazard", TransCov = TransCov, BASELINE = "weib") ### Note: This will take a moment
OUT = VarianceMCEM_NOBOOT(fit,datWIDE, CovImp = Proper[[1]], GImp = Proper[[2]], YRImp = Proper[[3]], deltaRImp = Proper[[4]],  ASSUME = "SameHazard", TransCov, BASELINE = "weib")

lbeesleyBIOSTAT/MultiCure documentation built on April 18, 2018, 11:08 p.m.