TP.idm-package: Estimation of Transition Probabilities for the Illness-Death...

Description Details Author(s) References Examples

Description

Estimation of transition probabilities for the illness-death model. Both the Aalen-Johansen estimator for a Markov model and a novel non-Markovian estimator by de Uña-Álvarez and Meira-Machado (2015) <doi:10.1111/biom.12288>, see also Balboa and de Uña-Álvarez (2018) <doi:10.18637/jss.v083.i10>, are included.

Details

The DESCRIPTION file: This package was not yet installed at build time.

Index: This package was not yet installed at build time.

This package incorporates the function TPidm which can be used to compute a nonparametric transition probability matrix for the illness-death model. For a complete list of functions use library(help=TP.idm).

Author(s)

Vanesa Balboa-Barreiro, Jacobo de Una-Alvarez <[email protected]> and Luis Meira-Machado <[email protected]>

Maintainer: Vanesa Balboa-Barreiro <[email protected]>

References

Andersen P. K., Borgan O., Gill R. D. and Keiding N. (1993). Statistical models based on counting processes. Springer Series in Statistics. New York, NY: Springer.

Datta S. and Satten G. A. (2001). Validity of the Aalen-Johansen estimators of stage occupation probabilities and Nelson-Aalen estimators of integrated transition hazards for non-Markov models. Statistics and Probability Letters 55(4), 403–411.

de Uña-Álvarez J. and Meira-Machado L. (2015). Nonparametric estimation of transition probabilities in a non-Markov illness-death model: a comparative study. Biometrics 71, 364–375.

Meira-Machado L. F., de Uña-Álvarez J. and Cadarso-Suárez C. (2006). Nonparametric estimation of transition probabilities in a non-Markov illness-death model. Lifetime Data Annal 12(3), 325–344.

Balboa-Barreiro V., de Uña-Álvarez J. (2018). Estimation of Transition Probabilities for the Illness-Death Model: Package TP.idm. Journal of Statistical Software 83(10), 1–19.

Examples

1
2
3
4
5
6
7
data(colonTP)

# create a TPidm object with s = 0, t = "last" (default),
# CI = TRUE (default) and method = "AJ":
aj0ci <- TPidm(colonTP, s = 0, method = "AJ")
summary(aj0ci) # summarizing the results
plot(aj0ci) # plotting all occupation probabilities

Example output

Call:
TPidm(data = colonTP, s = 0, method = "AJ")

Parameters: 

Time s:
[1] 0

Time t:
[1] 2910

Method:
[1] "AJ"

CI:
[1] TRUE

CI transformation:
[1] "linear"

Possible transitions:
[1] "1 1" "1 2" "1 3"

Occupation probabilities at time t: 
 

 transition      probs      lower     upper     variance
        1 1 0.42892004 0.39697971 0.4608604 2.655721e-04
        1 2 0.03249031 0.02117202 0.0438086 3.334767e-05
        1 3 0.53858966 0.50649029 0.5706890 2.682235e-04

TP.idm documentation built on April 19, 2018, 1:04 a.m.