key: cord-282097-a1pwq4fi authors: Hoertel, N.; Sanchez, M.; Vernet, R.; Beeker, N.; Neuraz, A.; Alvarado, J.; Daniel, C.; Paris, N.; Gramfort, A.; Lemaitre, G.; Salamanca, E.; Bernaux, M.; Bellamine, A.; Burgun, A.; Limosin, F. title: Dexamethasone use and Mortality in Hospitalized Patients with Coronavirus Disease 2019: a Multicenter Retrospective Observational Study date: 2020-10-27 journal: nan DOI: 10.1101/2020.10.23.20218172 sha: doc_id: 282097 cord_uid: a1pwq4fi Objective: To examine the association between dexamethasone use and mortality among hospitalized patients for COVID-19. Design: Multicenter observational retrospective cohort study. Setting: Greater Paris University hospitals, France. Participants: 12,217 adults hospitalized with COVID-19 between 24 January and 20 May 2020, including 171 patients (1.4%) who received dexamethasone orally or by intravenous perfusion during the visit. Data source: Assistance Publique-Hopitaux de Paris Health Data Warehouse. Main outcome measures: The primary endpoint was time to death. We compared this endpoint between patients who received dexamethasone and those who did not in time-to-event analyses adjusting for sex, age, obesity, current smoking status, any medical condition associated with increased COVID-19-related mortality, and clinical and biological severity of COVID-19 at admission, while stratifying by the need of respiratory support (i.e., oxygen or intubation). The primary analysis was a multivariable Cox model and the secondary analysis used a univariate Cox regression in a matched analytic sample. Results: Among patients who required respiratory support, the end-point event of death occurred in 10 patients (15.9%) who received dexamethasone and 298 patients (26.4%) who did not. In this group of patients, there was a significant association between dexamethasone use and reduced mortality in both the crude, unadjusted analysis (hazard ratio (HR), 0.40; 95% CI, 0.18 to 0.87, p=0.021) and the adjusted multivariable analysis (HR, 0.46; 95% CI, 0.22 to 0.96, p=0.039). In the sensitivity analysis, the univariate Cox regression model in the matched analytic sample yielded a same tendency, albeit non-significant (HR, 0.31; 95% CI, 0.08 to 1.14, p=0.077). Among patients without respiratory support, the end-point event of death occurred in 14 patients (13.0%) who received dexamethasone and 1,086 patients (10.0%) who did not. In this group of patients, there was no significant association between dexamethasone use and the endpoint. When examining the association between the cumulative dose of dexamethasone received during the visit and the endpoint, we found that the administration of a cumulative dose between 60 mg to 150 mg among patients who required respiratory support was significantly associated with a lower risk of death in the crude, unadjusted analysis (HR, 0.28; SE, 0.58, p=0.028), the adjusted multivariable analysis (HR, 0.24; SE, 0.65, p=0.030), and in the univariate Cox regression model in the matched analytic sample (HR, 0.32; SE, 0.58, p=0.048), whereas no significant association was observed with a different dose. Among patients without respiratory support, there was no significant association between the cumulative dose of dexamethasone and the endpoint in the crude and in the adjusted multivariable analyses. Conclusions: In this observational study involving patients with Covid-19 who had been admitted to the hospital, dexamethasone use administered either orally or by intravenous injection at a cumulative dose between 60 mg and 150 mg was associated with decreased mortality among those requiring respiratory support. contributed to database build process. AN, NP, CD, AG, GL, MB, and AB critically revised the manuscript for scientific content. This work did not receive any external funding. . CC-BY-NC-ND 4.0 International license It is made available under a perpetuity. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted October 27, 2020. ; https://doi.org/10.1101/2020. 10 .23.20218172 doi: medRxiv preprint (10.0%) who did not. In this group of patients, there was no significant association between dexamethasone use and the endpoint. When examining the association between the cumulative dose of dexamethasone received during the visit and the endpoint, we found that the administration of a cumulative dose between 60 mg to 150 mg among patients who required respiratory support was significantly associated with a lower risk of death in the crude, unadjusted analysis (HR, 0.28; SE, 0.58, p=0.028), the adjusted multivariable analysis (HR, 0.24; SE, 0.65, p=0.030), and in the univariate Cox regression model in the matched analytic sample (HR, 0.32; SE, 0.58, p=0.048), whereas no significant association was observed with a different dose. Among patients without respiratory support, there was no significant association between the cumulative dose of dexamethasone and the endpoint in the crude and in the adjusted multivariable analyses. Conclusions: In this observational study involving patients with Covid-19 who had been admitted to the hospital, dexamethasone use administered either orally or by intravenous injection at a cumulative dose between 60 mg and 150 mg was associated with decreased mortality among those requiring respiratory support. Key words: Covid-19; SARS-CoV-2; dexamethasone; treatment; efficacy; death; oxygen; ventilation. . CC-BY-NC-ND 4.0 International license It is made available under a perpetuity. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted October 27, 2020. ; On June 16 th , an interim analysis of the RECOVERY (Randomized Evaluation of COVid-19 thERapY) trial, a randomized clinical trial examining a range of potential treatments for COVID-19, indicated that low-dose dexamethasone could reduce mortality in patients with COVID-19 requiring oxygen or mechanical ventilation support. 1 In that study, a total of 2,104 patients were randomized to receive dexamethasone 6 mg once per day for ten days, administered either orally or by intravenous injection, and were compared with 4,321 patients randomized to usual care alone. Dexamethasone was significantly associated with reduced 28day mortality in ventilated patients (rate ratio (RR), 0.65; 95% confidence interval (CI), 0.48 to 0.88; p=0.0003) and in patients receiving oxygen only (RR, 0.80; 95% CI, 0.67 to 0.96; p=0.0021). No benefit was observed among patients who did not require respiratory support (RR, 1.22; 95% CI, 0.86 to 1.75; p=0.14). These preliminary findings are of utmost importance and highlight that research into dexamethasone use in patients with COVID-19 is a priority. In this report, we present results of a large multicenter retrospective observational study of patients admitted for Covid-19 to any of the 39 Greater Paris University hospitals. We examined whether oral or intravenous administration of dexamethasone to hospitalized adult patients with COVID-19 was associated with reduced mortality (i) among those who required respiratory support, i.e. mechanical ventilation or oxygen, and (ii) in those who did not. Following results of the RECOVERY trial interim analysis, 1 we hypothesized that dexamethasone administration would be associated with reduced mortality in patients with COVID-19 who required respiratory support and not in those who did not. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted October 27, 2020. ; https://doi.org/10.1101/2020.10.23.20218172 doi: medRxiv preprint We conducted this study at AP-HP, which comprises 39 hospitals, 23 of which are acute, 20 adult and 3 pediatric hospitals. We included all adults aged 18 years or over who have been admitted with COVID-19 to these medical centers from the beginning of the epidemic in France, i.e. January 24 th , until May 20 th . For all patients, COVID-19 was ascertained by a positive reverse-transcriptase-polymerase-chain-reaction (RT-PCR) test from analysis of nasopharyngeal or oropharyngeal swab specimens. This observational study using routinely collected data received approval from the Institutional Review Board of the AP-HP clinical data warehouse (decision CSE-20-20_COVID19, IRB00011591). AP-HP clinical Data Warehouse initiative ensures patients' information and consent regarding the different approved studies through a transparency portal in accordance with European Regulation on data protection and authorization n°1980120 from National Commission for Information Technology and Civil Liberties (CNIL). We used data from the AP-HP Health Data Warehouse ('Entrepôt de Données de Santé (EDS)'). This warehouse contains all the clinical data available on all inpatient visits for COVID-19 to any of the 39 Greater Paris University hospitals. The data obtained included patients' demographic characteristics, vital signs, laboratory test and RT-PCR test results, medication administration data, current medication lists, current diagnoses, oxygen and ventilator use data, and death certificates. We obtained the following data for each patient at the time of the hospitalization: sex; age, which was categorized based on the OpenSAFELY study results 2 (i.e. 18-50, 51-70, 71+); obesity, defined as having a body-mass index higher than 30 kg/m 2 or an International is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted October 27, 2020. ; https://doi.org/10.1101/2020.10.23.20218172 doi: medRxiv preprint Statistical Classification of Diseases and Related Health Problems (ICD-10) diagnosis code for obesity (E66.0, E66.1, E66.2, E66.8, E66.9); self-reported current smoking status; any medical condition associated with increased COVID-19-related mortality 2,3 based on ICD-10 diagnosis codes, including diabetes mellitus (E11), diseases of the circulatory system (I00-I99), diseases of the respiratory system (J00-J99), neoplasms (C00-D49), and diseases of the blood and blood-forming organs and certain disorders involving the immune mechanism (D5-D8); clinical severity of COVID-19 at admission, defined as having at least one of the following criteria: 4 (i) respiratory rate > 24 breaths/min or < 12 breaths/min, (ii) resting peripheral capillary oxygen saturation in ambient air < 90% , (iii) temperature > 40°C, or (iv) systolic blood pressure < 100 mm Hg; and biological severity of COVID-19 at admission, defined as having at least one of the following criteria: 4,5 (i) high neutrophil-to-lymphocyte ratio or (ii) low lymphocyte-to-C-reactive protein ratio (both variables were dichotomized at the median of the values observed in the full sample), or (iii) plasma lactate levels higher than 2 mmol/L. All medical notes and prescriptions are computerized in Greater Paris University hospitals. Medications and their mode of administration (i.e., dosage, frequency, date, condition of intake) were identified from medication administration data or scanned handwritten medical prescriptions, through two deep learning models based on BERT contextual embeddings, 6 one for the medications and another for their mode of administration. The model was trained on the APmed corpus, 7 a previously annotated dataset for this task. Extracted medications names were then normalized to the Anatomical Therapeutic Chemical (ATC) terminology using approximate string matching. The endpoint was the time from study baseline to death. Patients without an end-point event had their data censored on May 20 th , 2020. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted October 27, 2020. ; https://doi.org/10.1101/2020.10.23.20218172 doi: medRxiv preprint Study baseline was defined as the date of hospital admission. Dexamethasone use was defined as receiving this medication orally or by intravenous injection at any time during the follow-up period, from study baseline to the end of the hospitalization or death. Dexamethasone cumulative dose received was calculated and considered as a categorical variable with the following categories defined a priori: (i) 60 mg to 150 mg based on usual prescribing practice for acute respiratory distress syndrome (ARDS) in AP-HP hospitals (corresponding to 10 mg per day for 6 days to 20 mg per day for 5 days followed by 10 mg per day for 5 days), (ii) other doses (i.e., more than 150 mg or less than 60 mg), and (iii) missing data. All analyses were stratified by the need for respiratory support, i.e. oxygen or mechanical ventilation, at any time during the follow-up. We calculated frequencies and means (± standard deviations (SD)) of each baseline characteristic described above in patients receiving and not receiving dexamethasone, and compared them using chi-square tests or Welch's t-tests. To examine the association of dexamethasone use and the endpoint, we performed Cox proportional-hazards regression models. To help account for the nonrandomized prescription of dexamethasone and reduce the effects of confounding, the primary analysis used a multivariable Cox regression model including as covariates sex, age, obesity, current smoking status, any medical condition associated with increased COVID-19-related mortality, and clinical and biological severity of COVID-19 at admission. Weighted Cox regression models were used when the proportional hazards assumption was not met. Kaplan-Meier is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted October 27, 2020. ; https://doi.org/10.1101/2020.10.23.20218172 doi: medRxiv preprint curves were performed and their pointwise 95% confidence intervals were estimated using the nonparametric bootstrap method. 8 As a sensitivity analysis, we performed a univariate Cox regression model in a matched analytic sample using a 1:10 ratio, based on the same variables used for the multivariable Cox regression analysis. To reduce the effects of confounding, optimal matching was used in order to obtain the smallest average absolute distance across all clinical characteristics between exposed patient and non-exposed matched controls. We also examined whether the cumulative dose of dexamethasone received during the visit was associated with the endpoint. Finally, we performed multivariable Cox regression models including interaction terms to examine whether the association between dexamethasone use and the primary endpoint significantly differed across baseline characteristics. For all associations, we performed residual analyses to assess the fit of the data, check assumptions, including proportional hazards assumptions, and examined the potential influence of outliers. To improve the quality of result reporting, we followed the recommendations of The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) Initiative. 21 Statistical significance was fixed a priori at two-sided p-value<0.05. All analyses were conducted in R software version 2.4.3 (R Project for Statistical Computing). Of the 16,170 with a positive COVID-19 RT-PCR test consecutively admitted to the hospital from January 24 th to May 20 th , a total of 3,953 patients (24.4%) were excluded because of missing data or their age (i.e. less than 18 years old of age). Of the remaining is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint were excluded because the route of administration was either ophthalmic, through aerosol, or unknown. Among the 171 remaining patients exposed to dexamethasone, 63 (36.8%) needed respiratory support (i.e. mechanical ventilation or oxygen) and 108 (63.2%) did not (Figure 1 ). The mean cumulative dose administered was 107.8 mg (SD=63.5; median=100 mg; range: 10.0 mg to 320.0 mg) in those who required respiratory support and 101.8 mg (SD=89.8; median=90 mg; range: 10.0 mg to 600.0 mg) in those who did not. This treatment was administered orally in 97.1% of patients and by intravenous injection in 2.9% of them. COVID-19 RT-PCR test results were obtained after a mean delay of 3.7 days (SD=8.4, median=0.9 days) from the date of hospital admission in patients who required respiratory support. This delay was not significantly different between patients receiving or not receiving dexamethasone [mean delay in the exposed group=3.0 days (SD=6.2); mean delay in the nonexposed group=3.7 days (SD=8.5); Welch's t-test=0.79, p=0.430)]. In patients who did not require respiratory support, the mean delay was 4.7 days (SD=10.2, median=1.0 day), and this delay was not significantly different between patients receiving or not receiving dexamethasone [mean delay in the exposed group=3.5 days (SD=8.9); mean delay in the nonexposed group=4.8 day (SD=10.2); Welch's t-test=1.5, p=0.163)]. Among patients who required respiratory support, the mean follow-up was 27.9 days (SD=20.4; median=21 days; range: 1 day to 106 days) and 308 patients (25.8%) had an endpoint event at the time of data cutoff on May 20 th . Among those who did not require respiratory support, the mean follow-up was 18.5 days (SD=24.6; median=6; range: 1 day to 117 days), and 1,100 (10.0%) patients had an end-point event at the time of data cutoff. Associations between baseline characteristics and the endpoint are given in eTable 1. The distribution of the patients' characteristics by dexamethasone use is shown in Table 1 . Dexamethasone use significantly differed in clinical and biological severity at admission is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted October 27, 2020. ; https://doi.org/10.1101/2020.10.23.20218172 doi: medRxiv preprint among patients who required respiratory support, and in sex, age, obesity and biological severity among those who did not. In the matched analytic samples, there were no significant differences in patients' characteristics according to dexamethasone use (Table 1) . Among patients who required respiratory support, the end-point event of death occurred in 10 patients (15.9%) who received dexamethasone and 298 patients (26.4%) who did not ( Table 2 ). In this group, we found a significant association between dexamethasone use and reduced risk of the endpoint of death in both the crude, unadjusted analysis (hazard ratio (HR), 0.40; 95% CI, 0.18 to 0.87, p=0.021) and the adjusted multivariable analysis (HR, 0.46; 95% CI, 0.22 to 0.96, p=0.039) (Figure 2 ; Table 2 ). In the sensitivity analysis, the univariate Cox regression model in the matched analytic sample yielded a same tendency, albeit nonsignificant (HR, 0.31; 95% CI, 0.08 to 1.14, p=0.077) (Figure 3 ; Table 2 ). Among patients without respiratory support, the end-point event of death occurred in 14 patients (13.0%) who received dexamethasone and 1,086 patients (10.0%) who did not ( Table 2 ). In this group of patients, there was no significant association between dexamethasone use and the endpoint, neither in the crude, unadjusted analysis (HR, 0.73; 95% CI, 0.42 to 1.26, p=0.253) or in the adjusted multivariable analysis (HR, 0.59; 95% CI, 0.30 to 1.16, p=0.126) ( Figure 2 ; Table 2 ). In the sensitivity analysis, the univariate Cox regression model in the matched analytic sample yielded a similar result (HR, 1.06; 95% CI, 0.45 to 2.49, p=0.894) (Figure 3; Table 2) . A post-hoc analysis indicated that in the full sample of patients with respiratory support, we had 80% power to detect hazard ratios for dexamethasone treatment of at least 2.33/0.15 and of at least 2.37/0.14 in the matched analytic sample. In those without respiratory support, we had 80% power to detect hazard ratios for dexamethasone treatment of at least 1.99/0.30 in the full sample and 2.05/0.29 in the matched analytic sample. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted October 27, 2020. ; https://doi.org/10.1101/2020. 10.23.20218172 doi: medRxiv preprint When examining the association between the cumulative dose of dexamethasone received during the visit and the endpoint, we found that the administration of a cumulative dose between 60 mg to 150 mg among patients who required respiratory support was significantly associated with a lower risk of death in the crude, unadjusted analysis (HR, 0.28; SE, 0.58, Finally, the association between dexamethasone use and the outcome did not significantly differ across subgroups defined by baseline characteristics in both groups with and without respiratory support, except for patients with obesity who did not require respiratory support, for whose dexamethasone use was significantly associated with a higher risk of death compared to obese patients without dexamethasone (HR, 3.90; 95% CI, 1.13 to 13.44, p=0.031) (eTable 3). However, none of the 16 patients with obesity who received a cumulative dose of 60 to 150 mg of dexamethasone died. In this multicenter retrospective observational study involving a large sample of patients admitted to the hospital with COVID-19, we found that dexamethasone use, administered either orally or by intravenous injection at a cumulative dose between 60 mg and . CC-BY-NC-ND 4.0 International license It is made available under a perpetuity. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted October 27, 2020. ; https://doi.org/10.1101/2020.10.23.20218172 doi: medRxiv preprint 150 mg, was significantly and substantially associated with reduced mortality among patients with COVID-19 requiring oxygen or mechanical ventilation support. This association did not significantly differ according to baseline clinical characteristics. No significant association between dexamethasone use and mortality was observed among patients with COVID-19 without respiratory support, except in the univariate Cox regression model in the matched analytic sample where the administration of a cumulative dose between 60 mg to 150 mg was significantly associated with reduced mortality. Although these findings should be interpreted with caution due to the observational design, they are in line with the results of the RECOVERY trial interim analysis, 1 which indicated that dexamethasone 6 mg once per day for ten days, administered either orally or by intravenous injection, was significantly associated with reduced 28-day mortality in ventilated patients and in patients receiving oxygen only, and extend them by suggesting that cumulative dose between 60 mg and 150 mg might be efficient in reducing mortality among patients with COVID-19 requiring respiratory support. Our study has several limitations. First, some amount of unmeasured confounding may remain. However, our analyses adjusted for numerous potential confounders, including sex, age, obesity, current smoking status, any medical condition associated with increased COVID-19-related mortality, and clinical and biological severity of COVID-19 at admission. Second, there are missing data for some variables and potential for inaccuracies in the electronic health records, such as the possible lack of documentation of illnesses or medications, or the misidentification of treatment mode of administration (e.g., dose, frequency), especially for hand-written medical prescriptions. Finally, despite the multicenter design, our results may not be generalizable to other settings or regions. In this observational study involving patients with Covid-19 who had been admitted to the hospital, dexamethasone use administered either orally or by intravenous injection at a . CC-BY-NC-ND 4.0 International license It is made available under a perpetuity. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted October 27, 2020. ; https://doi.org/10.1101/2020.10.23.20218172 doi: medRxiv preprint cumulative dose between 60 mg and 150 mg was associated with decreased mortality among patients requiring respiratory support. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted October 27, 2020. ; https://doi.org/10.1101/2020.10.23.20218172 doi: medRxiv preprint is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted October 27, 2020. ; https://doi.org/10.1101/2020.10.23.20218172 doi: medRxiv preprint Figure 3 . Kaplan-Meier curves for death in the matched analytic samples of hospitalized patients with Covid-19 who required respiratory support (i.e., mechanical ventilation or oxygen) (N=693) (A) and of those who did not (N=1,188) (B) according to dexamethasone use. The shaded areas represent pointwise 95% confidence intervals. For each exposed case, ten non-exposed controls were selected. . CC-BY-NC-ND 4.0 International license It is made available under a perpetuity. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted October 27, 2020. ; https://doi.org/10.1101/2020.10.23.20218172 doi: medRxiv preprint Table 1 . Characteristics of patients receiving or not receiving dexamethasone in patients with or without respiratory support (oxygen or intubation). Exposed to dexamethasone Not exposed to dexamethasone Non-exposed matched group Exposed to dexamethasone vs. ). β Current smoking status was self-reported. γ Assessed using ICD-10 diagnosis codes for diabetes mellitus (E11), diseases of the circulatory system (I00-I99), diseases of the respiratory system (J00-J99), neoplasms (C00-D49), and diseases of the blood and blood-forming organs and certain disorders involving the immune mechanism (D5-D8) µ Defined as having at least one of the following criteria: respiratory rate > 24 breaths/min or < 12 breaths/min, resting peripheral capillary oxygen saturation in ambient air < 90%, temperature > 40°C, or systolic blood pressure < 100 mm Hg. κ Defined as having at least one of the following criteria: high neutrophil-to-lymphocyte ratio, low lymphocyte-to-C-reactive protein (both variables were dichotomized at the median of the values observed in the full sample), and plasma lactate levels higher than 2 mmol/L. * p-value is significant (p<0.05). 12,217 adult inpatients (178 exposed to dexamethasone and 12,039 not exposed) included in the analyses 12,210 adult inpatients (171 exposed to dexamethasone administered either orally or by intravenous injection and 12,039 not exposed) included in the analyses 7 patients with dexamethasone were excluded due to the route of administration: -Aerosol: n = 2 -Ophthalmology: n = 1 -Unknown: n = 4 1,192 with respiratory support (63 exposed to dexamethasone) 11 ,018 without respiratory support (108 exposed to dexamethasone) A B . CC-BY-NC-ND 4.0 International license It is made available under a perpetuity. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted October 27, 2020. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted October 27, 2020. ; https://doi.org/10.1101/2020.10.23.20218172 doi: medRxiv preprint News & events. Dexamethasone reduces death in hospitalised patients with severe respiratory complications of COVID-19 OpenSAFELY: factors associated with COVID-19-related hospital death in the linked electronic health records of 17 million adult NHS patients Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study. The Lancet Statement on the management at home or in a care facility of suspected or confirmed Covid-19 patients Neutrophil-to-lymphocyte ratio and lymphocyte-to-Creactive protein ratio in patients with severe coronavirus disease 2019 (COVID-19): A meta-analysis Pre-training of deep bidirectional transformers for language understanding MedExt: combining expert knowledge and deep learning for medication extraction from French clinical texts Nonparametric standard errors and confidence intervals The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement: guidelines for reporting observational studies. Annals of internal medicine Marginal structural models and causal inference in epidemiology Observational study of hydroxychloroquine in hospitalized patients with Covid-19