Abstract
Objectives The interest in interleukin-6 receptor antagonists (IL-6RA) and steroids have increased recently due to their potential role as immunomodulatory effect in critically ill coronavirus disease (COVID-19). The magnitude of this therapy in subgroups of patients with invasive mechanical ventilation (MV) remains to be fully clarified. We compared the clinical characteristics and outcomes of patients requiring iMV, and receiving IL-6RA and steroids with different steroids regimens.
Design International, multicenter, observational study derived from Viral Infection and Respiratory Illness University Study registry and conducted through Discovery Network, Society of Critical Care Medicine. Marginal structural modeling was used to adjust time-dependent confounders; observations were weighted using inverse probability of treatment weight. A sensitivity analysis was conducted for target trial design.
Setting 168 hospitals, 16 countries.
Patients Covid-19 ICU patients (≥18 years) requiring MV between March 01,2020, and January 10,2021.
Intervention None.
Measurements and Main Results Of 860 patients met eligibility criteria, 589 received steroids, 170 IL-6RA, and 101 combinations; groups were balanced after adjustment. Median daily steroid dose was 7.5 mg dexamethasone or equivalent (IQR:6-14 mg); 80.8% and 19.2% received low-dose and high-dose steroids, respectively. The median C-reactive protein level was >75 mg/L in majority of our cohort. The use of IL-6R antagonists alone or in combination was not associated with a significant difference in ventilator-free days (VFD) compared to steroids alone with different steroids regimens (adjusted incidence rate ratio [95% CI]): IL-6R antagonists (1.12 [0.88,1.4]), combination (0.83 [0.6,1.14]). Patients treated with low or high-dose steroids had non-significant differences in VFD compared to IL-6RA (ß=0.62, 95% CI −1.54,2.78 for low-dose steroid; ß=-1.19, 95% CI −3.85,1.47 for high-dose steroid). There was no difference in 28-day mortality and hospital mortality with IL-6RA alone or in combination compared to steroids alone (28-day mortality adjusted odds ratio [95% CI]): IL-6RA (0.68[0.44,1.07]), combination (1.07[0.67,1.70]). Sensitivity analysis findings were consistent with primary analysis. Liver dysfunction was higher in IL-6RA (p=0.04) while rate of bacteremia did not differ among groups.
Conclusions In adult ICU COVID-19 patients on iMV, we found no difference in outcomes between those who received IL-6RA, steroids, or combination therapy and those who received IL-6RA or low-or high-dose steroids. Further randomized trials are needed to enhance our understanding for IL-6RA safety with different steroids regimen and the magnitude of benefit in those subgroups of patients.
Background
Severe coronavirus disease 2019 (COVID-19) can manifest as respiratory failure with elevated inflammatory markers, resulting in cytokine storms, for which interleukin-6 receptor antagonists (IL-6RA) are approved treatments (1–3). Interest in IL-6RA and corticosteroids has increased recently due to their potential role as immunomodulators (4–7). In view of the results from International Randomized, Embedded, Multi-factorial, Adaptive Platform Trial for Community-Acquired Pneumonia (REMAP-CAP) and Randomized Evaluation of COVID-19 Therapy (RECOVERY), guidelines by National Institutes of Health and Infectious Disease Society of America conditionally suggest tocilizumab in combination with steroids (low-dose dexamethasone, 6 mg daily for 10 days) for ICU patients exhibiting rapid respiratory failure progression or high inflammatory markers (8–11). Data with other steroid regimens are limited (high-dose dexamethasone, 20 mg IV daily) which are commonly used for non-COVID-19 acute respiratory distress syndrome (ARDS) and may be utilized for COVID-19 ARDS (12, 13). Moreover, a limited number of patients on invasive mechanical ventilation (MV) were included in the previous trials (29 % in REMAP-CAP, 14 % in RECOVERY, 16 % in TOCIBRAS, and 37% in COVACTA) (13). Therefore, the magnitude of this combination therapy in subgroups of patients with iMV remains to be fully clarified and was identified as a research priority to enhance personalized clinical decision-making (14, 15). We sought to describe clinical characteristics and outcomes of ICU patients requiring iMV, and receiving IL-6RA and steroids using the data from the largest COVID-19 registry and provided further insights into sarilumab and utilization pattern of multiple doses of steroids.
Materials and Methods
Data Source
Data for this ancillary study were derived from Viral Infection and Respiratory Illness Universal Study (VIRUS) registry, an international, multicenter, observational study conducted through Discovery Network, Society of Critical Care Medicine (SCCM). The study was approved by SCCM Scientific Review Committee, SCCM Discovery Steering Committee, and King Faisal Specialist Hospital and Research Center (KFSH&RC) Institutional Review Board (IRB) (IRB# 2201053, registered at ClinicalTrials.gov NCT04486521). In some participating hospitals, verbal consent was obtained whereas in other countries consent was waived by local research ethics committees. Local investigators were responsible for obtaining local approval in line with applicable regulations. VIRUS is one of largest registries that consecutively collects data on COVID-19 patients. Details of this database described elsewhere (16–19) (Supplemental Table 1). The study is reported following the STrengthening the Reporting of Observational Studies in Epidemiology (STROBE) and the Risk Of Bias In Nonrandomized Studies of Interventions (ROBINS-I) guidelines (16–19).
Study Participants and Exposures
Adult patients (≥ 18 years) were eligible if admitted to ICU from March 01, 2020 to January 10, 2021, required invasive mechanical ventilation (iMV), had positive polymerase chain reaction (PCR) SARS-CoV-2, and received IL-6RA (tocilizumab or sarilumab), corticosteroid (dexamethasone, hydrocortisone, methylprednisolone, or prednisone), or combination. We excluded patients who were < 17 years, repeatedly admitted to ICU during the same hospital visit, on chronic systemic steroid at home or taking steroids for indications other than COVID-19, or died before receiving IL-6RA or steroid. We stratified population into three groups: IL-6RA, corticosteroid, and combination. Corticosteroid group was stratified further based on dexamethasone equivalent (mg) into high and low-dose. High-dose was defined as > 15 mg/day of dexamethasone or equivalent (12, 13). Patients were followed until hospital discharge, death, 28 days, or January 13th, 2021-date on which database was locked for analysis.
Data Collection
The following variables were collected on days 0 (ICU admission), 1-3, 7, 14, 21, and 28 from electronic medical records according to standard operating procedure: demographics, biomarkers (ferritin, interleukin-6 (IL-6), D-dimer, fibrinogen, C-reactive protein [CRP], lactate dehydrogenase, and lymphocyte count), microbiology, concomitant medications, MV duration, discharge status, PaO2(mm Hg)/FiO2 (arterial partial pressure of oxygen over fractional inspired oxygen concentration; PF ratio), and timing of drug initiation and dose. Data were collected online and stored on a secure data server using Mayo Clinic Research Electronic Data Capture application (REDCap). Data quality was assessed routinely to ensure completeness and accuracy.
Study Outcomes and Variable Definitions
Primary outcome was ventilator-free days score (VFD) from iMV at day 28 which was chosen as a patient-centered outcome and highly influenced by mortality (20). Secondary outcomes included: ICU and hospital mortality, 28-day mortality, hospital and ICU length of stay (LOS), median change in PF ratio, adverse events (AEs), and biomarkers concentration. Details on outcome definitions are provided in Supplemental Table 2.
Statistical Analysis
Sample size was determined pragmatically, based on all available ICU patients in VIRUS database who met eligibility criteria. Statistical analyses were performed using R software, V3.6.3 (Vienna, Austria). Counts and percentages were used to represent categorical variables. Continuous data were summarized using means ± standard deviations (SDs) or medians (interquartile ranges). Chi-square test was used to compare distribution of categorical variables and Kruskal–Wallis test and ANOVA were used to compare distributions of continuous non-normal and normal variables, respectively.
We evaluate three statistical models for this study. Model one was for Marginal structural models (MSM) in the primary analysis which included the overall treatment strategy not solely restricted to the first 48 hours. Model two was exploratory analyses with added Sequential Organ Failure Assessment (SOFA) score, and highest FiO2 to the base model. Model three is a sensitivity analysis to emulate target trial approach restricted to the treatment received in the first 48 hours.
Model 1
Marginal structural models (MSM) were used to compare three regimens after adjusting for non-time and time-varying covariates (21). Observations were weighted by inverse probability of treatment weight (IPTW). Nonparametric modeling using generalized boosted model (GBM) estimate IPTWs with automatic handling of missing data (22). TWANG package in R was used for analysis. Covariates distributions were compared before and after applying propensity weights. Standardized mean difference (SMD) was used to assess balance of pseudo-population using Kolmogorov–Smirnov stop rule. Average treatment effect (ATE) estimation method was used throughout analysis on entire sample. Covariates included in model were judged as likely to influence outcomes and have been identified in several studies: age; sex; ethnicity; asthma/chronic obstructive pulmonary disease (COPD); ARDS grade; history of diabetes, hypertension, coronary heart disease, or congestive heart failure; the lowest FiO2; therapeutic anticoagulation; hydroxychloroquine; azithromycin; and antivirals (including remdesivir); vasopressors; and paralytics (Supplemental File 1). Analyses were performed using robust variance estimators taking into account clustering within hospitals and patients. Quasi-Poisson generalized linear modeling was used to compare VFD after weighting. Patients who died before 28-day were assigned zero VFD. Quasi-Poisson regression coefficients were exponentiated to obtain incidence rate ratio (IRR), which is expected change in VF-days (ratio), compared to reference category. Quasibinomial generalized linear modeling was used to compare mortality between three groups. Weighted Kaplan–Maier estimator was used to compare hospital and ICU LOS. Patients who were still in ICU after 28-day study period were censored at 28 days. Log-rank test was used to compare survival before and after weighing. Cause-Specific Cox regression analysis was used to assess factors associated with mortality at 28 days. Discharge and mortality were included as competing risks. Linear mixed modeling was used to compare change in laboratory values and PF ratio. Hypothesis testing was performed at 5% level of significance.
Model 2
In exploratory analyses of potential factors associated with variation in VFD and mortality, we added Sequential Organ Failure Assessment (SOFA) score, and highest FiO2 to base model 1.
Model 3
We performed a sensitivity analysis to emulate target trial to reduce immortal time bias. Eligible patients for target trial approach were those on MV, received IL-6RA or steroids within ICU days 0 or 1. Patients who died within first 48 hours were excluded. Those who get treatment of interest after day 1 were classified by treatment exposure in first 48 hours (similar to intention-to-treat analysis in randomized trial) (23). The outcomes of interest in target trial approach were VFD and 28-day mortality.
Results
Overall Characteristics
A total of 23,783 patients were screened, 860 met eligibility criteria and classified as follows: 170 received IL-6RA, 589 received steroids, and 101 received both therapies. In sensitivity analysis to emulate target trial, 562 patients fulfilled inclusion criteria: 406 received steroids (72.2%), 121 received IL6-RA (21.5%), and 35 received both (6.23%) (Figure 1). In unadjusted analysis, baseline characteristics were not balanced. IL-6RA patients were younger (median 57.5 years vs. 61.1 in steroid vs. 61.8 in combination; p=0.009) and chronic pulmonary diseases, excluding asthma (COPD, bronchiectasis, and interstitial lung disease), were more prevalent in steroid group (59 [10.0%] vs. 7 [4.12%] in IL-6RA vs. 16 [15.8%] in combination, p=0.005). The ARDS was prevalent in 35.5% of steroid-treated patients, 14.7% in the IL-6R antagonist, and 36.6 % in combination therapy; the majority of these patients had severe ARDS grades (PF ratio < 100). There was no significant difference in the incidence of acute liver failure, and bacteremia among the three arms. After adjustments for baseline covariates, SMD was ∼0.1 for majority of the covariates, indicating three groups were more likely balanced (Table 1).
The median daily steroid dose was 7.5 mg dexamethasone equivalent (IQR 6-14 mg). Of these, 80.8% received low-dose steroids, and 19.2% received high-dose steroids. Of the patients who received IL-6RA, the majority (84.5%) received only one dose of tocilizumab while second and third doses were administered to 13.1% and 2.3%, respectively. Similarly, most patients (84.7%) in the sarilumab group received a single dose. The dosage range of tocilizumab was 600–800 mg and sarilumab was 200-400 mg (Table 2). Matched patients received comparable amounts of sedatives, paralytics, and anticoagulation. The CRP level (assessed in 413) was consistently > 75 mg/L across all groups while IL-6 concentrations (assessed in 24) were highest in combination group. Additional comorbidities, ventilator and radiology characteristics, ICU supports, and biomarker levels are included in Supplemental Tables 3-11.
Primary Outcomes
The results for outcomes before adjusting for covariates is presented in Supplemental Tables 12. After adjustment (model 1), we found no significant difference in VFDs between the use of IL-6R antagonists alone or in combination compared to using steroids alone with different steroids regimen (adjusted IRR [95% CI]): IL-6RA (1.12 [0.88,1.40]), combination (0.83 [0.6,1.14]). Exploratory analyses findings (model 2) were comparable to primary analysis. Sensitivity analysis as target trial design (model 3) showed consistent results when compared to steroids alone (aIRR [95% CI]): IL-6RA (1.13 [0.87,1.45]), combination (0.64 [0.34,1.23]) [Figure 2A]. Linear regression analysis was performed stratified by steroid dose and showed no significant differences between patients received low and high steroid doses compared to IL-6RA (ß = 0.62, 95% CI −1.54,2.78; p=0.57 for low-dose steroid; ß = −1.19, 95% CI −3.85,1.47; p=0.38 for high-dose steroid). Factors associated with higher odds of VFD score were non-use of paralytics, the use of therapeutic anticoagulation, and low FiO2 (Table 3).
Secondary Outcomes
There was no significant difference in 28-day mortality between IL-6RA alone or in combination compared to steroids alone with different steroids regimen in model 1 (aOR [95% CI]): IL-6RA (0.68 [0.44,1.07]), combination (1.07 [0.67,1.70]). Similarly, exploratory analyses (model 2), sensitivity analysis findings as target trial approach (model 3), and hospital mortality were consistent with primary analysis when compared to standalone steroids (hospital mortality aOR [95% CI]: IL6-RA (0.68, [0.43,1.09]), combination (1.23, [0.72,2.11]) [Figure 2B]. Cause-specific Cox proportional hazard regression showed that the hazard of mortality was not significantly different between IL-6RA (HR = 0.9, p>0.05) or combination (HR = 0.84, p>0.05) compared to steroids (Figure 3). Detailed analysis of secondary outcomes is available in Supplemental Tables 13-17, Figures 1-3. Binary logistic regression was used to assess 28-day-mortality stratified by steroid dose and showed no significant association between use of steroids (low-or high-dose) and mortality compared to IL-6RA (alone or combined). Weighted Kaplan-Maier estimates showed no statistically significant differences between the groups for ICU and hospital LOS (Log-rank p=0.52 and 0.18, respectively). Changes in PF ratio, and incidence of bacteremia were not significantly different as well. IL-6RA patients were more likely to experience liver dysfunction (23.9%) than steroids (14.6 %) or combination (11.9%) (p=0.04) while high troponin levels were reported in combination therapy.
Discussion
The IL-6RAs were previously investigated with conflicting results (24). Some trials found lower duration of ICU and hospital stays, lower MV and death composite rate, and increased number of organ support-free days in COVID-19 patients requiring respiratory support (8,9,25–29). A recent meta-analysis, trial sequential analysis (TSA), and meta-regression of RCTs by Snow TAC et al showed that the use of tocilizumab may be associated with a short-term mortality benefit in patients with COVID-19: OR 0.87 [0.74–1.01]; p=0.07; TSA adjusted CI 0.66–1.14 (30). There is an ongoing need exists to better identify subgroups of patients who will benefit the most from immunomodulation therapy according to the degree of respiratory support. Herein, we reported characteristics and outcomes of adult ICU COVID-19 patients, who required iMV, and received IL-6RA or steroids within large, multinational registry. We observed no significant difference in outcomes among IL6RA, steroids, and combination therapy and among those who received IL-6R antagonists compared to low-and high-dose steroids.
We proposed several reasons may explain our findings. First, the information on IL6-RA monotherapy use was collected prior to June, 2020 before announcement of RECOVERY steroids domain results. Before this time, we observed steroids use was limited due to concern of delayed viral clearance as shown in MERS coronavirus (31). This may explain the findings and point estimates that were seen here with standalone IL6-RA arm. Whether the effect of IL6-RA monotherapy is beneficial remains unclear and unlikely to be investigated in adequately powered RCT because the equipoise for withholding corticosteroids is no longer present given the convincing evidence for steroids which became standard care for ICU COVID-19 patients requiring MV (32).
Second, and most importantly, our study included only patients with iMV while RECOVERY or REMAP-CAP included varying degrees of respiratory support (iMV at baseline constitute 29 % in REMAP-CAP, and 14 % in RECOVERY) [8,9]. Therefore, the beneficial effect observed in RECOVERY or REMAP-CAP with combined IL6-RA and steroid therapy may be due to patients’ enrollment early in their critical care trajectory while the patients in our study at more severe end of the spectrum. This is in line with the recent Bayesian reanalysis of RECOVERY trial which showed COVID-19 hospitalized patients on non-invasive ventilation have a high probability of a clinically meaningful outcome benefit from tocilizumab and uncertainty for its effect in patients on iMV (15). Taken all together, we think the immunomodulatory effect of IL6-RA and steroid appears to be most beneficial for patients on high-flow nasal cannula, or non-invasive ventilation soon after the time of clinical deterioration at systemic hyperinflammation onset and before intubation or multi-system organ failure (14, 15).
Third, different types and doses of steroids including high-dose steroids in our analysis which could have muted the beneficial effect that IL6-RA might otherwise have had. The REMAP-CAP included fixed duration steroids and shock dependent hydrocortisone and RECOVERY used low-dose dexamethasone 6 mg IV daily. The Italian National Institute for Infectious Diseases recommends 10-day regimen: 5-day of methylprednisolone 1 mg/kg or dexamethasone 20 mg daily (defined as high-dose steroids in our study) (33). This regimen is higher than fixed dexamethasone dose in RECOVERY and that described in non-COVID ARDS literature (12, 13). It was also noted in 19.2% of our study cohort. Prior studies suggested moderate-to-high dose steroid regimens resulted in greater reduction in mortality, organ dysfunction, MV requirement, and no increase in medical or infectious complications compared to low-dose regimens in ICU COVID-19 patients (34, 35). Whether the beneficial immunomodulatory effects of IL-6R blockade could be achieved more easily and with less cost by using a different steroid regimen (higher dose) and its safety need to be investigated in adequately powered RCT (14). Currently, there is an ongoing trial for higher vs. lower doses of dexamethasone for COVID-19 and severe hypoxia (COVIDSTEROID2, NCT04509973) which will include IL6-RA subgroup analysis to provide more information on whether the dose of steroid will have an impact on clinical outcomes (36).
Lastly, racial and ethnic differences in the populations studied, and differences in the study conduction time may provide additional explanation given that the decline in in-hospital mortality from Covid-19 occurred in later phase of the pandemic (24).
Strengths of this study lie in multinational nature. We used data from the largest COVID-19 registry and represented spectrum of intensive care with racially and ethnically diverse cohorts, thereby maximizing generalizability. Moreover, we included large numbers of patients on iMV at baseline which was a limited subset of the population included in the previous trials. Additionally, we included laboratory-confirmed COVID-19 PCR tests, thereby minimizing selection or surveillance bias at each center. Lastly, prior observational studies could be biased by immortal time and indication bias, leading to possible overestimation of positive treatment effect (physicians were more likely to use IL6RA in patients who were more likely to survive or sicker as a last-resort salvage attempt) (14, 38). We used a novel statistical analysis to overcome these limitations and application of target trial design which is likely better approach to reduce immortal time bias and best resembling clinical practice (21–23). Although COVID-19 treatments were not used uniformly before matching and patients did not receive up-to-date standards of care, especially in early pandemic phases, we utilized rigorous methods to match the three groups concerning anticoagulation use, hydroxychloroquine, antivirals, and azithromycin. Importantly, this study reinforces risk of liver injury with IL-6 inhibition. In REMAP-CAP, safety assessment was limited by low AEs number. Clinicians should be aware of these effects and should weigh treatment risks and benefits. Whether this risk is intensified with repeated IL-6RA administration or related to patients’ critical presentation warrants investigation in future studies.
We acknowledge several limitations. First, due to the observational nature of this study, these results should be understood as possibly inconclusive due to insufficient power and frequentist statistical framework. Therefore, there is a possibility of type II errors. Adequately powered RCTs are still needed to further explore the causal treatment effect and clinical outcomes in ICU patients with iMV. Second, fair number of patients were excluded due to missing outcomes and incomplete data for time-dependent confounders to be used in MSM. It is possible that we may have missed patients that would have met eligibility criteria but were not included because of those reasons (possible elimination for non-random missingness). Although we used best available methods to compare well-balanced groups, controlling for confounders in observational study may remain incomplete despite all efforts. Notably, some imbalances in potential confounders (e.g., sex, cardiovascular disease, treatments offered) were observed after weighting. Such potential residual confounding may have a role in the observed non-difference between treatment strategies. Moreover, there was large proportion of incomplete data for SOFA and APACHE II scores, because of heavy burden experienced by participating clinicians during pandemic. We also did not report steroid and IL-6RA administration time from symptom onset and did not address potential role of IL-6 level in predicting IL-6RA responses, all-important questions warranting further investigation. Steroid duration was also highly variable with inconsistent data entry, making it impossible to fully minimize its effect. Finally, while our follow-up is limited to 28 days, additional follow-up may help characterize long-term sequelae, especially mucormycosis, aspergillosis, pneumocystis pneumonia, and multi-drug resistant organisms. Future adequately powered RCT should focus on assessing the survival as primary endpoint which is important to reflect patients, family members, clinicians, and regulators perspectives. Treating steroid dose as a continuous variable and its impact on 28-day mortality will be explored in the subsequent iterations of the registry through ancillary studies by conducting a meta-regression analysis. Going forward, we propose to include research findings, including ours, into the VIRUS web-based decision aid to evaluate risk, benefit, and burden of specific interventions considering prior knowledge, clinical context, preference, and how aforementioned parameters are modified by new findings.
Conclusions
In adult ICU COVID-19 patients on iMV, we observed no significant difference in outcomes between IL-6RA, steroids, or combination therapy. The lessons learned from this study may suggest a one-size-fits-all for the immunomodulatory effect of IL6-RA and steroid therapy is unlikely to be the answer. Their effect appears to be less pronounced if administered for subsets of patients with iMV. While observational studies cannot suggest causality, further RCTs are needed to enhance the clinicians’ personalized clinical decision-making for COVID-19-associated ARDS.
Data Availability
The datasets used and analyzed during the current report are available from the corresponding author on reasonable request.
Author contributions
MA: proposed the study, led the manuscript writing and did the data collection. AK: made the statistical analysis of this study. All authors made substantial contributions to the concept, design, data collection, and have further participated in the analysis and interpretation of the data. All authors participated in revising the manuscript for important intellectual content and approved the final version to be published. We confirm that the authorship followed the uniform requirements for manuscripts submitted to biomedical journals.
Availability of data and material
The datasets used and analyzed during the current report are available from the corresponding author on reasonable request.
Supplemental Content
Supplementary file 1. Additional details for statistical analysis
International System of Units (SI) were used to represent laboratory data unless otherwise specified. The data were checked for integrity and outliers in the laboratory data (numbers > the 99th percentile) were removed based on clinical judgment and whether the extreme values were clinically plausible.
Marginal structural modeling (MSM)
allows for proper adjustment of time-dependent confounding variables, immortal time bias, and indication bias. Immortal time bias refers to the requirement that patients survive long enough to receive the intervention of interest, leading to a potential incorrect overestimation of a positive treatment effect. Indication bias from time-varying confounding variables refers to having an association related to the indication of the intervention that evolves throughout the course of an illness.
MSM using propensity scores were used to compare the primary and secondary outcomes between the three regimens after adjusting for non-time and time-varying covariates [1]. All patients had baseline data on ICU admission. However, the follow-up time varied between patients. Moreover, some variables, such as therapeutic anticoagulation, varied daily. MSMs were used to adjust for variability between patients who received the three regimens. Thus, handling time-varying covariates when weighing observations may have the following limitations:
1. the longitudinal outcome and time-varying covariate may not be available at the same time points: 2. the longitudinal outcome at a particular time point, t, may depend not only on the value of the covariate at the same time point but also at those at other time points. Thus, weighted regressions may be used to account for time-varying confounders using MSMs.
For dichotomous outcomes, the corresponding mixed effects model, namely a mixed effects logistic regression gives fixed effects coefficients that have an interpretation conditional on the random effects. Most often, this is not the interpretation we want. The Generalized Estimating Equations (GEE) approach does give coefficients with a marginal / population-averaged interpretation.
However, an additional practical point that also need to consider is missing data and almost always we have to deal with incomplete data. With regard to this point, mixed models give valid results under the less stringent missing at random assumption compared to the (standard not weighted) GEE that give valid results under the less realistic missing completely at random assumption.
Taking both points (i.e., interpretation and missing data) into account, would most often like to fit a mixed model to be more protected for the missing data but want to obtain parameters that have a population averaged interpretation. An early solution towards this direction was the marginalized mixed models proposed by Heagerty (2), but, in general, these are computationally intensive to fit. A more recent approach that seems to solve the problem has been proposed by Hedeker et al (3). This is implemented in the function marginal_coefs() of the R package GLMMadaptive.
The fixed effects estimates in mixed models with nonlinear link functions have an interpretation conditional on the random effects. However, often we wish to obtain parameters with a marginal / population averaged interpretation, which leads many researchers to use generalized estimating equations, and dealing with potential issues with missing data. Nonetheless, recently Hedeker et al. have proposed a nice solution to this problem. Their approach is implemented in function marginal_coefs(). For example, for model fm we obtain the marginalized coefficients.
The function calculates the marginal log odds ratios in our case (because we have a binary outcome) using a Monte Carlo procedure with number of samples determined by the M argument. Standard errors for the marginalized coefficients are obtained by setting std_errors = TRUE in the call to marginal_coefs(), and require a double Monte Carlo procedure for which argument K comes also into play.
Generalized estimating equations (GEE) were constructed to evaluate the impact of the treatment (t) on the outcome (d) and to handle the time-varying covariate (v) in the panel dataset. Since the outcome was a continuous variable, a generalized linear Gaussian family with identity link was used. Auto-regressive (AR1) correlation structure was selected since this was a time series (panel) data; we expected the correlation to decay as the outcome values were farther away from the time of interest.
Covariates included in the model were age, sex, ethnicity, asthma/COPD, ARDS grade (none, mild, moderate or severe) at admission, history of CVD (diabetes, hypertension, coronary heart diseases or congestive heart failure), and the baseline lowest FiO2. Other covariates included the use of any anticoagulation, therapeutic anticoagulation, hydroxychloroquine, azithromycin, and antivirals, (including remdesivir); vasopressors; and neuromuscular blockers. The effect of adding steroids, IL-6 antagonists was assessed by including each variables as a dummy coded variable in the marginal model. In addition, an interaction term was added to the model to assess whether the combination of both drugs was associated with a statistically significant effect size.
inverse probability weighting (IPTW) are weights assigned to each observation across time conditioned on the previous exposure history, which are then multiplied to generate a single weight for a subject (4). Similar to conventional propensity score estimation, IPTW is generated using either a logit or probit model that regresses covariates to a treatment group (exposure) variable. With IPTW, the previous exposure history is incorporated to the propensity score estimation, which is time-varying. Standardized weights in a longitudinal setting are estimated. The numerator contains the probability of the observed exposure at each time point conditioned on the observed exposure history of the previous time point and the observed non-time varying covariates. The denominator contains the probability of the observed exposure at each time point conditioned on the observed exposure history of the previous time point, the observed time-varying covariates history at the current time point, and the non-time varying covariates. In standardized weights, the time-varying confounders are captured in the denominator but not in the numerator. However, the non-time varying (also known as fixed-time) covariates are captured in both the numerator. Due to the longitudinal nature of the data, the GLMMadaptive package was used for model fitting. The Package GLMMadaptive provides a suit of functions for fitting and post-processing mixed effects models for grouped/clustered outcomes which have a distribution other than a normal distribution. In particular, let yi denote a vector of grouped/clustered outcome for the i-th sample unit (i=1,…,n).
Standardized mean difference (SMD) were reported for comparison between cohorts before and after matching
Handling of missing data
The generalized boosted model (GMB) creates a separate level during matching called NA and creates weights based on the distribution of missing values. So basically it considers NA as a level of its own. The weights are created for all respondents and are taking the missing values into consideration. Imputation requires that data is missing at random which is not the case here. Moreover, the variables we used had few missing values so there is no problem which would not have been the case if we used the SOFA and APACHE which had a lot of missing data.
The pseudo population
The populations are weighed by the probability of receiving treatments based on other confounders. The weights were produced using generalized boosted models and these weights were used when doing further analysis. The pseudo population is the sum of weights for each group rather than considering each observation as 1. Thus, the numbers can increase or decrease compared to the initial number of the group depends on the assigned weights. The aim to emulate randomized trial to create a pseudo study cohort, where the weighted version can balance off the covariate bias and mimic a randomized treatment assignment situation
Supplementary File 2: Collaborative Co-Author List
¥ On behalf of the Society of Critical Care Medicine Discovery Viral Infection and Respiratory Illness Universal Study (VIRUS): COVID-19 Registry Investigator Group
Belgium
Centre Hospitalier Jolimont: Jean-Baptiste Mesland, Pierre Henin, Hélène Petre, Isabelle Buelens, Anne-Catherine Gerard
The Brugmann University hospital, Bruxelles: Philippe Clevenbergh
Bosnia and Herzegovina
University Clinical Hospital, Mostar, Bosnia and Herzegovina: Dragana Markotić, Ivana Bošnjak
Columbia
Clinica Medical SAS: Oscar Y Gavidia, Felipe Pachon, Yeimy A Sanchez
Croatia
Clinical Hospital Center Rijeka, Rijeka, Croatia: Danijel knežević
Egypt
Helwan University: Mohamed El Kassas, Mohamed Badr, Ahmed Tawheed, Ahmed Tawheed, Hend Yahia
Hondurus
Honduras Medical Center: Jose Luis Ramos Coello, Guillermo Perez, Ana Karen Vallecillo Lizardo, Gabina María Reyes Guillen, Helin Archaga Soto
Hungary
Uzsoki Teaching Hospital: Csaba Kopitkó, Ágnes Bencze, István Méhész, Zsófia Gerendai,
India
Jawaharlal Institute of Postgraduate Medical Education and Research, Pondicherry: Anusha Cherian, Sreejith Parameswaran, Magesh Parthiban, Menu Priya A.
Maulana Azad Medical College and Lok Nayak hospital: Mradul Kumar Daga, Munisha Agarwal, Ishan Rohtagi
Japan
Center Hospital of the National Center for Global Health and Medicine: Wataru Matsuda, Reina Suzuki
Sapporo City General Hospital: Yuki Itagaki, Akira Kodate, Reina Suzuki, Akira Kodate, Yuki Takahashi, Koyo Moriki
Tokyo Medical and Dental University: Hidenobu Shigemitsu, Yuka Mishima, Nobuyuki Nosaka, Michio Nagashima
Hiroshima University: Michihito Kyo
Mexico
Hospital Universitario, Universidad Autonoma de Nuevo León: Rene Rodriguez-Gutierrez, Jose Gerardo Gonzalez-Gonzalez, Alejandro Salcido-Montenegro, Adrian Camacho-Ortiz
Pakistan
Dow University Hospital: Muhammad Sohaib Asghar, Mashaal Syed, Syed Anosh Ali Naqvi
Russia
Kuban State Medical University with affiliation Territorial Hospital #2: Igor Borisovich Zabolotskikh, Konstantin Dmitrievich Zybin, Sergey Vasilevich Sinkov, Tatiana Sergeevna Musaeva
Saudi Arabia
King Saud University: Mohammed A Almazyad, Mohammed I Alarifi, Jara M Macarambon, Ahmad Abdullah Bukhari, Hussain A. Albahrani, Kazi N Asfina, Kaltham M Aldossary
King Saud bin Abdulaziz University for Health Sciences and King Abdullah International Medical Research Center: Yaseen M Arabi, Sheryl Ann Abdukahil
King Fahad Armed Forces Hospital: Razan K Alamoudi, Hassan M. AlSharif, Sarah A. Almazwaghi, Mohammed S Elsakran, Mohamed A Aid, Mouaz A Darwich, Omnia M Hagag, Salah A Ali, Alona rocacorba, Kathrine Supeña, Efren Ray Juane, Jenalyn Medina, Jowany Baduria
King Faisal Specialist Hospital & Research Centre-Riyadh: Marwa Ridha Amer, Mohammed Abdullah Bawazeer, Khalid Maghrabi, Abid Shahzad Butt, Talal I. Dahhan, Eiad Kseibi, Syed Moazzum Khurshid, Muath Rabee, Mohammed Abujazar, Razan Alghunaim, Maal Abualkhair, Abeer Turki AlFirm, Ali Al-Janoubi
Serbia
UMC Zvezdara, Belgrad: Bojan Kovacevic, Jovana Bojicic
Institute for Pulmonary Diseases of Vojvodina, Sremska Kamenica: Ana Andrijevic, Srdjan Gavrilovic, Vladimir Carapic
Spain
Hospital Universitario La Paz: Santiago Y. Teruel, Belen C. Martin, Santiago Y. Teruel
United States
Boston University School of Medicine, Boston, MA: Allan J. Walkey, Sushrut S. Waikar, Michael A. Garcia, Mia Colona, Zoe Kibbelaar, Michael Leong, Daniel Wallman, Kanupriya Soni, Jennifer Maccarone, Joshua Gilman, Ycar Devis, Joseph Chung, Munizay Paracha, David N. Lumelsky, Madeline DiLorenzo, Najla Abdurrahman, Shelsey Johnson
Albany Medical Center: Suzanne Barry, Christopher Woll, Gregory Wu, Erin Carrole, Kathryn Burke, Mustafa Mohammed
AnMed Health: Abhijit A Raval, Andrea Franks
Ascension St.Vincent Hospital: Anmol Kharbanda, Sunil Jhajhria, Zachary Fyffe
Ascension/St. Thomas Research Institute West Campus: Stephen Capizzi, Bethany Alicie, Martha
Green, Lori Crockarell, Amelia Drennan, Kathleen Dubuque, Tonya Fambrough, Nikole Gasaway, Briana Krantz, Peiman Nebi, Jan Orga, Margaret Serfass, Alina Simion, Kimberly Warren, Cassie Wheeler, CJ Woolman
Augusta Health: Andrew S. Moyer, George M. Verghese
Augusta University Medical Center: Andrea Sikora Newsome, Christy C. Forehand, Rebecca Bruning, Timothy W. Jones
Banner University Medical Center-Tucson: Jarrod M Mosier, Karen Lutrick, Beth Salvagio Campbell, Cathleen Wilson, Patrick Rivers, Jonathan Brinks, Mokenge Ndiva Mongoh, Boris Gilson
Baylor College of Medicine, Baylor St. Lukes Medical Center: Christopher M Howard, Cameron McBride, Jocelyn Abraham, Orlando Garner, Katherine Richards, Keegan Collins, Preethi Antony, Sindhu Mathew
Baylor Scott & White Health: Valerie C. Danesh, Gueorgui Dubrocq, Amber L. Davis, Marissa J Hammers, ill M. McGahey, Amanda C. Farris, Elisa Priest, Robyn Korsmo, Lorie Fares, Kathy Skiles, Susan M. Shor, Kenya Burns, Corrie A Dowell, Gabriela “Hope” Gonzales, Melody Flores, Lindsay Newman, Debora A Wilk, Jason Ettlinger, Jaccallene Bomar, Himani Darji, Alejandro Arroliga, Alejandro C Arroliga, Corrie A. Dowell, Gabriela Hope Conzales, Melody Flores, Lindsay Newman, Debora A. Wilk, Jason Ettlinger, Himani Darji, Jaccallene Bomar
Beth Israel Deaconess Medical Center: Valerie M. Banner-Goodspeed, Somnath Bose, Lauren E. Kelly, Melisa Joseph, Marie McGourty, Krystal Capers, Benjamin Hoenig, Maria C. Karamourtopoulos, Anica C. Law, Elias N. Baedorf Kassis
Cedars Sinai Medical Center: Pooja A. Nawathe, Isabel Pedraza, Jennifer Tsing, Karen Carr, Anila Chaudhary, Kathleen Guglielmino
Chambersburg Hospital: Raghavendra Tirupathi, Alymer Tang, Arshad Safi, Cindy Green, Jackie Newell
Children’s Hospital Colorado, University of Colorado Anschutz Medical Campus: Katja M. Gist, Imran A Sayed, John Brinton, Larisa Strom
Children’s Hospital of Philadelphia: Kathleen Chiotos, Allison M. Blatz, Giyoung Lee, Ryan H. Burnett, Guy I. Sydney, Danielle M. Traynor
Cox Medical Center Springfield: Steven K. Daugherty, Sam Atkinson, Kelly Shrimpton
Detar Family Medicine residency: Sidney Ontai, Brian Contreras, MD, Uzoma Obinwanko, Nneka Amamasi, Amir Sharafi
Detroit Medical Centre: Sarah Lee, Zahia Esber, Chetna Jinjvadia
George Washington University: David P. Yamane, Ivy Benjenk, Nivedita Prasanna
Howard University Hospital: Orma Smalls
Jacobs Medical Center UC San Diego Health – La Jolla: Atul Malhotra, Abdurrahman Husain, Qais Zawaydeh
JPS Health Network: Steven Q. Davis, Valentina Jovic, Valentina Jovic, Max Masuda, Amanda Hayes
Lakes Region General Hospital: Michael Smith, William Snow, Riley Liptak, Hannah Durant, Valerie Pendleton, Alay Nanavati, Risa Mrozowsk
LifeBridge Health/Sinai and Northwest Hospitals: Namrata Nag, Jeff Brauer, Ashwin Dharmadhikari, Sahib Singh, Franco Laghi, Ghania Naeem, Andrew Wang, Kevin Bliden, Amit Rout, Jaime Barnes, Martin Gesheff, Asha Thomas, Melbin Thomas, Alicia R. Liendo, Jovan Milosavljevic, Kenan Abbasi, Nicholas B. Burley, Nicole Rapista, Samuel Amankwah, Sanjay K Poudel, Saroj Timilsina, Sauradeep Sarkar, Oluwasayo Akinyosoye, Shashi K. Yalamanchili, Sheena Moorthy, Sonia Sugumar, Jonathan Ford, Martin C. Taylor, Charlotte Dunderdale, Alyssa Henshaw, Mary K. Brunk, Jessica Hagy, Shehryar Masood, Sushrutha Sridhar
Loyola University Medical Center: Yuk Ming Liu, Sarah Zavala, Sarah Zavala, Esther Shim
M Health-Fairview, University of Minnesota: Ronald A. Reilkoff, Julia A. Heneghan, Sarah Eichen, Lexie Goertzen, Scott Rajala, Ghislaine Feussom, Ben Tang
MacNeal Hospital Loyola Medicine: Christine C. Junia, Robert Lichtenberg, Hasrat Sidhu, Diana Espinoza, Shelden Rodrigues, Maria Jose Zabala, Daniela Goyes, Ammu Susheela, Buddhi Hatharaliyadda, Naveen Rameshkumar, Amulya Kasireddy, Genessis Maldonado, Lisseth Beltran, Akshata Chaugule, Hassan Khan
Mayo Clinic Arizona: Rodrigo Cartin-Ceba, Ayan Sen, Amanda Palacios, Giyth M. Mahdi
Mayo Clinic Rochester: Rahul Kashyap, Ognjen Gajic, Vikas Bansal, Aysun Tekin, Amos Lal, John C. O’Horo, Neha N. Deo, Mayank Sharma, Shahraz Qamar, Cory J. Kudrna, Juan Pablo Domecq
Mayo Clinic, Eau Claire: Abigail T. La Nou, Marija Bogojevic
Mayo Clinic, Florida: Devang Sanghavi, Pramod Guru, Karthik Gnanapandithan, Hollie Saunders, Zachary Fleissner, Juan Garcia, Alejandra Yu Lee Mateus, Siva Naga Yarrarapu
Mayo Clinic, Mankato: Syed Anjum Khan, Nitesh Kumar Jain, Thoyaja Koritala
Medical Center Health System, Odessa: Alexander Bastidas, Gabriela Orellana, Adriana Briceno Bierwirth, Eliana Milazzo, Juan Guillermo Sierra, Thao Dang
Medical Center Navicent Health: Amy B. Christie, Dennis W. Ashley, Rajani Adiga
Mercy Hospital and Medical Center, Chicago: Travis Yamanaka, Nicholas A. Barreras, Michael Markos, Anita Fareeduddin, Rohan Mehta
Mercy Hospital, Saint Louis: Chakradhar Venkata, Miriam Engemann, Annamarie Mantese
Millard Fillmore Suburban Hospital: Anna Eschler, Mary Hejna, Emily Lewandowski, Kristen Kusmierski, Clare Martin
Montefiore Medical Center The Bronx: Jen-Ting Chen, Aluko Hope, Zoe Tsagaris, Elise Ruen, Aram Hambardzumyan
OSF Saint Francis Medical Center: Bhagat S. Aulakh, Sandeep Tripathi, Jennifer A. Bandy, Lisa M. Kreps, Dawn R. Bollinger, Jennifer A. Bandy
OSF Saint Francis Medical Center: Bhagat S. Aulakh, Sandeep Tripathi, Jennifer A. Bandy, Lisa M. Kreps, Dawn R. Bollinger, Jennifer A. Bandy
Parkview Health System, Fort Wayne: Roger Scott Stienecker, Andre G. Melendez, Tressa A. Brunner, Sue M Budzon, Jessica L. Heffernan, Janelle M. Souder, Tracy L. Miller, Andrea G. Maisonneuve
Samaritan Health Services: Brian L. Delmonaco, Anthony Franklin, Mitchell Heath
Sarasota Memorial Hospital: Antonia L. Vilella, Sara B. Kutner, Kacie Clark, Danielle Moore
St. Joseph’s Candler Health System: Howard A. Zaren, Stephanie J. Smith, Grant C. Lewis, Lauren Seames, Cheryl Farlow, Judy Miller, Gloria Broadstreet
St.Mary Medical Center, Langhorne: Umang Patel, Jordesha Hodge, KrunalKumar Patel, Shivani Dalal, Himanshu Kavani, Sam Joseph
Stamford Health: Michael A. Bernstein, Ian K. Goff, Matthew Naftilan, Amal Mathew, Deborah Williams, Sue Murdock, RN, Maryanne Ducey, Kerianne Nelson
Stanford Hospital and Clinics: Paul K Mohabir, Connor G O’Brien, Komal Dasani
The Children’s Hospital at OU Medicine: Neha Gupta, Tracy L Jones, Shonda C Ayers, Amy B Harrell, Dr. Brent R Brown
The University of Tennessee Medical Center: Megan Edwards, Caleb Darby, Kristy Page, Amanda Brown, Jessie McAbee
Thomas Jefferson University Hospital: Katherine A. Belden, Michael Baram, Devin M. Weber, Rosalie DePaola, Yuwei Xia, Hudson Carter, Aaron Tolley, Mary Barletta
Truman Medical Centers: Mark Steele, Laurie Kemble
Tulane University Medical Center and University Medical Center New Orleans: Joshua L. Denson, A. Scott Gillet, Margo Brown, Rachael Stevens, Andrew Wetherbie, Kevin Tea, Mathew Moore
UNC Medical Center: Benjamin J Sines, Thomas J Bice
University Medical Center of Southern Nevada Las Vegas; University of Nevada, Las Vegas: Rajany V. Dy, Alfredo Iardino, Jill Sharma, Julia Christopher, Marwan Mashina, Kushal Patel
University of Alabama at Birmingham: Erica C. Bjornstad, Nancy M. Tofil, Scott House, Isabella Aldana
University of Arkansas for Medical Sciences: Nikhil K. Meena, Jose D. Caceres, Nikhil K Meena, Sarenthia M. Epps, Harmeen Goraya, Kelsey R. Besett, MD, Ryan James, Lana Y. Abusalem, Akash K. Patel, Lana S Hasan
University of Cincinnati: Dina Gomaa B.S., Michael Goodman, Devin Wakefield, Anthony Spuzzillo, John O. Shinn II
University of Iowa Carver College of Medicine: Patrick W. McGonagill, Colette Galet, Janice Hubbard, David Wang, Lauren Allan, Aditya Badheka, Madhuradhar Chegondi
University of Kansas Medical Center: Usman Nazir, Garrett Rampon, Jake Riggle, Nathan Dismang
University of Louisville Hospital: Ozan Akca, Rainer Lenhardt, Rodrigo S. Cavallazzi, Ann Jerde, Alexa Black, Allison Polidori, Haily Griffey, Justin Winkler, Thomas Brenzel
University of Michigan Health System: Pauline Park, Andrew Admon, Sinan Hanna, Rishi Chanderraj, Maria Pliakas, Ann Wolski, Jennifer Cirino
University of Missouri, Columbia: Dima Dandachi, Hariharan Regunath, Maraya N. Camazine, Grant. E. Geiger, Abdoulie O. Njai, Baraa M. Saad
University of Vermont Larner College of Medicine: Renee D. Stapleton, Anne E. Dixon, Olivia Johnson, Sara S. Ardren, Stephanie Burns, Anna Raymond, Erika Gonyaw, Kevin Hodgdon, Chloe Housenger, Benjamin Lin, Karen McQuesten, Heidi Pecott-Grimm, Julie Sweet, Sebastian Ventrone
Wake Forest University School of Medicine; Wake Forest Baptist Health Network: Ashish K. Khanna, Lynne Harris, Bruce Cusson, Jacob Fowler, David Vaneenenaam, Glen McKinney, Imoh Udoh, Kathleen Johnson
Yale New Haven Health New Haven: Kevin Sheth, Abdalla Ammar, Mahmoud Ammar, Victor Torres Lopez, Charles Dela Cruz, Akhil Khosla, Samir Gautam
Acknowledgments
We thank the SCCM Discovery VIRUS data registry and the Collaborative Co-authors listed in Supplementary file 2 for providing and maintaining the database of hospitalized COVID-19 patients which made this study possible. We are particularly grateful to Professor Yaseen M Arabi at King Saud bin Abdulaziz University for his inputs and review of this manuscript. We thank the ICU physicians, nurses, respiratory therapists, physiotherapists, pharmacists across the globe for their efforts to save the lives of ICU patients and their braveness in fighting the COVID-19 pandemic.
Footnotes
¥ On behalf of the Society of Critical Care Medicine Discovery Viral Infection and Respiratory Illness Universal Study (VIRUS): COVID-19 Registry Investigator Group
Financial support and sponsorship: None
The VIRUS Registry is supported, in part, by the Gordon and Betty Moore Foundation, and Janssen Research & Development, LLC
We added revised section for the statistical analysis, results and discussion section based on reviewers comments