Abstract
Objective To estimate the prevalence of disability and anxiety in Covid-19 survivors at discharge from hospital and analyze relative risk by exposures.
Design Multi-center retrospective cohort study.
Setting Twenty-eight hospitals located in eight provinces of China.
Methods A total of 432 survivors with laboratory-confirmed SARS CoV-2 infection participated in this study. At discharge, we assessed instrumental activities of daily living (IADL) with Lawton’s IADL scale, dependence in activities of daily living (ADL) with the Barthel Index, and anxiety with Zung’s self-reported anxiety scale. Exposures included comorbidity, smoking, setting (Hubei vs. others), disease severity, symptoms, and length of hospital stay. Other risk factors considered were age, gender, and ethnicity (Han vs. Tibetan).
Results Prevalence of at least one IADL problem was 36.81% (95% CI: 32.39-41.46). ADL dependence was present in 16.44% (95% CI: 13.23-20.23) and 28.70% (95% CI: 24.63-33.15) were screened positive for clinical anxiety. Adjusted risk ratio (RR) of IADL limitations (RR 2.48, 95% CI: 1.80-3.40), ADL dependence (RR 2.07, 95% CI 1.15-3.76), and probable clinical anxiety (RR 2.53, 95% CI 1.69-3.79) were consistently elevated in survivors with severe Covid-19. Age was an additional independent risk factor for IADL limitations and ADL dependence; and setting (Hubei) for IADL limitations and anxiety. Tibetan ethnicity was a protective factor for anxiety but a risk factor for IADL limitations.
Conclusion A significant proportion of Covid-19 survivors had disability and anxiety at discharge from hospital. Health systems need to be prepared for an additional burden resulting from rehabilitation needs of Covid-19 survivors.
Introduction
With rapidly increasing numbers of confirmed cases and deaths caused by the Coronavirus Disease-2019 (Covid-19) pandemic, clinical outcomes research has almost exclusively focused on disease progression and mortality[1]. With the case fatality rate estimated to lie somewhere between 0.56 and 9.38 percent globally (95% prediction interval, I2 = 100%),[2] it is obvious however that a vast majority of people actually survive Covid-19. A substantial proportion of survivors have experienced severe disease episodes requiring hospitalization and - in a considerable number of cases - intensive care including mechanical ventilation[3]. In addition, contracting a potentially fatal disease is stressful and can cause anxiety and other psychiatric manifestations in many patients.[4] It would be naïve to assume that health problems just disappear upon discharge from acute care or when viral RNA is no longer detectable, and equally naïve to take planned and systematic follow up of patients for granted.[5, 6] Yet, the burden of disease resulting from physical and psychological sequelae of Covid-19 in patients surviving the disease as well as potential long-term adverse effects from treatments including steroids such as Dexamethasone are just beginning to catch attention[7, 8]. Globally, respiratory infections and tuberculosis assume already the fourth rank of all diseases in terms of disability adjusted life years (DALYs) and the 17th rank in years lived with disability (YLD). In infectious disease, they sadly rank first and second, respectively.[9, 10] Covid-19 will increase this burden considerably and cause longer term mental and physical health problems, work disability, and reduced quality of life in survivors, frontline health professionals, and those quarantined.[11, 12] What we know from the 2002/2003 Severe Acute Respiratory Syndrome (SARS) epidemic is not encouraging. In Hong Kong for example, about 80 percent of SARS survivors still needed follow-up treatment two years after the outbreak had been contained.[13] The number of patients needing follow-up care and rehabilitation due to Covid-19 will be unprecedented. Countries need to consider and plan for this additional long-term challenge to the health system. It is important to understand risk factors for disability and mental health problems now in order to direct scarce resources and meet windows of opportunity for treatment. This study makes a first attempt to address this issue. We aimed to (1) estimate the prevalence of disability and anxiety in Covid-19 survivors from eight Provinces/centrally governed municipalities of the PR China at discharge from acute inpatient treatment, and (2) investigate relative risk of adverse outcomes by various determinants including gender, age, comorbidity, setting, ethnicity and disease severity. Prevalence estimates are important to gain insight in the potential number of patients needing follow-up care. Analyzing risk factors is instrumental in identifying vulnerable groups and allocating resources for early interventions.
Methods
Study design
We conducted a multi-center retrospective cohort study. The study protocol was approved by the ethics committee of West China Hospital, Sichuan University (2020-273).
Subjects
Four hundred and thirty-two Covid-19 survivors were surveyed who had received treatment in 28 designated hospitals in Hubei, Sichuan, Guizhou, Henan, Neimenggu (inner Mongolia), Shandong, Hainan, and Chongqing from January 18 to March 15, 2020 and gave informed consent. Included were adults (≥ 16 years) who had laboratory-confirmed SARS coronavirus-2 (SARS CoV-2) infection and met the diagnostic criteria of the Novel Coronavirus Pneumonia Diagnosis and Treatment Program (Trial Version 7) issued by the National Health Commission of the PR China (Third edition).[14] Survivors further needed to be in stable medical condition and have been discharged or be about to be discharged from the participating centers. Excluded were patients with pneumonia caused by other types of coronaviruses, seasonal influenza, bacteria, or other not SARS CoV-2 etiology.
Patients were contacted by health professionals who were in charge of their treatment. All health professionals were instructed on implementation of the study protocol, informed consent procedures, and possible sources of bias by videoconference.
Measures
Outcomes
Outcomes were disability and anxiety. Cut off date for outcome assessment was April 7, 2020. Disability was measured in two ways: (1) Limitations in instrumental activities of daily living (IADL) were assessed with the Lawton IADL scale[15] and classified into no limitations vs. one or more limitations (primary outcome). An index of the number of reported IADL limitations was also created. (2) Dependence in activities of daily living (ADL) was measured with the Barthel Index[16] and classified into severe or moderate dependence vs. mild dependence/independence.[17] Anxiety was evaluated with Zung’s Self-Reported Anxiety Scale[18] (SAS) and classified into probable clinical anxiety disorder vs. not according to a conservative cutoff recommended for research.[19] Outcome data were patient-reported and collected with an online questionnaire.
Risk factors
Clinical and demographic data were extracted from hospital records by participating health professionals. Current smoking status was self-reported by participants. Demographic data included age, gender, province (Hubei vs. others), and ethnicity (Han vs. Tibetan). Having received treatment in Hubei was considered a risk factor because of the particular situation in Hubei which was the epicenter of the epidemic in China representing about 80 percent of confirmed cases. Ethnicity was included as a significant proportion of the sample from Sichuan was from the Tibetan ethnic group, and belonging to a minority ethnic group is known to be associated with health outcomes in other disease. Clinical data included comorbidity, disease severity, symptoms at admission, scope of pneumonia, and length of inpatient stay. Recorded comorbid conditions included COPD, chronic bronchitis, pulmonary fibrosis, hypertension, diabetes, cardiovascular disease, kidney disease, hepatitis, rheumatism, and gout. Degree of Covid-19 severity was classified according to the Chinese standard.[14] Patients were defined as severe cases when they met one of the following criteria at any time during hospitalization: acute respiratory distress, respiratory rate ≥30 breath/min; pulse oxygen saturation (SpO2) ≤93% at rest; arterial blood partial pressure of oxygen/fraction of inspired oxygen (PaO2/FiO2) ≤300 mmHg (1 mmHg = 0.133 kPa); respiratory failure requiring mechanical ventilation; septic shock; failure of other organs requiring ICU treatment. Symptoms recorded at admission included fever, cough, diarrhea, fatigue and pain. Scope of pneumonia was noted as uni- or bilateral per radiographic findings. Data on length of hospital stay were available for 408 patients (94.44%).
Sample size calculation
Sample size to detect a two-fold relative risk for reporting one or more IADL limitations in the severe disease group as compared to the non-severe group with a power of 80% and alpha error of 5% was estimated under the following assumptions: 0.2 ratio of severe to non severe group (based on patients received at our own center), 15 percent prevalence of outcome in non-severe group. This yielded a minimal sample size of 401 (67 severe, 334 non-severe).
Analysis
Prevalence of limitations in IADL, ADL dependence, and anxiety for the full sample and by disease severity were provided with 95% confidence intervals (CIs) estimated from logistic regression. Risk ratios (RR) and 95% CIs for gender, age group (<= 50, 50-60, > 60 years), province (others vs. Hubei), ethnicity (Han vs. Tibetan), comorbid health conditions (none, one, multiple), disease severity (non-severe vs. severe), symptoms (<=1, >1) at admission, and infection scope (unilateral vs. bilateral) were estimated with log-linear Poisson regression with robust standard errors.[20] Since longer time of being hospitalized could increase anxiety, length of hospital stay was included in the anxiety-model in addition. For each outcome unadjusted RR and RR adjusted for all other potential risk factors evaluated were provided. Zero-inflated Poisson (ZIP) regression was used to estimate the effect of risk factors on the extent of IADL limitations. ZIP regression assumes that there are two groups in the population from which the data are sampled, one which always has zero counts, and one which may have zero or higher counts. Young people with non-severe disease who have no comorbidity may never develop IADL limitations, while other groups may develop (additional) IADL limitations due to factors associated with disease progression. ZIP regression simultaneously performs logistic and Poisson regression making two types of predictions: (1) prediction of excess zeros assuming a binary dependent variable (no IADL limitations vs. one or more), (2) prediction of the count portion (number of experienced IADL limitations). Predictors included in the logistic part were the following background factors: age group and interaction of age and disease severity, province, ethnicity, smoking, and comorbid conditions. Predictors included in the count part were potential drivers of disease consequences: disease severity, interaction of age and disease severity, symptoms, and infection scope. Vuong’s test confirmed fit of the ZIP model.
All analyses were performed with Stata 14 (Stata corporation, Texas, USA).
Results
Demographic and clinical characteristics and prevalence of outcomes of the study participants are provided in Table 1. Median age was 49 (IQR 35-60) and 49% were female. About one third of the study participants were classified as severe cases. Fever, cough, and fatigue were the most common symptoms at admission. A majority had no pre-existing health condition and were non-smokers. Patients with severe Covid-19 were older, were more often from Hubei and belonged to the Han majority, presented with more symptoms, had more often bilateral pneumonia, and stayed in the hospital for a longer time. More pre-existing health conditions were present in severe cases, but the percentage of smokers was lower. Over one third of Covid-19 survivors had at least one IADL limitation, about 15 percent had at least moderate ADL dependence, and six percent were severely dependent. Probable clinical anxiety disorder was found in about 29 percent of the population. Prevalence of all outcomes was more than four times higher in survivors of severe Covid-19 than in those with non-severe disease.
Results on risk factors for disability and anxiety are displayed in Fig 1 and 2 (for detailed estimates see S1-3 Tables). Risk ratios for reporting at least one IADL limitation (Fig 1, upper panel) were increased in the older age groups, particularly in survivors aged older than 60 years, in patients from Hubei, and in survivors with severe Covid-19. The adjusted model also showed an elevated relative risk in survivors from the Tibetan ethnic group. Only in unadjusted models was the relative risk of having one or more IADL limitations clearly elevated in survivors with comorbidities, in patients with more symptoms at admission, and in those with bilateral infections. Both in unadjusted and adjusted models, the relative risk for ADL dependence (Fig 1 lower panel) was consistently increased in Covid-19 survivors aged older than 60 years, and in participants with severe disease course and more symptoms. All individual risk ratios were higher in unadjusted models. Setting, comorbidity, and infection site played a role in unadjusted analysis only.
Risk ratios and 95% confidence intervals. Estimated from log-linear Poisson regression with robust standard errors; unadjusted risk ratios are estimated from univariable model for the effect of the respective predictor on the outcome, adjusted risk ratios are estimated from multivariable model containing all potential risk factors. IADL = Instrumental Activities of Daily Living. ADL = Activities of Daily Living.CI=confidence interval. Note: In the case of ADL dependence, the upper limit of the confidence interval for the unadjusted estimate of age > 60 years has been truncated for better readability of the other effects (see S1-3 Tables for details). Adjusted model and unadjusted model for age have been estimated for 431 cases (because of one missing date of birth).
Risk Ratios and 95% confidence intervals. Estimated from log-linear Poisson Regression with robust standard errors; unadjusted risk ratios are estimated from univariable model for the effect of the respective predictor on the outcome, adjusted risk ratios are estimated from multivariable model containing all potential risk factors. The adjusted model has been estimated for 407 cases (because of missing information regarding length of stay for 24 cases and one case with missing date of birth). &Unadjusted model estimated for 407 cases (because of missing information regarding length of stay for 24 cases and one case with missing date of birth). CI=confidence interval.
Having severe Covid-19 was the strongest risk factor for probable clinically relevant anxiety, followed by having received treatment in Hubei province (Fig 2). In turn, belonging to the Tibetan group was associated with a largely decreased relative risk of anxiety. There was a trend for an increased relative risk of anxiety in survivors who had stayed in the hospital for more than 14 days. Older age, comorbidity, symptoms, and scope of pneumonia played a role in unadjusted analysis only.
Table 2 show the results of the ZIP regression of the number of IADL limitations on risk factors. Hospitalization in Hubei, Tibetan ethnicity, and multi-comorbidity were associated with lower odds of having a zero count of IADL limitations. For all age groups, odds for having zero IADL limitations were decreased in people with severe Covid-19. This effect became however smaller with increasing age. Disease severity, its interaction with age, and symptoms were associated with a greater number of health problems. The joint effect of age and degree of severity of Covid-19 on the number of IADL limitations from both parts of the model is illustrated in Fig 3. In all age groups having a severe course of Covid-19 was associated with a higher number of self-reported IADL limitations. Point estimates for survivors with non-severe and those with severe disease were however closer in the oldest age group. CIs for the number of IADL limitations in participants with non-severe and those with severe Covid-19 moved closer in both older age groups and overlapped in survivors aged 61 years and older.
Estimated from zero-inflated Poisson regression. Predictions are adjusted for other factors in the model. Predictors included in the logistic part (excess zeros) were the following: age group, interaction of age and disease severity, setting (Hubei vs. other province), Tibetan ethnicity (reference: Han), smoking history, and comorbid conditions. Predictors included in the count part were: disease severity, interaction of age and disease severity, symptoms, and infection scope. The model has been estimated for 431 cases (because of one missing date of birth).
Discussion
This retrospective cohort study found a considerable prevalence of functional limitations and dependence as well as anxiety in Covid-19 survivors at about time of discharge from acute inpatient treatment. Disease severity was a major independent risk factor for all outcomes. Age was an additional risk factor for both disability outcomes, and setting (Hubei) for IADL limitations and anxiety. Tibetan ethnicity was a strong protective factor for anxiety but a risk factor for IADL limitations when it was adjusted for covariates. Demographic and clinical characteristics of the study population were comparable to other reports on Covid-19 patients from China [1, 21] with two exceptions: a higher percentage was classified as severe cases and current smoking was more prevalent in the non-severe disease group. An explanation for the higher percentage of severe cases in the present study is that the applied criteria for classifying Covid-19 cases as severe by the National Health Commission of the PR China are more liberal than the criteria for severe community-acquired pneumonia by the American Thoracic Society and Infectious Disease Society of America[22] used by Guan et al.[21] The lower prevalence of smokers in participants with severe disease was unexpected and contradicts previous evidence in this regard.[21] A possible explanation is that current smoking status was self-reported and participants with more severe disease outcomes did not reveal their true smoking status, perhaps for fear of being stigmatized and held responsible for contracting SARS CoV-2 or the severity of Covid-19.
We found a high prevalence of disability and anxiety, particularly in severe cases. Almost three quarters of severe Covid-19 survivors experienced at least one IADL problem and more than one third had at least moderate ADL dependence. The six percent of the overall sample with severe ADL dependence are of particular concern since severe functional dependence has been a reliable predictor of mortality in populations with other health conditions.[23] No other original study directly investigated disability outcomes related to Covid-19 so far. Follow up studies of SARS and Middle East Respiratory Syndrome (MERS) survivors have however shown that decreased exercise capacity and distance walked, fatigue, sleeping problems, shortness of breath, reduced self-rated health, work disability, and structural damage (e.g. osteonecrosis) from corticosteroids were prevalent in survivors even years after onset of symptoms.[11, 24, 25] In addition, access to health care, quality of services, and employment opportunities of former SARS patients were compromised by enduring stigma[13, 26]. Our findings imply that a large proportion of Covid-19 survivors, particularly those who experienced severe disease episodes are in need of follow up care and rehabilitation services. This proportion may be reduced by the integration of early rehabilitation interventions into acute care[27, 28].
Probable clinical anxiety disorder was found in one third of the overall sample and its prevalence was almost four-fold increased in severe as compared to non-severe cases. Previous studies on mental health in the Covid-19 pandemic have been conducted in the general population[29, 30] and health professionals[31] but not in the patient population. Studies in SARS survivors have however reported highly elevated long-term psychiatric morbidity[32] and demonstrated adverse effects on the mental health of caregivers.[11] Our findings point to the need of psychological counseling during acute care and mental health follow-up in many Covid-19 survivors to mitigate suffering and prevent the manifestation of psychiatric disorders.
Major risk factors for mortality due to Covid-19[1, 33] were also major risk factors of disability: older age and severe disease course both independently increased the relative risk of IADL limitations and ADL dependence. However, comorbidity[34] had an adverse effect on disability outcomes in unadjusted analysis only. Though attenuated, disease severity continued to play a major role when it was adjusted for age, co-morbidity, and other potential confounders demonstrating an additional effect of Covid-19 when other known risk factors for disability were held statistically constant. In all age groups having severe Covid-19 was associated with a higher number of self-reported IADL limitations. For survivors aged over 60 years we found that the extent of IADL limitations was more similar in those with non-severe and severe Covid-19 than in the other age groups. Accordingly, the oldest participants with non-severe Covid-19 were more similar to survivors with severe disease in the younger age groups. There are two possible explanations for this phenomenon that deserve further scrutiny: First, older people have already more activity limitations before they acquire Covid-19. Second, a non-severe disease course has a greater effect on IADL limitations for them.
Persistence of this effect when it was adjusted for comorbidity speaks for the latter explanation. Having been treated in hospitals of Hubei province where Covid-19 spread first and overwhelmed an unprepared health system was furthermore associated with adverse disability outcomes, particularly IADL limitations. Moreover, belonging to the Tibetan ethnic minority predicted worse IADL outcomes when it was controlled for other potential risk factors. Given that the Tibetan population of Sichuan province mainly lives in rural, mountainous areas, both findings point at the role that environmental factors including the availability and accessibility of targeted services play in the disablement process.[35] Our findings indicate that priority in rehabilitation resource allocation should be given to patients with severe disease course, older people, and those living in environments with restricted accessibility of services.
Regarding anxiety, outstanding independent risk factors were severity of Covid-19 and setting, i.e. having been hospitalized in Hubei province. Apart from the specific situation in hospitals in Hubei, insecurity about causes, transmission patterns, and prognosis when Covid-19 first occurred likely contributed to this finding. Moreover, we found a trend for longer length of hospital stay being independently associated with increased relative risk of anxiety. In contrast to disability outcomes, Tibetan ethnicity was a consistent protective factor in unadjusted as well as adjusted models of anxiety. Possible explanations are cultural factors including larger families (due to the non application of China’s one child policy to this group) and greater social cohesion and support. This finding needs further exploration. Our results point towards heightened need for psychological support in patients who experienced severe disease episodes and people from areas hardest hit by the Covid-19 pandemic.
Study Limitations
The present study has limitations. First, baseline data for the outcomes analyzed here were not available. Because of this we cannot rule out that patients had disability or anxiety before they contracted Covid-19, indeed having disability may have made people more vulnerable to the disease [36-38]. However, most patients first presented already with severe symptoms and baseline information collected before admission is rare in any disease. Baseline information on comorbidity is thus as close as we get to baseline disability and the effects of disease severity on outcomes remained consistent when adjusted for baseline comorbidity in our analysis. Second, we may be criticized for not having employed some kind of control group. An appropriate control group for the effects and population under investigation is however difficult to define and to establish. Age and gender matched general population controls are surely not appropriate because of the lack of hospitalization. Investigating patients who have been hospitalized for some other disease would not help disentangle the effects of disease and hospitalization. Quarantined persons suspected for Covid-19 or those treated at home may make up an interesting comparator but are difficult to access. We believe that our comparison of patients with severe and non-severe Covid-19 nonetheless provides initial evidence for the potential impact of the disease and is an important starting point in investigating this population’s rehabilitation needs. Third, outcome data were self-reported by patients due to the lack of options for objective measurement. Survivors thus may have under- or overestimated their physical abilities. Patient-reported outcome measures have however not been used in this population before to the best of our knowledge, making it difficult to appraise if and to what degree this kind of problem existed. Fourth, psychiatric confirmation of anxiety in those screened positive was not possible. Fifth, this study evaluated outcomes at discharge and conclusions about longer-term prevalence, spontaneous recovery, and chronic manifestations cannot be drawn at this point and further follow-up of the cohort is therefore indicated.
Recommendations
In spite of these limitations, there are several recommendations that can be derived from our findings. First, awareness about possible functional and psychological consequences of Covid-19 needs to be raised in patients, care providers, and health policy makers[28]. Second, early psychological[39] and pulmonary rehabilitation interventions including mobilization and exercise[27, 28, 40] which have been investigated in other health conditions could be effective and important to meet therapeutic windows; consequently their application (besides pharmacological treatments) needs to be studied in the population of Covid-19 patients, particularly those with severe disease. Third, this involves a better integration of primary care and psychiatric and rehabilitation services which is unfortunately lacking in many lower as well as higher resourced countries including China.[41] Fourth, systematic community follow up and two-way referral systems between community health centers and specialized hospital units, another weak point of many health systems, need to be strengthened.[41] Fifth, we recommend that future prospective cohort studies in Covid-19 survivors include functional and psychiatric outcomes.
Conclusion
A significant proportion of Covid-19 survivors had disability and anxiety at discharge from hospital. Disease severity was the only independent risk factor with consistent adverse effects on all outcomes. Health systems need to be prepared for an additional long-term burden due to Covid-19. This includes raising awareness about mental and physical health problems of survivors, early interventions, strengthening medical follow up, and increasing physical and psychological rehabilitation capacity.
Author contributions
CH, SZ and JR had full access to all of the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis. CH, SZ, QG, LY, QW, YY and JR designed the study. WX, XC, YH, DZ, YZ, GZ, SW, YW, ZZ, HL, CZ, BZ, JY, MF, ZN, BC, CD, and HH collected and processed the data. JR, SZ, QG, CH, and YY interpreted the data. JR performed the statistical analysis and wrote the first draft of the manuscript. All other authors revised the manuscript for critical content. All authors reviewed and approved the final submitted version.
Funding
This study was funded by the National Natural Science Foundation (81972146), the Department of Science and Technology of Sichuan Province (20YYJC3320), China Postdoctoral Science Foundation (2020M673251), Health Commission of Sichuan Province (20PJ034), and West China Hospital of Sichuan University (HX-2019-nCoV-011 to Chengqi He, and 2019HXBH058 to Siyi Zhu). The funders played no role in the design, conduct, or reporting of this study.
Competing interests
The authors declare that they have no competing interests.
Supporting information
Table S1. Outcome: One or more IADL problems (Lawton IADL scale, at least one item endorsed)&. & n= 431 for adjusted model and unadjusted model for age, n= 432 for all other models. RR = Risk Ratio. CI = Confidence Interval.
Table S2. Outcome: At least moderate ADL dependence (Barthel Index, < 75)&. &n=431 for adjusted model and unadjusted model for age, n= 432 for all other models. RR = Risk Ratio. CI = Confidence Interval.
Table S3. Outcome: Anxiety (Zung’s self-reported anxiety scale, raw score < 40)&. &n=407 for adjusted model and n=407 for unadjusted model for length of stay, n= 431 for unadjusted model for age, n= 432 for all other models. RR = Risk Ratio. CI = Confidence Interval.
Acknowledgment
We thank all the hospital staff members and coordinators under the CSPRM Covid-19 network for their efforts in collecting the information that was used in this study; Liangjiang Huang (Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology), Xu Qin (No.3 People’s Hospital of Chengdu), Minqing Li (Dazhou Central Hospital), Kehui Hu (Suining Central Hospital), and Jinling Zhang (West China Second University Hospital, Sichuan University) for their dedication to data entry and verification; Prof. Deying Kang (Epidemiologist, West China Hospital, Sichuan University) and Mrs. Rouyue Fang (Statistician, Industrial and Commercial Bank of China) for their assistance in initial data extraction and processing; all the patients who consented to provide their data for analysis and the medical staff members who are at the frontline caring for patients with Covid-19. We would like to thank Prof. Dr. Gerold Stucki and Cristiana Baffone for constructive comments on an earlier version of the manuscript.
Footnotes
We have discovered a mistake in the Figure 1, upper panel, adjusted results. Due to an ommission in the code, the figure called upon the wrong data and the adjusted results were not displayed correctly in the Figure as reported in the Appendix and manuscript. All data and parameter estimates in the tables and manuscript are correctly reported. The Figure 1 was replaced by the corrected version. Sincerely apologize for the mistake.