Home | About | Journals | Submit | Contact Us | Français |

**|**HHS Author Manuscripts**|**PMC3069508

Formats

Article sections

- Abstract
- 1. Introduction
- 2. REGRESSION MODELS
- 3. APPLICATION TO DATA FROM THE WOMEN’S INTERAGENCY HIV STUDY
- 4. DISCUSSION
- References

Authors

Related links

Stat Med. Author manuscript; available in PMC 2012 March 15.

Published in final edited form as:

PMCID: PMC3069508

NIHMSID: NIHMS247529

The publisher's final edited version of this article is available at Stat Med

In the analysis of survival data, there are often competing events that preclude an event of interest from occurring. Regression analysis with competing risks is typically undertaken using a cause-specific proportional hazards model. However, modern alternative methods exist for the analysis of the subdistribution hazard with a corresponding subdistribution proportional hazards model. In this paper, we introduce a flexible parametric mixture model as a unifying method to obtain estimates of the cause-specific and subdistribution hazards and hazard ratio functions. We describe how these estimates can be summarized over time to give a single number that is comparable to the hazard ratio that is obtained from a corresponding cause-specific or subdistribution proportional hazards model. An application to the Women’s Interagency HIV Study is provided to investigate injection drug use and the time to either the initiation of effective antiretroviral therapy, or clinical disease progression as a competing event.

With survival data, a competing risk situation occurs when there exists an event which precludes the event of interest from occurring. The canonical example is in the analysis of specific causes of death. Death due to cancer removes the possibility of the individual experiencing death due to cardiovascular disease, for example (for reviews of competing risks, see [1]). The components of survival analysis (e.g. risk sets, hazard functions, and hazard ratios) and modeling strategies are similar but altered in construction when competing risks exist [2]. Two semi-parametric methods have become prominent approaches for competing risk analyses, a cause-specific hazard [3, 4] and a subdistribution hazard approach [6]. The cause-specific hazard approach and accompanying cause-specific hazards ratio (_{cs}HR) was introduced in 1978 [3, 4] and only more recently in pioneering work by Fine and Gray [6], has the subdistribution approach for estimating the subdistribution HR (_{sd}HR) been introduced. The _{sd}HR was motivated by the fact that the cause-specific HR (_{cs}HR) may not reflect an association on the cumulative incidence function (CIF) [6]. Recently Beyersmann, et al. [7] illustrated a case in which the _{cs}HR suggested a significant association, but the _{sd}HR had a non-significant association in the opposite direction suggesting that the covariate does not affect the CIF. Beyersmann and Schumacher [8] suggested that it is for this reason that Fine and Gray’s method [6] has been gaining attention in the medical literature.

An alternative to the semi-parametric methods for modeling data with competing risks is to use a parametric mixture model. Initially suggested by Cox [11], Larson and Dinse [12] expanded on the idea by decomposing the CIF into a distribution conditioned on the outcome and the probability of that outcome. Further extensions were made by Maller and Zhou [13]. Previous applications of the mixture model have focused almost exclusively on obtaining an estimate of the CIF. In previous work, the mixture model was suggested as a unifying approach between the cause-specific and subdistribution proportional hazards models, but details were not given [2]. The purpose of this paper is to give details on obtaining the _{cs}HR(*t*) and _{sd}HR(*t*) functions over follow-up time from a parametric mixture model that may be summarized into estimates corresponding to their proportional hazards model counterparts. While we provide an approach to summarize the _{cs}HR(*t*) and _{sd}HR(*t*) functions over follow-up from a parametric mixture model, it is recognized that a summary measure of the hazard ratio may have two shortcomings whether it comes from a parametric or semi-parametric model [14]. First, if the hazard ratio changes over time, then the summary association (i.e. average HR(*t*)) that an observational study will find may be dependent on the length of follow-up. The second issue is that the HR(*t*) has a built in selection bias in that individuals at time *t* are those remaining at risk for the outcome which may consist of those who are at less prone to develop the event [14]. Therefore, an advantage of the proposed parametric mixture model is that either the HR(t) or the summary HR for either the cause-specific or subdistribution hazards may be presented, whereas the semi-parametric proportional hazards models must either present the summary HR or be complicated by time-exposure interactions. In Section 2, we review the semi-parametric approaches and detail how the parametric mixture model can be used for hazards-ratio estimation under competing risks. In Section 3, we apply our approach to data from the Women’s Interagency HIV Study consisting of a subsample of 1,164 women followed over 10 years. Finally in Section 4, we discuss the strengths and limitations of estimating either the _{cs}HR or _{sd}HR from a parametric mixture model.

Let the observed survival data be represented by (*T _{i}*,

Proportional hazards models using either the cause-specific hazard or the subdistribution hazard have been proposed [3, 6]. The cause-specific hazard is defined as [3]:

(1)

where **X** is a vector of covariates and *S*(*t*) is the net survival function such that in which *h*(*t*) is the net hazard function for having any of the competing events. This can be used to construct a cause-specific proportional hazards model [3]:

(2)

where *h*_{0j}(*t*) is an arbitrary baseline cause-specific hazard function and *β _{j}* is the vector of unknown coe cents for covariates

Alternatively, the subdistribution hazard may be used directly to construct a different proportional hazards model. The subdistribution hazard is defined as [5, 6]:

(3)

where is the subsurvivorship function. The corresponding proportional hazards model is:

(4)

where *λ*_{0j}(*t*) is an arbitrary baseline subdistribution hazard function and * _{j}* is the vector of unknown coefficients for covariates

Larson and Dinse [12] describe an alternative competing risks model represented by a parametric mixture model rather than the semi-parametric proportional hazards model. The CIF corresponding to the *j*th type of failure is decomposed as:

with the corresponding subdensity function defined as:

where **X**_{j} **X** and **X**_{ω} **X**. The distribution functions *F _{j}*(

(5)

where **X**_{i1} and **X**_{i2} are the *i*th individual’s corresponding set of covariates being modeled for the events *J* = 1 and *J* = 2 respectively; *t _{i}* is the

such that individuals are right censored if *δ _{i}* =

Since , , and the net survivorship function (*S*(*t*) = *π*(*S*_{1}(*t*)) + (1 − *π*)*S*_{2}(*t*)) are determined by the parameter estimates, it is easy to see that the cause-specific hazard functions and the subdistribution hazard functions can also be determined. These building blocks allow for the estimation of the ratio of cause-specific hazards, as a function of time *t*:

where is the corresponding subdensity function for the arbitrary baseline hazard function and *S*_{0j}(*t*) is the net survival function for individuals with the covariate values that correspond to the reference population that defines the baseline arbitrary hazard function. Similarly the subdistribution hazard ratio curves can be estimated as:

where and are the corresponding subdensity function and CIF for the arbitrary baseline hazard function, respectively.

The likelihood function in (5) can be modified to allow for left truncation (Appendix). This would allow estimation of the HR(*t*) while avoiding bias due to late entry into the study or for incorporating time-varying covarites into the survival models [15]. While it is straight forward to incorporate left truncation into the cause-specific proportional hazards model, it is not easily done for the subdistribution proportional hazards model [16]. Additionally uncertainty in exact timing of the event could be incorporated through interval censoring (Appendix), which is not easily accounted for in the semi-parametric proportional hazards models.

To facilitate appropriate inferences, the parametric distributions to be included in the mixture model need to accurately reflect the data. To estimate *f*_{1}(*t _{i}*) and

The above approach provides an evaluation of how exposures impact the hazard rates of different competing events over time. In many settings, it may be desirable to summarize these hazards over follow-up and provide a single summary estimate comparable to the HR typically obtained from a proportional hazards model. This is often the case in the medical setting in which clinicians want to quickly disseminate information regarding the relationship between an exposure and outcome. A single number is much easier to convey compared to either two hazard functions or a HR curve in which each is a function of time. The summary HR from a proportional hazards model where the proportionality assumption is unmet essentially weights the ratio of the time-dependent hazards at each failure time-point equivalently assuming few or no tied events [25, 26]. Such weighting may not be the most appropriate in summarizing the time-dependent hazard ratios [26, 27, 28, 29]. Briefly, Kalbfleisch and Prentice [27] introduced a summary measure of the hazard ratio when proportionality was not met by allowing the hazards to be weighted by some survivor function, for (*t* > 0) in which *S*_{0}(*t*) and *S*_{1}(*t*) are the survivor functions for the two sample problem and *α* > 0 such that *α* = 0.5 would be the geometric average of the two survivor functions and *α* < .5 or α > .5 would emphasize earlier or later periods, respectively [27]. Additional work in this area has been done providing potentially other metrics, e.g. weighted by number at risk providing a Wilcoxon-type estimator [26, 28, 29]. However, Grambauer et al. have shown that the subdistribution hazard ratio when mis-specified is still useful as it is the least false parameter which may be interpreted as the time-averaged effect on the cumulative event probabilities [10]. Nevertheless, we present here a simple outline in which all points are equivalently weighted. However, our methods may be appropriately adapted to include such alternative weights. We outline two methods. The first estimates the hazard curves for an “average” individual with and without exposure. The ratio of these cuves are then determined at each of the *w* unique failure time points.

• “Average” individual approach (Method 1)

- Create a record with the mean covariate level for each covariate such that the record has a covariate vector of means .
- Let
*w*=*w*_{1}, … ,*w*correspond to the ordered unique failure times,_{K}*T*such that*w*_{1}= min(*T*_{1}, … ,*T*) and_{N}*w*= max(_{K}*T*_{1}, … ,*T*)._{N} - Duplicate the mean covariate record and assign the exposure variable to be “unexposed” to the original record and the exposure variable as “exposed” for the duplicate record. Should the exposure of interest be continuous then the difference between exposed and unexposed values should be 1 unit of scientific interest centered around the study population mean.
- Duplicate the pair of records
*K*times giving the*k*th replicate the*w*failure time._{k} - Estimate either the cause-specific or subdistribution hazard for each of the 2
*W*records and obtain HR(*t*=*w*) =_{k}*g*_{1}(*t*)/*g*_{0}(*t*) where*g*_{1}(*t*) and*g*_{0}(*t*) correspond to the exposed and unexposed hazard functions, respectively. - The mean of these time-specific HRs (e.g. ) corresponds to the
_{cs}HR or_{sd}HR obtained from the proportional hazards models depending on whether*g*(*t*) is the cause-specific or subdistribution hazard function, respectively.

This approach is easily accomplished in SAS by creating the average unexposed and exposed records and appending this to the original data file. These additional records can be excluded from the fitting process by utilizing the replicate statement with a weight variable to identify the original records from the created records. Additional programming statements to the NLMIXED procedure are then added to estimate the both the cause-specific and subdistribution hazards and the predict statement is used to output each of the hazards.

The second method entails generating two hazard curves for each individual (one for the observed exposure status, and one for the unobserved exposure status).

• Empirically weighted approach (Method 2)

- Let
*w*=*w*_{1}, … ,*w*correspond to the ordered unique failure times,_{K}*T*such that*w*_{1}= min(*T*_{1}, … ,*T*) and_{N}*w*= max(_{K}*T*_{1}, … ,*T*)._{N} - For each of the
*N*individuals, create a duplicate record that matches the original on all covariates with the exception of the variable of interest. Instead, the duplicate record is given the unobserved (i.e., opposite) exposure value. - For each
*N*pairs of records, replicate the pair*K*times giving the*k*th pair of records the*w*failure time. Therefore the final database will have 2 *_{k}*N***K*records. - Estimate either the cause-specific or subdistribution hazards at the corresponding failure time for each of the records by using the estimated model parameters and the cause-specific hazard function and subdistribution hazard function . Thus the HR
_{i}(*t*=*w*) =*g*_{1i}(*t*)/*g*_{0i}(*t*) where*g*_{1i}(*t*) and*g*_{0i}(*t*) correspond to the exposed and unexposed hazard functions, respectively, for the*i*th individual at the*w*failure time. - Then the corresponds to either the
_{cs}HR and_{sd}HR depending on whether*g*(*t*) is the cause-specific or subdistribution hazard function, respectively.

We applied these methods to data from the Women’s Interagency HIV Study (henceforth, the Study) to examine whether women with a history of injection drug use were more or less likely to have one of two competing events: either initiating highly active antiretroviral therapy (henceforth, therapy) or the onset of clinical disease progression as defined by the occurrence of either an AIDS defining illness or death [2]. Briefly, the Study was established in August of 1993 to comprehensively examine the impact of HIV on women at six US sites [30]. For this application, we are interested in whether individuals with a history of injection drug use were less likely to initiate therapy prior to AIDS or death as compared to those without a history of injection drug use. Prior studies have shown that HIV-infected individuals with past injection drug use are less likely to initiate therapy [31, 32, 33, 34] and are more likely to die in the era of effective treatment for HIV [34, 35, 36]. Yet the comparison of treatment initiation prior to disease progression (defined by AIDS or death) has not been well explored in a competing risk framework.

The study sample consisted of 1,164 women that were alive and free of clinical AIDS as of December 6, 1995, when the first protease inhibitor (Saquinavir mesylate) was approved by the US Federal Drug Administration. This date marks the first time-point in which individuals could be placed on a highly active antiretoviral drug regimen. Women were followed until the first of: initiation of therapy, an AIDS diagnosis, death, or administrative censoring defined as September 28, 2006. Covariates included history of injection drug use upon study enrollment, whether an individual was African-American, age on December 6 1995, and CD4 nadir prior to December 6 1995.

The study sample consisted of 439 women with and 725 without a history of injection drug use. Overall, 679 (58%) women initiated therapy prior to either AIDS or death; 359 (31%) had the competing event and either died or had an AIDS diagnosis before initiation of effective therapy. The remaining 11% were censored (5% loss to follow-up, 5% administratively censored). Reflecting the demographics of HIV in women and in the cohort, the majority of the study population was African-American (58%) and had a median age of 36 (interquartile range 31-41) as of December 6, 1995. A greater proportion of injection drug users tended to be African-American and of older age (p= 0.017 and *p*< 0.001, respectively).

The maximum likelihood estimates for the mixture model are shown in Table 1. A lognormal distribution was used to model initiation of therapy (*f*_{1}) as a sufficiently small enough gradient could not be achieved for the shape parameter when modeling a generalized gamma distribution. However, the shape parameter indicated that the distribution was tending towards a lognormal distribution which is a limiting case for the generalized gamma distribution [20, 37]. A generalized gamma distribution was used to model the incidence of AIDS/Death prior to therapy (*f*_{2}). Under this model, an African-American woman with a history of injection drug use that is 36 years of age with a CD4 nadir of 349 cells/mm^{3} has a probability of 0.46 for ever initiating therapy (*π _{i}* =

The cumulative incidence function (panels a and b), cause-specific hazards (panels c and d, dotted and dash-dot-dash curves), and subdistribution hazards (panels e and f, dotted and dash–dot–dash curves) from the multivariate mixture model **...**

The Cox proportional hazards model (and similarly the subdistribution proportional hazards model) is often used in research despite potential violation of the proportional hazards assumption. While the subdistribution proportional hazards model still provides the least false parameter when mis-specified [19, 10], the summary hazard ratios from any proportional hazards model still has shortcomings such as the dependency on total length of follow-up and the potential for selection bias [14]. Despite concerns about the appropriateness of the proportional hazards assumption (global test for cause-specific proportional hazards assumption, therapy initiation: p< 0.001; AIDS/Death: p=0.29, subdistribution proportional hazards assumption, therapy initiation: p< 0.001; **AIDS/Death: p**< 0.001), we constructed both cause-specific and subdistribution hazard ratios (Table 2). Both the _{cs}HR and _{sd} HR their 95% confidence intervals (obtained from bootstrap methods with 1000 replicates) were similar between the summarized mixture estimates and proportional hazards models. To obtain the bootstrap confidence intervals individuals were resampled all 1164 women with replacement and the model was re-estimated based upon each sample. The 0.025 and 0.975 quantiles of the HR estimates were then obtained to provide the confidence intervals. These results show that the parametric mixture model allows for the estimation of both the summarized cause-specific and subdistribution HR that is equivalent to the semi-parametric proportional hazards models. However, an advantage of the parametric mixture model over the semi-parametric approaches is the ablity to handle left truncation (subdistribution proportional hazards model [16]) and interval censoring (both proportional hazards models) easily. Furthermore, both semi-parametric models may fail to meet the proportional hazards assumption and the assumption cannot hold simultaneously for both types of proportional hazards models [19, 8]. Figure 2 shows the results in which age is made the time scale for the analysis and therefore requires left truncation. The model is adjusted for CD4 nadir, race, and time on study. The _{cs}HR(*t*) and _{sd}HR(*t*) are not proportional as evident by the non-constant HR over the entire range of ages. While a standard cause-specific proportional hazards could still be used, the interpretation would be complicated by the inclusion of time-interactions. Alternatively the average _{cs}HR could be calculated using an appropriate weighting scheme [26, 27, 28]. However, note that while this may be appropriate for initiation of therapy it appears less appropriate for progression to AIDS or death as the HR(*t*) is qualitatively different before (HR> 1) and after (HR< 1) age 47. This qualitative change in _{cs}HR for AIDS/Death such that history of injection drug use is protective at older ages unlikely to reflect a true protective effect (figure 2b), but is more likely to be a reduction in the number of susceptible individuals among those with history of injection drug use at least relative to those without a history of injection drug use. This illustrates the potential for selection bias in the HR [14]. Indeed the _{sd}HR remains elevated and increasing, indicating that the cumulative probability for disease progression to AIDS/death prior to therapy initiation remains higher among those with a history of injection drug use (Figure 2d).

The _{cs}HR(*t*) (panels a and b, left axis) and the _{sd}HR(*t*) (panels c and d, left axis) in which the time scale was age requiring left truncation in the analysis. The HRs compare those with a history of injection drug use for the initiation of therapy (left **...**

The cause-specific proportional hazards model for competing risk analysis has been used since the 1980s [3, 4]. Only with the relatively recent introduction of the subdistribution proportional hazards model by Fine and Gray [6] have regression methods evolved to estimate a hazard ratio that is directly interpretable with the cumulative incidence function (which may be estimable with left truncated data in the R package etm [41]). Several recent papers have discussed the effect of a covariate on either type of hazard may be quite different [7, 8, 19]. Our work demonstrates how a unifying parametric mixture model can be used to estimate both the cause-specific or subdistribution hazard functions and the corresponding _{cs}HR and _{sd}HR. Our example illustrates that the proposed parametric method may result in similar inferences and 95% confidence intervals. There are several advantages of the parametric mixture model over the semi-parametric approaches. First the parametric approach is not limited by the strong proportional hazards assumption. While this can be overcome by inclusion of time-covariate interactions this can complicate the interpretation. While incorporation of appropriate weights at each failure time may alleviate some of the interpretation difficulties in obtaining an average HR [26, 27, 28, 29], this approach may obscure qualitative differences in the HR(*t*). Second, the parametric mixture model allows for left truncation for estimating the _{sd}HR(*t*), while current semi-parametric methods for estimating the _{sd}HR(*t*) cannot handle left truncation. Finally, an advantage of all parametric methods over semi- and non-parametric approaches is that if an appropriate distribution has been chosen then parametric methods will be statistically more efficient.

Two methods for summarizing the hazard curves were presented. The “average” individual approach is convenient in that it is less computationally intensive. However, it can be counterintuitive to estimate the hazard function the “average covariate values” when covariates are binary (e.g. sex), substituting the mean value may not be consistent with the meaning of the data. Therefore, while the averages for each covariate are used, there is no guarantee that these are attained by any individual. The second empirically weighted approach uses the original data as it occurred with the alternative level of exposure. It overcomes the disadvantages of the average approach but requires more time in the construction of the estimate, with the replication of each individual to have a set of failure time points with and without exposure. However should the proportional hazard assumption be violated, it is preferable to present the hazard functions and HR curves as funtions of time rather than the summarized estimate. This will then convey information about how the association between exposure and outcome changes over time similar to methods we have previously proposed comparing two competing events [40]. The results from the semi-parametric and mixture model approaches applied to the data from the Women’s Interagency HIV Study yielded similar results. This is as one would expect given that a proportional hazards model essentially weights the ratio of the time-dependent hazards at each failure time-point equivalently and the mixture model was also similarly weighted. However, besides the advantage of providing the HR(*t*), the mixture model approach can easily be adapted to include a more appropriate weighting scheme that would provide a different summary result than the standard proportional hazards model [26, 27, 28, 29].

The mixture model approach is not the only parametric method for competing risks in the literature. Jeong and Fine [42, 43] provided a method that relies on the Gompertz distribution to directly parametrize the cumulative incidence function. The Gompertz distribution was chosen as it is an improper distribution and therefore a natural distribution for modeling the cumulative incidence function and may or may not provide a better fit than the mixture model [42]. However, as the authors state, the two-parameter Gompertz distribution only allows monotone hazard shapes. In addition, the cause-specific hazards are not estimable.

Furthermore, like many survival analysis methods, the mixture model presented assumes that eventually one of the competing events will occur (i.e. ). This may be a reasonable assumption for some situations (e.g. where all cause mortality is a competing event), but this assumption may be inappropriate in other cases in which indicating that some individuals are immune to the risks under investigation. In this situation, Maller and Zhou [13] have extended the mixture model to situations in which there may be individuals “immune” to the outcomes and may be of interest. Furthermore, a concern with all parametric methods is whether the chosen distribution is appropriate to describe the data. The use of a mixture of generalized gamma distributions allows the model to be flexible and includes some common distributions (such as exponential, log-normal, and Weibull distributions) as special cases of the generalized gamma family [20]. Additionally, because of the parametric nature of the model one can extrapolate beyond the point in which there is no additional data. However, one must be cautious in doing so to avoid over interpreting or over extending the inferences. Neverthless, semi-parametric and non-parametric methods are not able to do this. Moreover, the mixture model may be less robust than the semi-parametric proportional hazards model counterparts. This may be the case in which there are only a small number of the two competing events in which the addition of a few more of either of the events may change the inferences from the model. Of course this should not be a problem if there are a substantial number of events.

In conclusion, we have extended the use of a mixture model for competing risks. The mixture model provides a unifying approach to estimate both the cause-specific and subdistribution hazard functions as well as a summary estimate of the _{cs}HR and _{sd}HR. Estimates of the latter two are comparable to those obtained from proportional hazards models based upon either the cause-specific or subdistribution hazards. However, the mixture model has the advantage of being able to be flexible and provide an estimate the hazard ratios as a function of time or as a summary measure regardless of whether or not left truncation or interval censoring is present in the data.

We receive funding from the National Institutes of Health (U01-AI-42590 for the Women’s Interagency HIV Study; U01-AI069918 for the North American AIDS Cohort Collaboration on Research and Design, which is a part of the International Epidemiologic Databases to Evaluate AIDS (IeDEA); and K01-AI071754 (to B.L.)). The funding sources have had no involvement with this manuscript.

Contract/grant sponsor: National Institutes of Health (U01-AI-42590, U01-AI069918, K01-AI071754)

As previously stated, the likelihood function in equation (5) can be modified to allow for left truncation in order to allow estimation of the HR(*t*) while avoiding bias due to late entry into the study or for incorporating time-varying covariates into the survival models. We have discarded the covariate vector, **X**_{i1} and **X**_{i2}, from the equations to simplify the notation.

where *v _{i}* > 0 and

where *t*_{1i} < *t*_{2i} and *t*_{1i} and *t*_{2i} indicate the two boundary times of the interval in which the event *J _{i}* = 1 or

We declare that we have no conflict of interest.

1. Putter H, Fiocco M, Geskus RB. Tutorial in biostatistics: Competing risks and multi-state models. Statistics in Medicine. 2007;26:2389–2430. DOI: 10.1002/sim.2712. [PubMed]

2. Lau B, Cole SR, Gange SJ. Competing risk regression methods for epidemiologic data. American Journal of Epidemiology. 2009;170:244–256. DOI:10.1093/aje/kwp107. [PMC free article] [PubMed]

3. Prentice RL, Kalbleisch JD, Peterson AV, Jr., et al. The analysis of failure times in the presence of competing risks. Biometrics. 1978;34:541–554. [PubMed]

4. Kalbfleisch JD, Prentice RL. The Statistical Analysis of Failure Time Data. John Wiley and Sons, Inc.; New York, NY: 1980. pp. 247–277.

5. Gray RJ. A class of K-sample tests for comparing the cummulative incidence of a competing risk. The Annals of Statistics. 1988;16:1141–1154.

6. Fine JP, Gray RJ. A Proportional Hazards Model for the Subdistribution of a Competing Risk. Journal of the American Statistical Association. 1999;94:496–509.

7. Beyersmann J, Dettenkofer M, Bertz H, et al. A competing risks analysis of bloodstream infection after stem-cell transplantation using subdistribution hazards and cause-specific hazards. Statistics in Medicine. 2007;26:5360–5369. DOI: 10.1002/sim.3006. [PubMed]

8. Beyersmann J, Schumacher M. Letter to the editor: Misspecificed regression model for the subdistribution hazard of a competing risks. Statistics in Medicine. 2007;26:1649–1652. DOI: 10.1002/sim.2727. [PubMed]

9. Beyersmannn J, Latouche A, Buchholz A, Schumacher M. Simulating competing risks data in survival analysis. Statistics in Medicine. 2009;28:956–971. DOI: 10.1002/sim.3516. [PubMed]

10. Grambauer N, Schumacher M, Beyersmann J. Proportional subdistribution hazards modeling offers a summary analysis, even if misspecified. Statistics in Medicine. 2010;29:875–884. DOI: 10.1002/sim.3786. [PubMed]

11. Cox DR. The Analysis of Exponentially Distributed Lifetimes with 2 Types of Failure. Journal of the Royal Statistical Society Series B-Statistical Methodology. 1959;21:411–21.

12. Larson MG, Dinse GE. A mixture model for the regression analysis of competing risks data. Applied Statistics. 1985;34:201–211.

13. Maller RA, Zhou X. Analysis of parametric models for competing risks. Statistica Sinica. 2002;12:725–750.

14. Hernán M. The Hazards of Hazard Ratios. Epidemiology. 2010;21:13–15. DOI:10.1097/EDE.0b013e3181c1ea43. [PubMed]

15. Klein JP, Moeschberger ML. Survival Analysis: Techniques for Censored and Truncated Data. 2nd edition Springer Verlag; New York, NY: 2003.

16. Geskus RB. Cause-Specific Cumulative Incidence Estimation and the Fine and Gray Model Inder both Left Truncation and Right Censoring. Biometrics. 2010 in press. DOI: 10.1111/j.1541-0420.2010.01420.x. [PubMed]

17. Klein JP, Andersen PK. Regression Modeling of Competing Risks Data based on Pseudovalues of the Cumulative Incidence Function. Biometrics. 2005;61:223–229. DOI: 10.1111/j.0006-341X.2005.031209.x. [PubMed]

18. Andersen PK, Abildstrom SZ, Rosthøj S. Competing risks as a multi-state model. Statistical Methods in Medical Research. 2002;11:203–215. DOI: 10.1191/0962280202sm281ra. [PubMed]

19. Latouche A, Boisson V, Chevret S, Porcher R. Misspecified regression model for the subdistribution hazard of a competing risk. Statistics in Medicine. 2007;26:965–974. DOI: 10.1002/sim.2600. [PubMed]

20. Cox C, Chu H, Schneider MF, Muñoz A. Parametric survival analysis and taxonomy of hazard functions for the generalized gamma distribution. Statistics in Medicine. 2007;26:4352–4374. DOI: 10.1002/sim.2836. [PubMed]

21. Cox C. The generalized F distribution: an umbrella for parametric survival analysis. Statistics in Medicine. 2008;27:4301–4312. DOI: 10.1002/sim.3292. [PubMed]

22. Bates DM, Watts DG. Nonlinear Regression Analysis and Its Applications. Wiley; New York, NY: 1998.

23. Efron B. The Jacknife, the Bootstrap, and Other Resampling Plans. Society for Indiustrial and Applied Mathmatics; Philadelphia, PA: 1982.

24. Efron B, Tibshirani R. An introduction to the Bootstrap. Chapman Hall; London: 1993.

25. Lagakos SW, Schoenfeld DA. Properties of Propostional-Hazards Score Tests under Misspecified Regression Models. Biometrics. 1984;40:1037–1048. [PubMed]

26. Schemper M. Cox analysis of survival data with non-proportional hazard functions. The Statistician. 1992;41:455–465.

27. Kalbfleisch JD, Prentice RL. Estimation of the average hazard ratio. Biometrika. 1981;68:105–112.

28. Sasieni P. Maximum Weighted Partial Likelihood Estimators for the Cox Model. Journal of the American Statistical Association. 1993;88:144–152.

29. Schemper M, Wakounig S, Heinze The estimation of average hazard ratios by weighted Cox regression. Statistics in Medicine. 2009;28:2473–2489. DOI: 10.1002/sim.3623. [PubMed]

30. Barkan SE, Melnick SL, Preston-Martin S, et al. The Women’s Interagency HIV Study. WIHS Collaborative Study Group. Epidemiology. 1998;9:117–25. [PubMed]

31. Celentano DD, Vlahov D, Cohn S, et al. Self-reported antiretroviral therapy in injection drug users. JAMA. 1998;280:544–546. [PubMed]

32. Celentano DD, Galai N, Sethi AK, et al. Time to initiating highly active antiretroviral therapy among HIV-infected injection drug users. AIDS. 2001;15:1707–1715. [PubMed]

33. Lucas GM, Cheever LW, Chaisson RE, et al. Detrimental effects of continued illicit drug use on the treatment of HIV-1 infection. J Acquir Immune Defic Syndr. 2001;27:251–259. [PubMed]

34. Rodriguez-Arenas MA, Jarrin I, del Amo J, et al. Delay in the initiation of HAART, poorer virological response, and higher mortality among HIV-infected injecting drug users in Spain. AIDS Res Hum Retroviruses. 2006;22:715–723. DOI: 10.1089/aid.2006.22.715. [PubMed]

35. Keiser O, Taffe P, Zwahlen M, et al. All cause mortality in the Swiss HIV Cohort Study from 1990 to 2001 in comparison with the Swiss population. AIDS. 2004;18:1835–1843. [PubMed]

36. Lau B, Gange SJ, Moore RD. Risk of non-AIDS-related mortality may exceed risk of AIDS-related mortality among individuals enrolling into care with CD4 counts greater than 200 cells/mm3. J Acquir Immune Defic Syndr. 2007;44:179–187. DOI: 10.1097/01.qai.0000247229.68246.c5. [PubMed]

37. Leemis LM, McQueston JT. Teacher’s Corner: Univariate Distribution Relationships. The American Statistician. 2008;62:45–53. DOI: 10.1198/000313008X270448.

38. Cole SR, Hernán MA. Adjusted survival curves with inverse probability weights. Computer Methods and Programs in Biomedicine. 2004;75:45–49. DOI: doi:10.1016/j.cmpb.2003.10.004. [PubMed]

39. Xie J, Liu C. Adjusted Kaplan-Meier estimator and log-rank test with inverse probability of treatment weighting for survival data. Statistics in Medicine. 2005;24:3089–110. DOI: 10.1002/sim.2174. [PubMed]

40. Lau B, Cole SR, Moore RD, Gange SJ. Evaluating competing adverse and beneficial outcomes using a mixture model. Statistics in Medicine. 2008;27:4313–4327. SOI: 10.1002/sim.3293. [PMC free article] [PubMed]

41. Allignol A, Schumacher M, Beyersmann J. A Note on Variance Estimation of the Aalen-Johansen Estimator of the Cumulative Incidence Function in Competing Risks, with a View towards Left-Truncated Data. Biometrical Journal. 2010;52:126–137. DOI: 10.1002/bimj.200900039. [PubMed]

42. Jeong JH, Fine J. Direct parametric inference for the cumulative incidence function. Journal of the Royal Statistical Society: Series C (Applied Statistics) 2006;55:187–200. DOI: 10.1111/j.1467-9876.2006.00532.x.

43. Jeong JH, Fine J. Parametric regression on cumulative incidence function. Biostatistics. 2007;8:184–196. DOI: 10.1093/biostatistics/kxj040. [PubMed]

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's Canada Institute for Scientific and Technical Information in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |