Application of additive hazards models for analyzing survival of breast cancer patients
Parisa Ataee Dizaji1, Mahtab Vasheghani Farahani1, Ayeh Sheikhaliyan2, Akbar Biglarian3
1 Department of Biostatistics, University of Social Welfare and Rehabilitation Sciences, Tehran, Iran
2 Department of Industrial Engineering, Malek Ashtar University of Technology, Tehran, Iran
3 Department of Biostatistics, Social Determinants of Health Research Center, University of Social Welfare and Rehabilitation Sciences, Tehran, Iran
|Date of Submission||03-Nov-2019|
|Date of Decision||28-Jan-2020|
|Date of Acceptance||29-Jun-2020|
|Date of Web Publication||28-Oct-2020|
Dr. Akbar Biglarian
Department of Biostatistics, Social Determinants of Health Research Center, University of Social Welfare and Rehabilitation Sciences, Tehran
Source of Support: None, Conflict of Interest: None
Background: Survival rates for breast cancer (BC) are often based on the outcomes of this disease. The aim of this study was to compare the performance of three survival models, namely Cox regression, Aalen's, and Lin and Ying's additive hazards (AH) models for identifying the prognostic factors regarding the survival time of BC patients. Materials and Methods: This study was a historical cohort study which used 1025 females' medical records that underwent modified radical mastectomy or breast saving. These patients were admitted to Besat and Chamran Hospitals, Tehran, Iran, during 2010–2015 and followed until 2017. The Aalen's and Lin and Ying's AH models and also traditional Cox model were applied for analysis of time to death of BC patients using R 3.5.1 software. Results: In Aalen's and also Lin and Ying's AH models, age at diagnosis, history of disease, number of lymph nodes, metastasis, hormonal therapy, and evacuation lymph nodes were prognostic factors for the survival of BC patients (P < 0.05). In addition, in the Lin and Ying's AH model tumor size (P = 0.048) was also identified as a significant factor. According to Aalen's plot, metastasis, age at diagnosis, and number of lymph nodes had a time-varying effect on survival time. These variables had a different slope as the times go on. Conclusion: AH model may yield new insights in prognostic studies of survival time of patients with BC over time. Because of the positive slope of estimated cumulative regression function in Aalen's plot, metastasis, higher age at diagnosis, and high number of lymph nodes are important factors in reducing the survival BC, and then based on these factors, the therapists should consider a special therapeutic protocol for BC patients.
Keywords: Additive hazards model, breast cancer, metastasis, survival analysis
|How to cite this article:|
Ataee Dizaji P, Vasheghani Farahani M, Sheikhaliyan A, Biglarian A. Application of additive hazards models for analyzing survival of breast cancer patients. J Res Med Sci 2020;25:99
|How to cite this URL:|
Ataee Dizaji P, Vasheghani Farahani M, Sheikhaliyan A, Biglarian A. Application of additive hazards models for analyzing survival of breast cancer patients. J Res Med Sci [serial online] 2020 [cited 2020 Nov 23];25:99. Available from: https://www.jmsjournal.net/text.asp?2020/25/1/99/299489
| Introduction|| |
Nowadays, cancer is one of the most important health problems in the world. It is a serious and prevalent problem that all of the countries of the world are faced with., In 2018, 18.1 million cancer cases and 9.6 million cases of death due to cancer were estimated from around the world. In 2015, breast, lung, and colorectal cancers were the most common cancers,, and in 2018, breast, colorectal, and lung cancers were the most commonly diagnosed cancers.
Breast cancer (BC) is the most common type of cancer diagnosed in women and is the main cause of cancer-related death in females in 2015, and also in 2018. BC was ranked second and fifth in terms of incidence and mortality worldwide, respectively. It alone comprises 11.6% of all cancers and 6.6% of cancer-related deaths, the leading cause of cancer-related death in over 100 countries worldwide.
BC has a high incidence rate in the world. This cancer is growing in South America, Africa, and Asia. BC is increasing in incidence among Southeast Asian women and also West Europe.
In Iran, cancer is the third leading cause of deaths after injuries and cardiovascular diseases, accounting for more than 53,000 annual deaths. Current evidences suggest that stomach, breast, prostate, leukemia, and lung are the most common incident cancers in both sexes in Iran. BC is the most common cancer occurring among Iranian women as well as in those from other countries. Similar to that in the Middle East countries, the mean age of BC diseases in Iran is 10 years lower than that in the developed countries.
The survival of patients with BC was dependent on access to medical facilities and BC screening. In the recent decade, prevention, early diagnosis, and proper treatment had improved survival rate of cancer.,,
The main goal of survival analysis is to investigate and model the probable relationship between risk factors and the death time of a patient. When the response variable usually is time to death, disease recurrence, or metastasis, there are three major approaches to analysis of the data: (a) Cox proportional hazards as a semiparametric; (b) parametric functions such as Weibull, Gompertz, and log-normal; and (c) nonparametric methods such as Aalen's additive risk model and Kaplan–Meier method. Survival analysis has traditionally been performed using the Cox model and preferred by most researchers. This model is simple to fit, and the results are easy to explain, but there are two potential drawbacks: proportionality and ignoring the effect of time-varying covariate.,, In contrast, the Aalen's additive model is nonparametric and propounds time-varying covariates., Unlike the PH model which estimates hazard ratios, in the additive model, the difference in hazards is used to describe the association between the risk factors and disease occurrence.,
This study aimed to compare the performance of three survival models, namely Cox regression, Aalen's, and Lin and Ying's additive hazards (AH) models for identifying the prognostic factors regarding the survival time of BC patients.
| Materials and Methods|| |
This study was a historical cohort study which used data gathered from medical records of 1025 BC patients who were admitted to Besat and Chamran Hospitals, Tehran, Iran, during 2010–2015 and followed until 2017. Females with a confirmed diagnosis of BC who underwent either modified radical mastectomy (MRM) or breast saving (BS) from 2010 to 2015 were enrolled in the study. This study approved by the Ethics Committee, affiliated to University of Social Welfare and Rehabilitation Sciences (IR.USWR.REC.1399.050).
Data were gathered from patients' medical records and histopathology reports by trained data collectors. Age at diagnosis, number of lymph nodes, tumor size, family history (no, yes), Stage (I, II, III, IV), Grade (I, II, III), metastasis (no, yes), human receptor of epidermal growth factor 2 (minus, plus), excapsular (no, yes), evacuation lymph nodes, history of disease (no, yes), estrogen receptors (minus, plus), progesterone receptors (minus, plus), pathology report (ductal carcinoma in situ, invasive carcinoma, both), hormonal therapy (no, yes), and surgery (MRM, BS) variables were gathered through patient's medical records. Finally, time to death was considered from the date of surgery (MRM/BS) to the date of current status (death/censoring).
At first, preprocessing of the data was made on Excel software, and the cases with missingness in these variables or unknown current status were excluded. After that, almost all of the patients had chemotherapy and radiotherapy. Hence, these variables exclude from the analysis. After that, the sample size was decreased to 933. Subsequently, by performing a univariate analysis, all variables with P < 0.2 were selected for modeling. By this strategy, age at diagnosis, number of lymph nodes, tumor size, metastasis, excapsular, evacuation lymph nodes, history of disease, hormonal therapy, and kind of surgery variables were selected for multiple analysis.
Finally, the Aalen's and Lin and Ying's AH models and also Cox proportional hazards model, as a traditional model, were applied for data analysis. In this way, the Schoenfeld residual plot was used to evaluate the proportionality assumption of Cox model. A nonzero slope is an indication of a violation of the proportional hazard assumption. Moreover, 95% pointwise confidence intervals were used to estimate cumulative regression functions based on Aalen's additive model. An estimated cumulative regression function plot (Aalen's plot) is obtained to see the effect of covariates over time. The smoothed line in the Aalen's plot with a slope was approximately equal to zero, suggesting that there may be no time-varying effect, and this is in agreement with the test. The slope of an estimated cumulative regression function is positive when covariate increases, and this fact corresponds to an increasing hazard rate. On the other hand, if the slope is negative while the covariate increases, then this fact points to a decreasing hazard rate. Data analysis was carried out with survival and addhazard packages in R 3.5.1 software. In this sense, coxph, survfit, ah, and aareg functions were used to fit the models. P < 0.05 was considered statistically significant.
Cox proportional hazards model
In the analysis of censored failure time observations, the standard Cox proportional hazards model assumes that the regression coefficients are time independent. This model, also known as the Cox model, has the form:
where λ (t | Z) is the hazard at time t given predictor values Z = (Z1,Z1,…., Zp) and λ0(t) is an arbitrary baseline hazard function and β is unknown time-independent parameters.
In the Cox proportional hazards model, the effect of the covariates acts multiplicatively on the baseline hazard rate, and the hazard for each covariate is assumed to be constant over time.
| Lin and Ying's Additive Hazards Model|| |
This model is semiparametric AH, and the effect of covariates is additive on the baseline hazard rate. According to this model, hazard function for failure time (Ti) has the form:
Where λ0(t) is baseline hazard function and β is unknown time-independent additive effects.,
Aalen's additive hazards model
This model is nonparametric AH, and the effect of covariates is additive on the baseline hazard rate. According to this model, hazard function for failure time (Ti) has the form:
Where β(t) is unknown time-dependent additive effects and Z (t) is a (P × 1) vector of covariates. This model is useful when the main interest is risk difference rather than relative risk, and the model allows covariate effects to vary with time. In addition, this model provides a cumulative regression function plot to display the change of covariates over time.,
| Results|| |
The mean age of patients at diagnosis was 50.71 (standard deviation [SD] = 11.30) years, and the mean of tumor size was 3.03 (SD = 1.85) cm. Of the 933 patients with BC, 90.2% were without metastasis, 78.8% had hormonal therapy, and 54.1% had mastectomy surgery. Other characteristics were reported in [Table 1].
|Table 1: Demographic and clinical characteristics of breast cancer patients|
Click here to view
According to Cox proportional hazards model, the age at diagnosis (P< 0.001), history of disease (P< 0.001), number of lymph nodes (P< 0.001), metastasis (P< 0.001), hormonal therapy (P< 0.001), and evacuation lymph nodes (P = 0.004) were identified as significant factors [Table 2]. For example, tumor size coefficient was positive, which indicates a worse prognosis effect, and this variable is associated with increased risk of BC (hazard = 1.10). Based on the scaled Schoenfeld residual plots, it seems that the proportionality is not exactly satisfied for metastasis, age at diagnosis, and hormonal therapy [Figure 1]. Hence, Aalen's additive model was used to analyze the data. The age at diagnosis (P< 0.001), history of disease (P< 0.001), number of lymph nodes (P = 0.030), metastasis (P< 0.001), hormonal therapy (P = 0.001), and evacuation lymph nodes (P = 0.011) had a significant effect on survival time. The other covariates had no effect on BC's lifetime. For example, the tumor size coefficient was 4.54e-4, which indicates the absolute difference in the BC rate per-unit change in the tumor size. In Lin and Ying's AH model, the age at diagnosis (P< 0.001), history of disease (P< 0.001), number of lymph nodes (P = 0.026), metastasis (P< 0.001), hormonal therapy (P = 0.002), tumor size (P = 0.048), and evacuation lymph nodes (P = 0.012) were identified as significant factors [Table 2]. For example, the coefficient of tumor size was 4.54e-4, which indicates that patients with per-unit change in the tumor size had an increase in hazard of 0.0005.
|Table 2: Results of fitting Aalen's additive hazards, Cox PH, and Lin and Ying's additive hazards models for breast cancer data|
Click here to view
|Figure 1: Smoothed scaled Schoenfeld residual plots with 95% pointwise confidence intervals for covariates|
Click here to view
[Figure 2] indicates that the estimates of cumulative regression function for all variables used in this data set. For patients who had history of disease, or gave hormonal therapy, or evacuation Lymph node; the risk decreases as the time goes on, i.e., the slope of the graph is negative for this patient during 80 month. For patients with metastasis, patients with higher age, and patients with high number of lymph nodes, the slope of an estimated cumulative regression function was positive. This shows that the risk of BC increases over time for these patients.
|Figure 2: Estimated cumulative regression functions with 95% pointwise confidence intervals based on Aalen's additive model; which x label is time in month and y label is cumulative coefficients (i.e., risk at time)|
Click here to view
| Discussion|| |
As we know, statistical models can describe the nature of data from any scientific field such as medicine, clinical studies, and health care. In this sense, one can make statistical inferences about the parameters in order to predict goals and also decision-making for some interventions in human health., Survival analysis is used to model the relationship between the risk factors in time-to-death cancer studies.,, In this field, the use of the Cox PH model is more popular. In this model, hazard ratios are estimated which are useful in understanding the magnitude of association between exposure and survival. Comparing to this model, AH models provide the difference in hazards, i.e., the change in the event rate due to the exposure variable (Xie, 2013).
Many studies were made on survival of Iranian BC patients.,,,, The overall relative survival in Iranian BC patients was reported higher than some Asian and eastern European countries and was lower than the U.S. and western European countries.
In this study, the Cox PH, the Lin and Ying's AH, and the Aalen's AH models were applied to BC data. The interpretation of additive models and Cox proportional hazards model is very different, so that the exponential of coefficients in Cox proportional hazards model is relative hazards, but those in additive models are the attributable risk.
Results of fitting the BC data with these models were similar with little difference in Lin and Ying's AH model. The Standard Error (SE) in Aalen's additive model was better than the Cox model and Lin and Ying's AH model, so Aalen's AH model is better than to two other mentioned models. Based on Aalen's AH model, age at diagnosis, number of lymph nodes, and evacuation lymph nodes had a significant effect on the hazard of the event. In addition, history of disease, hormonal therapy, and metastasis had a significant effect on the hazard of the event.
In the present study, age at diagnosis of BC was an important factor affecting the survival of patients. Some researches demonstrate that a relationship exists between age and the proportion of BC. Wei et al. showed that in the group of very young patients (age < 35 years), there were more triple-negative tumors than older patients. Largillier et al. reported a poor prognosis for BC patients over 50 years of age. In China, the prognostic mathematical model of lymph node-negative BC was firstly established in 2003. Compared with them, our study showed that the number of lymph nodes was a risk factor of BC. Our results showed that the patients who had not history of disease, the risk of disease decreases for them as the time goes on. This is close to results of some studies in Iran which they showed that history of disease was a significant factor of BC.,,, Metastasis was the significant risk factor in Aalen's AH model which we used in this study. This factor was also reported as a significant factor in other studies. The overall survival of BC after metastasis was shorter reported and also was reported as a factor which increases the hazard of the event of death.
Some studies suggested a linear and others suggested a nonlinear effect of tumor size. It is mentioned that some studies reported that tumor size increased the risk of metastasis in patients with BC,, but in our study, the relationship between tumor size and survival of BC patients only in Aalen's AH model has been controversial.
One of the limitations of this study was missingness in some clinical information such as grade and nodal involvement. In addition, there was not the survival time of some patients because of changing in contact information. So some samples be excluded from the analysis because of this missingness. Moreover, it was assumed that the patient censoring was not informative and thus was independent of the BC death. Another limitation of this study was the lack of time-varying scheme data, and hence, the extended Cox model was not used. On the other hand, our data could not support the use of the extended Cox model. It is suggested that future studies be conducted longitudinally to allow for data collection over time which makes a useful data source.
| Conclusion|| |
The current study showed that the age at diagnosis, number of lymph nodes, and evacuation lymph nodes had increased the effect on the hazard of the event. In addition, history of disease, hormonal therapy, and metastasis also increased the effects on the hazard of the event. For patients who had history of disease, or gave hormonal therapy, or evacuation Lymph node; the risk decreases as the time goes on, i.e., the slope of the graph is negative for this patient during 80 month. For patients with metastasis, patients with higher age, and patients with high number of lymph nodes, the slope of an estimated cumulative regression function was positive. This shows that the risk of BC increases over time for these patients. This study applied Aalen's AH and Lin and Ying's models beside Cox proportional hazards model to analyze BC patients. The results of these models were similar. The Aalen's AH model has a similar fit than other models but gives information about covariate effect when time goes on. This aspect of the association between the survival time and covariates is a useful and helpful interpreter for therapists and clinical researchers.
The authors are thankful to all managers and personnel of the Besaat and Chamran Hospitals in Tehran who cooperated in performing this study. We also wish many thanks to Mrs. Sheikh Alian for making few months of data collection. In addition, the authors wish to express their gratitude to the Deputy of Research and Technology of University of Social Welfare and Rehabilitation Sciences for financial support, with research ID number 2311.
Financial support and sponsorship
This study was supported by the Deputy of Research and Technology of University of Social Welfare and Rehabilitation Sciences in Tehran, Iran.
Conflicts of interest
There are no conflicts of interest.
| References|| |
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2015. CA Cancer J Clin 2015;65:5-29.
Torre LA, Bray F, Siegel RL, Ferlay J, Lortet-Tieulent J, Jemal A. Global cancer statistics, 2012. CA Cancer J Clin 2015;65:87-108.
Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin 2018;68:394-424.
Ghoncheh M, Pournamdar Z, Salehiniya H. Incidence and Mortality and Epidemiology of Breast Cancer in the World. Asian Pac J Cancer Prev 2016;17:43-6.
Ferlay J, Soerjomataram I, Dikshit R, Eser S, Mathers C, Rebelo M, et al
. Cancer incidence and mortality worldwide: Sources, methods and major patterns in GLOBOCAN 2012. Int J Cancer 2015;136:E359-86.
Nguyen L, Brunicardi FC, Dibardino DJ, Scott BG, Awad SS, Bush RL, et al
. Education of the modern surgical resident: Novel approaches to learning in the era of the 80-hour workweek. World J Surg 2006;30:1120-7.
Fallahzadeh H, Momayyezi M, Akhundzardeini R, Zarezardeini S. Five year survival of women with breast cancer in Yazd. Asian Pac J Cancer Prev 2014;15:6597-601.
Abedi G, Janbabai G, Moosazadeh M, Farshidi F, Amiri M, Khosravi A. Survival rate of breast cancer in Iran: A meta-analysis Asian Pac J Cancer Prev 2016;17:4615-21.
Ganz PA, Goodwin PJ. Breast cancer survivorship: Where are we today? Adv Exp Med Biol 2015;862:1-8.
Martinussen T, Vansteelandt S. On collapsibility and confounding bias in Cox and Aalen regression models. Lifetime Data Anal 2013;19:279-96.
Klein JP. Modelling competing risks in cancer studies. Stat Med 2006;25:1015-34.
Madadizadeh F, Ghanbarnejad A, Ghavami V, Zare Bandamiri M, Mohammadianpanah M. Applying additive hazards models for analyzing survival in patients with colorectal cancer in fars province, Southern Iran Asian Pac J Cancer Prev 2017;18:1077-83.
Atashgar K, Sheikhaliyan A, Biglarian A. Reviewing and analyzing breast cancer survival models. Iran J Surg 2016;24:26-40.
Therneau TM, Grambsch PM. Modeling Survival Data: Extending the Cox Model. New York: Springer; 2000.
Kleinbaum DG, Klein M. Survival Analysis: A Self-Learning Text. 3th
ed. New York: Springer-Verlag; 2012.
Breslow NE. and Day NE. Statistical methods in cancer research. In: The Analysis of Case-Control Studies. IARC Scientific Publication No. 32. Vol. 1. Geneva, Switzerland: World Health Organization; 1980.
Aalen OO. A linear regression model for the analysis of life times. Stat Med 1989;8:907-25.
Abadi A, Saadat S, Yavari P, Bajdik C, Jalili P. Comparison of Aalen's additive and Cox proportional hazards models for breast cancer survival: Analysis of population- based data from British Columbia, Canada. Asian Pac J Cancer Prev 2011;12:3113-6.
Lin DY, Ying Z. Semiparametric analysis of general additive-multiplicative hazard models for counting processes. Ann Statist 1995;23:1712-34.
Zhao S, Hu T, Ma L, Wang P, Sun J. Regression analysis of informative current status data with the additive hazards model. Lifetime Data Anal 2015;21:241-58.
Kutner M, Nachtsheim C, Neter J, Li W. Applied Linear Statistical Models.5th
ed. New York: McGraw-Hill/Irwin; 2004.
Friedman LM, Furberg CD, DeMets D, Reboussin DM, Granger CB. Fundamentals of Clinical Trials. 5th
ed.. Switzerland: Springer International Publishing; 2015.
Xie X, Strickler HD, Xue X. Additive hazard regression models: An application to the natural history of human papillomavirus. Comput Math Methods Med 2013;2013:796270.
Bakhshi E, Sheikhaliyan A, Atashgar K, Kooshesh M, Biglarian A. Survival analysing of the breast cancer patients using cure model. Iran Red Crescent Med J 2017;19:e55575.
Sharifian A, Pourhoseingholi MA, Emadedin M, Rostami Nejad M, Ashtari S, Hajizadeh N, et al
. Burden of breast cancer in iranian women is increasing. Asian Pac J Cancer Prev 2015;16:5049-52.
Wei XQ, Li X, Xin XJ, Tong ZS, Zhang S. Clinical features and survival analysis of very young (age& lt; 35) breast cancer patients. Asian Pac J Cancer Prev 2013;14:5949-52.
Largillier R, Ferrero JM, Doyen J, Barriere J, Namer M, Mari V, et al
. Prognostic factors in 1,038 women with metastatic breast cancer. Ann Oncol 2008;19:2012-9.
Scully OJ, Bay BH, Yip G, Yu Y. Breast cancer metastasis. Cancer Genomics Proteomics 2012;9:311-20.
Gohari MR, Khodabakhshi R, Shahidi J, Fard ZM, Foadzi H, Soleimani F, et al
. The impact of multiple recurrences in disease-free survival of breast cancer: An extended Cox model. Tumori 2012;98:428-33.
Bijker N, Peterse JL, Duchateau L, Julien JP, Fentiman IS, Duval C, et al
. Risk factors for recurrence and metastasis after breast-conserving therapy for ductal carcinoma-in-situ: Analysis of European Organization for Research and Treatment of Cancer Trial 10853. J Clin Oncol 2001;19:2263-71.
Dawood S, Broglio K, Esteva FJ, Ibrahim NK, Kau SW, Islam R, et al
. Defining prognosis for women with breast cancer and CNS metastases by HER2 status. Ann Oncol 2008;19:1242-8.
[Figure 1], [Figure 2]
[Table 1], [Table 2]