MI.ci(k, m, data, status, trans, cens.code, conf.int = F, alpha = 0.05)
est
A data frame with estimates…
Other objects
This function uses a multiple imputation approach to estimate a cumulative incidence function for interval censored competing
risks data.
Estimates are computed using Rubin's rules (Rubin (1987)). The cumulative incidence is computed as the mean of
cumulative incidences over imputations. The variance is computed at each point by combining the within imputation variance and the
between imputation variance augmented by an inflation factor to take into account the finite number of imputations.
At each iteration, the cumulative incidence is updated and multiple imputation is performed using the updated estimate.
If conf.int
is required, the log-log transformation is used to compute the lower confidence interval.
Print and plot methods are available to handle results.
The data
must contain at last three columns: left
, right
and status
. For interval censored data, the
left
and right
columns indicates lower and upper bounds of intervals, respectively. Inf
in the
right
column stands for right censored observations. When an observation is right censored, the status
column must
contain the censor indicator specified by cens.code
. The transition of interest must be specified by the trans
parameter.
PAN, Wei. A Multiple Imputation Approach to Cox Regression with Interval-Censored Data. Biometrics, 2000, vol. 56, no 1, p. 199-203.
Rubin, D. B. (1987). Multiple imputation for nonresponse in surveys.
Schenker, N. and Welsh, A. (1988). Asymptotic results for multiple imputation. The Annals of Statistics pages 1550-1566.
Tanner, M. A. and Wong, W. H. (1987). An application of imputation to an estimation problem in grouped lifetime analysis. Technometrics 29, 23-32.
Wei, G. C., & Tanner, M. A. (1991). Applications of multiple imputation to the analysis of censored regression data. Biometrics, 47(4), 1297-1309.
res <- MI.ci(k = 5, m = 5, status = 'status', trans = 1 , data = ICCRD,
conf.int = TRUE, cens.code = 0 , alpha = 0.05)
res
print(res)
plot(res)
Run the code above in your browser using DataLab