Description Usage Arguments Value Examples
The function VarianceEM performs variance estimation for the Multistate Cure Model when the model is fit using the EM Algorithm. Variances are estimated by fitting the multistate cure model (via EM) to bootstrap samples of the data.
1 2 | VarianceEM(fit, iternum, bootnum, datWIDE, Cov, ASSUME, TransCov, BASELINE,
PENALTY = "None")
|
fit |
Multistate cure model fit from MultiCure |
iternum |
Number of iterations of the EM algorithm to apply to each bootstrap sample |
bootnum |
Number of bootstrap samples |
datWIDE |
A data frame with the following columns:
|
Cov |
A matrix containing the covariates |
ASSUME |
This variables indicates what equality assumptions we are making regarding the 24 and 14 transitions. The possible options are:
|
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) |
OUT a matrix containing the following:
Estimate an estimate of the multistate cure model parameter
Variance an estimate of variance of the the multistate cure model parameter
1 2 3 4 5 6 7 | attach(SimulateMultiCure(type = "NoMissingness"))
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")
OUT = VarianceEM(fit,iternum=20, bootnum=50, datWIDE, Cov, ASSUME = "SameHazard", TransCov, BASELINE = "weib")
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.