Skip to main content

Survival analysis and influence of the surgical aggression of a cohort of orthopedic and trauma patients in a non-controlled spread COVID-19 scenario



Determining the infection rate and mortality probability in healthy patients who have undergone orthopedic and trauma surgeries (OTS) during a period of uncontrolled COVID-19 transmission may help to inform preparations for future waves. This study performed a survival analysis in a cohort of non-infected OTS patients and determined the effect of COVID-19 on mortality.


This observational study included 184 patients who underwent OTS in the month before surgical activities ceased and before the implementation of special measures. Four groups of surgery (GS) were established based on the location of the surgery and the grade of inflammation produced. Crude risk of infection and infection rates were assessed. Survival and failure functions by GS were analyzed. Comparison of the Kaplan-Meier survival curves by GS was assessed. Cox regression and Fine-Gray models were used to determine the effect of different confounders on mortality.


The crude risk of COVID-19 diagnosis was 14.13% (95% CI: 9.83–19.90%). The total incidence rate was 2.67 (1000 person-days, 95% CI: 1.74–3.91). At the end of follow-up, there was a 94.42% chance of surviving 76 days or more after OTS. The differences in K-M survivor curves by GS indicated that GS 4 presented a lower survival function (Mantel-Cox test, p = 0.024; Wilcoxon-Breslow test, p = 0.044; Tarone-Ware test, p = 0.032). One of the best models to determine the association with mortality was the age-adjusted model for GS, high blood pressure, and respiratory history, with a hazard ratio of 1.112 in Cox regression analysis (95% CI: 1.005–1.230) and a sub hazard ratio of 1.111 (95% CI: 1.046–1.177) in Fine-Gray regression analysis for competitive risk.


The infection risk after OTS was similar to that of the general population in a community transmission area; the grade of surgical aggression did not influence this rate. The survival probability was extremely high if patients had not previously been infected. With higher grades of surgical aggression, the risk of mortality was higher in OTS patients. Adjusting for age and other confounders (e.g., GS, high blood pressure and respiratory history) was associated with higher mortality rates.

Peer Review reports


On March 11, 2020, the World Health Organization (WHO) declared the novel coronavirus disease 2019 (COVID-19) as a global pandemic. By May 2020, COVID-19 had spread globally [1]. During the first and second (ongoing at the time of writing) waves, Spain became one of the countries with the highest number of cases and Madrid, one of the most affected regions [2].

During the first and second infection waves, elective and urgent orthopedic and trauma surgery (OTS) ceased, as many hospitals were operating with grade 3 COVID-19 occupancy scores [3]. However, normal orthopedic healthcare should be restored between infection waves and in the future. It is therefore essential to determine the risk of virus transmission and mortality following elective and urgent orthopedic and trauma surgery (OTS) in non-infected patients in a community with high COVID-19 transmission levels. Furthermore, it is important to determine whether surgical aggression grade effects these risks.

There are few reports in the literature that examine the risks of infection after performing elective and urgent OTS in a community with high levels of COVID-19 transmission. While several studies have analyzed mortality in this context, reporting high rates [3,4,5], the majority of patients in these studies had already been infected or were in the incubation period when the surgery was performed.

This study aimed to analyze the risk of infection and perform a survival analysis in a cohort of healthy patients who underwent OTS during a period of uncontrolled COVID-19 transmission – at the start of the first wave, before any control measures had been implemented. Infection and survival rates in OTS patients at this time can be compared to those during periods in which measures have been implemented. The results of this study could be used to guide decision-making between infection waves when elective or trauma surgeries need to be introduced again and to determine which patient and type of surgery has the lowest and highest risk.


Study design

This is an observational study with an ambispective cohort design [6]. It adhered to the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) checklist [7]. Patients who underwent OTS between February 13 and March 13, 2020 were included in this study. Once the study population was defined in a retrospective manner, the cohort was followed prospectively until May 5, 2020.


This study was carried out at the Hospital Universitario Infanta Leonor, located in the south-east side of Madrid City, Spain and covering a target population of 305262 patients. The hospital was at grade 3 COVID-19 occupancy at the time of the study [3]. On March 13, 2020, the occupancy rate reached 202%, and the hospital was declared a “100% COVID-19 hospital” – no other medical or surgical activity apart from the treatment of COVID-19 patients could be carried out.

Eligibility criteria

We included all patients who had undergone OTS in the month before surgical activities ceased. During this period, the state of “alert” had not yet been declared; thus, COVID-19-specific infection prevention and control measures had not been established and the rate of transmission was high. Surgical activity was therefore completely normal during this time, mimicking periods between waves or when a new one is about to begin. Patients were recruited by MC and DG via phone call. Prior to the phone interview, patients were asked for their consent to participate in the study.

Data collection

General patient data were obtained from the hospital database. Disease histories were obtained from the preoperative evaluation performed by the anesthesiology department. Only patients who were admitted to the hospital had their outcome data recorded on the electronic medical record. For the patients who were not admitted, outcomes were confirmed through a structured phone call interview, performed within 2 weeks for all patients to avoid any discrepancy in exposure (Additional file 1).


The study strictly followed the ethical principles of biomedical research. The use of clinical/personal patient data in the research was approved by the Ethics Committee of the hospital (R-026-20).


The inclusion criteria were pre-anesthetic assessment approval (no suspicion of COVID-19 or COVID-19 exposure during the preoperative anesthetic evaluation) and an age over 18 years. The exclusion criteria were the presence of symptoms compatible with COVID-19 and the absence of consent to participate in the study.

Variables and data sources

The primary outcomes of the study were diagnosis of COVID-19 and death from COVID-19. Case definitions were classified as recommended by the Spanish Ministry of Health [8] and the WHO [9]. The independent and potential confounders were: “Age”, “Sex”, “Body Mass Index” (BMI), “High Blood Pressure” (HBP), “Diabetes” (DM), “History of Respiratory Disease” (RespHist), “History of Cardiological Disease” (CardiolHist), “History of Kidney Disease” (RenalHist), “American Society of Anesthesiologists Physical Status Classification System” (ASA), and “Group of Surgery” (GS).

The “Group of Surgery” variable was based on the grade of aggression and inflammation. Four different groups were established based on the location of the surgery and the grade of the surgery taking the NICE guidelines into account [10]. Groups 1 and 2 included procedures performed in small joints (foot, ankle, hand and wrist) while Groups 3 and 4 included those performed on big joints (elbow, shoulder, knee and hip). The difference between Groups 1–2 and 3–4 was the grade of aggression of the procedure (1 and 3: soft tissue and arthroscopic procedures; 2 and 4 bone and arthroplasty procedures; Table 1).

Table 1 Groups of Surgery based on the location of the surgery and aggression and inflammation grade

Statistical analysis

Crude rates

The crude risk and risk ratio (RR) of COVID-19 diagnosis were determined by case definition and GS, and the corresponding 95% confidence intervals (CI) were calculated. The crude incidence rate (person-day) and crude incidence rate ratio (IRR) of COVID-19 diagnosis by GS and their 95% CI were obtained. RR and IRR were analyzed using GS 1 as a reference.

Survival and hazard function analysis

The survival function of the entire cohort was analyzed by Kaplan-Meier (K-M) [11] and life-table actuarial methods. Survival and failure functions by GS were analyzed by K-M and competing risk models [12]. Comparison of the K-M survival curves by GS was assessed using the Mantel-Cox log rank test (M-C) [13], the Wilcoxon-Breslow test [14], and the Tarone-Ware test (T-W) [15], which followed the chi-squared distribution under the null hypothesis (Ho) that the survival functions are equal and a level of significance (α) of 0.05.

Quality of follow-up

The distribution of the frequencies of the different states (Alive, Dead, Lost) at the end of the follow-up was calculated. The follow-up period for the different states at the last observation in the cohort was analyzed as another indicator of the quality of the follow-up. Finally, the minimum, maximum, mean, median, 25th, and 75th percentiles were assessed by the different states at the end of the follow-up.

Multivariate survival analysis for measuring the effect

Selection of independent variables

All independent variables (IVs) were assessed by univariate Cox regression analysis and those with p ≤ 0.20 were considered to be confounders [16]. BMI and ASA were excluded based on the univariate analysis (complete data in Additional file 2). “RespHist” (p = 0.989) with a non-significant hazard ratio (HR) was included in the multivariate analysis because of its theoretical importance. The final IVs selected were as follows: “Age,” “GS,” “HBP,” “DM,” “RespHist,” “CardiolHist,” and “RenalHist.”

Maximum model (MMax)

Based on previous publications [17], the variable “Age” was selected as the main predictor variable and the rest of the IVs were considered to be possible confounders. The final reference model (MMax) was built up based on the confounders and their interaction with the predictor variable.

Assessment of the interactions

The interactions were assessed using a chunk test, which compared the reference model (MMax) to the model without the interactions (MMaxNoInteract). It was based on the likelihood ratio statistic that followed the chi-squared distribution with the degrees of freedom equal to k (number of IVs) minus the number of interactions, and a level of significance of 0.05 (α = 0.05; complete data in Additional file 2).

Assessment of the confounders

In line with the Maldonado and Greenland recommendations [16], confounders were assessed based on the change in their effect or HR being less than 10% compared with the reference model. The final models were compared using the Stata user command “confound” [18] developed for modeling confounding in linear, logistic, and Cox regression (complete data in Additional file 2).

Final model assessment

The final selected models adjusted by confounders were analyzed using a Cox regression model and a Fine-Gray Regression model [19,20,21] for competing risks, which showed the effect based on the sub hazard ratio (SHR).

Diagnosis of the model

Two assumptions are required in the Cox proportional hazard model: the proportional assumption and the log-linear relationship. The proportionality assumption assumes that the effect of the predictors on the hazard rate is constant throughout the follow-up time. To verify this assumption, the interaction of each predictor with the survival time was added to the Cox model chosen to measure the effect. The null hypothesis was established, which assumed that the coefficients of interaction terms were statistically equal to zero. Finally, the relationship between the Schoenfeld residuals and survival time was analyzed using a chi-squared test, assuming proportionality when p > 0.05. The log-linear relationship assumption of the Cox model assumes that the relationship between the instantaneous IR and the explanatory variables is log-linear. The analysis of the squared linear predictor was used in this way to confirm that the squared predictor coefficient was not significant (p > 0.05).


Only one patient was excluded from this study because verbal consent was not obtained. In total, 184 patients were recruited for this study, none of whom dropped out at follow-up. Demographic data and the distribution of the IVs are shown in Table 2.

Table 2 Demographic data and independent variables distribution

The total number of missing values was 0.41% (9/184 missing BMI values); thus, it is unlikely that missing values affected the results [22].

The quality of the follow-up was classified as excellent. No patients had dropped out of the study by the end of the last observation. The mean follow-up time was 32.8 days.

The overall crude risk of COVID-19 diagnosis according to the four different GS is shown in Table 3. The relative risk of COVID-19 diagnosis by GS, taking GS 1 as reference, were not significant independently, both including or excluding the suspected cases (complete data in Additional file 3).

Table 3 Risks of COVID-19 diagnosis

The total IR and the crude IR (1000 person-days) of COVID-19 diagnosis by GS are shown in Table 4. There were no statistical differences in the IRR when using GS 1 as reference.

Table 4 Incidence rate and incidence rate ratio of COVID-19 diagnosis

Of the 26 patients with COVID-19 diagnosis, 16/184 (8.7%) presented with mild symptoms, stayed at home, and did not require any specific drug treatment while 5/184 (2.72%) were admitted to hospital due to more severe symptoms; all patients survived. There were 7/184 (3.80%) deaths in the entire cohort, 5 of them (2.72%) were deaths that fit the criteria of COVID-19 deaths and 2 (1.09%) were deaths not caused by COVID-19 (complete data in Additional file 4).

The cumulative survival probabilities of the entire cohort calculated using the K-M method and actuarial method are shown in Table 5 (complete data in Additional file 5). At the end of the follow-up, there was a 94.42% chance of surviving 76 days or more after OTS. The actuarial method showed that most of the COVID-19 deaths occurred during the first 1.5 months after surgery. The cumulative probability of survival between 60 and 75 days was 0.957.

Table 5 Survival function of entire cohort

The failure function (cumulative incidence of mortality) in the entire cohort calculated by using the K-M method and competing risk method are shown in Table 6 (complete data in Additional file 6). Our data suggest that there was approximately 5.5% chance of dying by COVID-19 in the first 71 days after OTS.

Table 6 Failure function. Entire cohort

The K-M failure function by GS showed that the risk of mortality at 76 days or earlier was 0.1563 (Table 7; complete data in Additional file 6). The comparison of K-M and life tables of survival (actuarial method) curves by GS are shown in Figs. 1 and 2. The differences in survival curves by GS were statistically significant (M-C test, p = 0.024; W-B test, p = 0.044; T-W test, p = 0.032) and determined that GS 4 had a lower survival function compared to the remaining GS (complete data in Additional file 7). One of the five COVID-19 patients who died had undergone arthroplasty for a shoulder fracture, while the remaining four had had hip fractures.

Table 7 Failure function by Group of Surgery
Fig. 1
figure 1

Kaplan-Meier survival by Group of Surgery

Fig. 2
figure 2

Proportion of survivors by Group of Surgery

The selected models, adjusted by confounders and analyzed using a Cox regression model (effect based on HR) and Fine-Gray Regression (effect based on SR) model, demonstrated that “Age” adjusted for “HBP” showed the highest association with mortality by COVID-19 (HR: 1.145; SHR: 1.144). The model “Age” adjusted for “GS,” “HBP,” and “RespHist” was also significantly associated with mortality by COVID-19, albeit to a lesser extent (HR and SHR: 1.111; Table 8). The diagnosis of the models demonstrated that the proportional hazard assumption was accepted based on Schoenfeld residuals (p = 0.405) and the test of interaction between survival time and the independent variables (p = 0.857). The log linear assumption was also accepted (p = 0.406; complete data in Additional file 8).

Table 8 Final selected models adjusted by confounders


This study demonstrated that the risk and IR of COVID-19 infection in patients who underwent OTS were similar for the different GSs. However, the cumulative survival probability in patients who had bone and arthroplasty procedures in big joints (GS 4) was significantly lower compared to the rest of the surgery groups; there was a 15.63% chance of dying by COVID-19 by day 76 following surgery. Age adjusted by GS, HBP, and RespHist was associated to the mortality. At the end of the follow-up, there was a 94.42% chance of surviving 76 days or more after OTS.

The risk of COVID-19 diagnosis (considering all types of diagnose) was 14.13% (95% CI: 9.83 to 19.90%) for the entire cohort. In Madrid, which has a population of 6663394 [23], there were 62989 confirmed cases (0.95%) on May 5, 2020; 40851 had been admitted to hospital (64.86% of all cases; 0.61% of the entire population) and 8420 had died (13.37% of all cases; 0.13% of the entire population) [24]. A seroprevalence study performed during the same period in Spain determined that the rate of infection was 11.3% (95% CI: 9.8 to 13.0) and 19.7%, with the consideration of suspected cases. Considering these data, the risk of diagnosis, including suspected cases, in our cohort is similar to the general population.

The relative risk of all COVID-19 diagnoses by the different GSs was not statistically significant when compared using group 1 as the reference. Thus, the grade of aggression and inflammation cannot be considered as a factor of a higher risk of infection in our cohort.

Notably, the most interesting variable to analyze is the probability of survival. In our study, the chance of survival at 76 days or more in a patient who had OTS was 94.42% (95% CI: 81.72–98.38). This high survival probability provides a more positive outlook when compared to results obtained in previous studies. Indeed, Lei et al. [4] reported a mortality rate of 20.5% after elective surgery; the COVID Surg Collaborative group showed a crude 30-day mortality rate of 23.8% in the entire cohort and of 28.8% in the orthopedic surgery cohort. In those studies, all patients had developed the disease by day 5 following surgery, or the COVID-19 infection was confirmed preoperatively thus, they had already been infected. This may be the most important difference; the patients in our study were healthy (not infected) before surgery. The mortality probability after having an OTS, even in an uncontrolled transmission area, is low if the patient has not been infected previously.

Regarding the influence of the aggression and inflammation of the OTS on mortality, there were no deaths in GS 1, 2 and 3. All deaths were concentrated in GS 4 and the lower survival function was statistically significant in this group. Lei et al. [4] concluded that a more complex surgery increases the risk of having the worst evolution; however, they detected several variables (older age and comorbidities) that are more frequent in this group, although they did not include a multivariate model to adjust for confounders. The higher mortality rate in older patients has been reported in several studies with a high number of patients [25,26,27,28], mirroring our model in which age was the most important variable. However, most of these studies considered age in addition to the risk factors from unadjusted studies for other confounders [29]. Age adjusted for the GS, HBP, and RespHist was one of the best models for explaining the association of mortality by COVID-19 in patients who underwent OTS, with a confirmed proportional hazard and log-linear assumptions that required the model diagnosis.

This study has some limitations. We could not define a nonsurgical subgroup in our cohort during the same recruitment period and not all types of OTS, such as spine surgery or pediatric surgery, were included as they are performed in another hospital. Representativeness of the general OTS population and changes in case definitions should be considered when making inferences from the risk and incidence rates of infection by COVID-19 observed in our cohort [30]. This is important because the rates of infection and mortality can be compared in different scenarios and with different infection prevention measures in place. Furthermore, this study is limited by the fact that COVID-19 case definition criteria and testing are still evolving.

The highlight of this observational study is the fact that the patients (who did not have COVID-19) underwent OTS during a period of uncontrolled COVID-19 spread; thus, these data could help clarify the risk of infection and survival probability and aid decision-making between infection waves.


The infection risk after OTS is similar to that of the general population in a community transmission area; the grade of surgical aggression did not influence this rate. The survival probability is extremely high if patients have not previously been infected with COVID-19. Mortality is higher if the group had a higher aggression and inflammation grade. Age adjusted for confounders such as GS, high blood pressure, and respiratory history is associated with mortality.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request. Most of the materials can be found in the additional material.



Orthopedic and trauma surgeries


Strengthening the Reporting of Observational Studies in Epidemiology


Body Mass Index


High Blood Pressure


Diabetes Mellitus


History of Respiratory Disease


History of Cardiological Disease


History of Kidney Disease


American Society of Anesthesiologists Physical Status Classification System


Group of Surgery


National Institute for Health and Care Excellence


Confidence intervals


Incidence rate


Relative risks


Incidence rate ratio




Mantel-Cox log rank test


Tarone-Ware test


Independent variables


Hazard ratio


Maximum model


Model without the interactions


Sub hazard ratio


  1. World Health Organization. Coronavirus disease (COVID-19) Situation Report– 106. Available from: Accessed 25 May 2020.

  2. Ministry of Health. Spain. COVID-19 evolution in Spain. Available from:ón-geográfica. Accessed 25 Sept 2020.

  3. Muñoz Vives JM, Jornet-Gibert M, Cámara-Cabrera J, Esteban P, Brunet L, Delgado-Flores L, et al. Mortality rates of patients with proximal femoral fracture in a worldwide pandemic. J Bone Joint Surg. 2020;102(13):e69.

    Article  PubMed  Google Scholar 

  4. Lei S, Jiang F, Su W, Chen C, Chen J, Mei W, et al. Clinical characteristics and outcomes of patients undergoing surgeries during the incubation period of COVID-19 infection. EClinicalMedicine. 2020;21:100331.

    Article  PubMed  PubMed Central  Google Scholar 

  5. COVID Surg Collaborative. Mortality and pulmonary complications in patients undergoing surgery with perioperative SARS-CoV-2 infection: an international cohort study. Lancet. 2020;396(10243):27–38.

    Article  Google Scholar 

  6. Kleinbaum DG, Kupper LL, Morgenstern H. Epidemiologic research: principles and quantitative methods. New York: Wiley; 1982.

    Google Scholar 

  7. Vandenbroucke JP, von Elm E, Altman DG, Gøtzsche PC, Mulrow CD, Pocock SJ, et al. Strengthening the reporting of observational studies in epidemiology (STROBE): explanation and elaboration. Int J Surg. 2014;12(12):1500–24.

    Article  PubMed  Google Scholar 

  8. Ministerio de Sanidad. Estrategia de diagnóstico, vigilancia y control en la fase de transición de la pandemia de covid-19 indicadores de seguimiento. Accessed 6 May 2020.

  9. World Health Organization. Global Surveillance for human infection with coronavirus disease (COVID-19). Accessed 2 May 2020.

  10. NICE guideline. Routine preoperative tests for elective surgery. Accessed 11 May 2020.

  11. Kaplan EL, Meier P. Nonparametric estimation from incomplete observations. J Am Stat Assoc. 1958;53(282):457–81.

    Article  Google Scholar 

  12. Gooley TA, Leisenring W, Crowley J, Storer BE. Estimation of failure probabilities in the presence of competing risks: new representations of old estimators. Stat Med. 1999;18(6):695–706.<695::AID-SIM60>3.0.CO;2-O.

    Article  CAS  PubMed  Google Scholar 

  13. Mantel N. Evaluation of survival data and two new rank order statistics arising in its consideration. Cancer Chemother Rep. 1966;50(3):163–70.

    CAS  PubMed  Google Scholar 

  14. Breslow N. A generalized Kruskal-Wallis test for comparing K samples subject to unequal patterns of censorship. Biometrika. 1970;57(3):579–94.

    Article  Google Scholar 

  15. Tarone RE, Ware J. On distribution-free tests for equality of survival distributions. Biometrika. 1977;64(1):156–60.

    Article  Google Scholar 

  16. Maldonado G, Greenland S. Simulation study of confounder-selection strategies. Am J Epidemiol. 1993;138(11):923–36.

    Article  CAS  PubMed  Google Scholar 

  17. Wynants L, Van Calster B, Bonten MMJ, Collins GS, Debray TPA, De Vos M, et al. Prediction models for diagnosis and prognosis of covid-19 infection: systematic review and critical appraisal. BMJ. 2020;369:m1328.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Doménech JM, Navarro JB. Find the best subset for linear, logistic and cox regression: user-written command confound for Stata [computer program]. V1.1.7. Bellaterra: Universitat Autònoma de Barcelona; 2020.

    Google Scholar 

  19. Li C. The Fine-Gray model under interval censored competing risks data. J Multivar Anal. 2016;143:327–44.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Fine JP, Gray RJ. A proportional hazards model for the subdistribution of a competing risk. J Am Stat Assoc. 1999;94(446):496–509.

    Article  Google Scholar 

  21. Fine JP. Regression modeling of competing crude failure probabilities. Biostatistics. 2001;2(1):85–97.

    Article  CAS  PubMed  Google Scholar 

  22. Rosales RS, Atroshi I. The methodological requirements for clinical examination and patient-reported outcomes, and how to test them. J Hand Surg Eur Vol. 2020;45(1):12–8.

    Article  PubMed  Google Scholar 

  23. Instituto Nacional de Estadística. Población por provincias y tamaño de los municipios españa.!tabs-tabla. Accessed 22 May 2020.

  24. Actualización no 96. Enfermedad por el coronavirus (COVID-19). Accessed 22 May 2020.

  25. Yang J, Zheng Y, Gou X, Pu K, Chen Z, Guo Q, et al. Prevalence of comorbidities and its effects in coronavirus disease 2019 patients: a systematic review and meta-analysis. Int J Infect Dis. 2020;94:91–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Zhou F, Yu T, Du R, Fan G, Liu Y, Liu Z, et al. Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study. Lancet. 2020;395(10229):1054–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Chen T, Wu D, Chen H, Yan W, Yang D, Chen G, et al. Clinical characteristics of 113 deceased patients with coronavirus disease 2019: retrospective study. BMJ. 2020;368:m1091.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Wu C, Chen X, Cai Y, Xia J, Zhou X, Xu S, et al. Risk factors associated with acute respiratory distress syndrome and death in patients with coronavirus disease 2019 pneumonia in Wuhan, China. JAMA Intern Med. 2020;180(7):934–43.

    Article  CAS  PubMed  Google Scholar 

  29. Jordan RE, Adab P, Cheng KK. Covid-19: risk factors for severe disease and death. BMJ. 2020;368:m1198. PMID: 32217618.

    Article  PubMed  Google Scholar 

  30. Tsang TK, Wu P, Lin Y, Lau EHY, Leung GM, Cowling BJ. Effect of changing case definitions for COVID-19 on the epidemic curve and transmission parameters in mainland China: a modelling study. Lancet. 2020;5(5):e289–96.

    Google Scholar 

Download references


Not applicable.


No funding was received.

Author information

Authors and Affiliations



DG and MC acquired the data and helped with the analysis. Participate in revising it critically for important intellectual content. FC and RSR designed the study, performed the analysis, interpreted the data, and drafted the manuscript. R L-G Made substantial contributions to the conception and design of the study, significantly revised the manuscript. All authors read and approved the final manuscript and Agreed to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.

Corresponding author

Correspondence to Fernando Corella.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the Research Ethics Commission of the Infanta Leonor University Hospital and Virgen de la Torre Hospital. (Reference number 026–20).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Interview guide. Document use by the researches to guide the telephone interview.

Additional file 2.

Multivariate analysis. The complete STATA data are shown for the following: A Selection of variables. Univariate Cox regression analysis. B Analysis of the confounders and interactions. C Confound assessment.

Additional file 3.

Relative risks of COVID-19 diagnosis by Group of Surgery. The complete STATA data are shown for the following: A COVID diagnosis (suspected + probable + confirmed cases). B COVID diagnosis (only probable + confirmed cases).

Additional file 4.

Distribution of cases in the entire cohort. Table showing the distribution of the evolution of cases in the entire cohort (non-infected; alive; COVID-19 deaths; non-COVID-19 deaths; severe disease).

Additional file 5.

Cumulative survival probabilities of the entire cohort according to the K-M method and the actuarial method. The complete STATA data are shown for the following: A Survival list. Entire cohort. B Actuarial method.

Additional file 6.

Failure function (cumulative incidence of mortality) in the entire cohort according to the K-M method and the competing risk. The complete STATA data are shown for the following: A Survival function by Group of Surgery. B K-M failure function. Entire cohort. C Failure function. Competing risk. D K-M failure function by Group of Surgery. E Comparison of life tables of survival (actuarial method) by Group of Surgery.

Additional file 7.

Differences in survival curves by Group of Surgery. The complete STATA data are shown for the following: Log rank, Breslow, and Tare tests for comparing survivor curves by Group of Surgery.

Additional file 8.

Model assessment. The complete STATA data are shown for the following: A Schonefeld residuals. B Test based on the interaction between survival-time and the independent variables. C Log linear assumption. D Quality of the follow-up.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Corella, F., Rosales, R.S., Guzman Domenech, D. et al. Survival analysis and influence of the surgical aggression of a cohort of orthopedic and trauma patients in a non-controlled spread COVID-19 scenario. BMC Musculoskelet Disord 22, 594 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: