We draw upon a series of articles in which a method based on pseudovalues is proposed for direct regression modeling of the survival function, the restricted mean, and the cumulative incidence function in competing risks with right-censored data. The models, once the pseudovalues have been computed, can be fit using standard generalized estimating equation software. Here we present Stata procedures for computing these pseudo-observations. An example from a bone marrow transplantation study is used to illustrate the method.
AndersenP. K., HansenM. G., and KleinJ. P.2004. Regression analysis of restricted mean survival time based on pseudo-observations. Lifetime Data Analysis10: 335–350.
2.
AndersenP. K., and KleinJ. P.2007. Regression analysis for multistate models based on a pseudovalue approach, with applications to bone marrow transplantation studies. Scandinavian Journal of Statistics34: 3–16.
3.
AndersenP. K., KleinJ. P., and RosthøjS.2003. Generalised linear models for correlated pseudo-observations, with applications to multistate models. Biometrika90: 15–27.
4.
AndersenP. K., and PermeM. P.2010. Pseudo-observations in survival analysis. Statistical Methods in Medical Research19: 71–99.
5.
CopelanE. A., BiggsJ. C., ThompsonJ. M., CrilleyP., SzerJ., KleinJ. P., KapoorN., AvalosB. R., CunninghamI., AtkinsonK., DownsK., HarmonG. S., DalyM. B., BrodskyI., BulovaS. I., and TutschkaP. J.1991. Treatment for acute myelocytic leukemia with allogeneic bone marrow transplantation following preparation with BuCy2. Blood78: 838–843.
6.
FineJ. P., and GrayR. J.1999. A proportional hazards model for the subdistribution of a competing risk. Journal of the American Statistical Association94: 496–509.
7.
GrawF., GerdsT. A., and SchumacherM.2009. On pseudovalues for regression analysis in competing risks models. Lifetime Data Analysis15: 241–255.
8.
KaplanE. L., and MeierP.1958. Nonparametric estimation from incomplete observations. Journal of the American Statistical Association53: 457–481.
9.
KleinJ. P.2006. Modeling competing risks in cancer studies. Statistics in Medicine25: 1015–1034.
10.
KleinJ. P., and AndersenP. K.2005. Regression modeling of competing risks data based on pseudovalues of the cumulative incidence function. Biometrics61: 223–229.
11.
KleinJ. P., GersterM., AndersenP. K., TarimaS., and PermeM. P.2008. SAS and R functions to compute pseudovalues for censored data regression. Computer Methods and Programs in Biomedicine89: 289–300.
12.
KleinJ. P., LoganB., HarhoffM., and AndersenP. K.2007. Analyzing survival curves at a fixed point in time. Statistics in Medicine26: 4505–4519.
13.
KleinJ. P., and MoeschbergerM. L.2003. Survival Analysis: Techniques for Censored and Truncated Data. 2nd ed. New York: Springer.
14.
LiangK.-Y., and ZegerS. L.1986. Longitudinal data analysis using generalized linear models. Biometrika73: 13–22.
15.
PermeM. P., and AndersenP. K.2008. Checking hazard regression models using pseudo-observations. Statistics in Medicine27: 5309–5328.
