Skip to main content

Exploration of different statistical approaches in the comparison of dopamine and norepinephrine in the treatment of shock: SOAP II

Abstract

Background

Exploring clinical trial data using alternative methods may enhance original study’s findings and provide new insights. The SOAP II trial has been published more than 10 years ago; but there is still some speculation that some patients may benefit from dopamine administration for shock management. We aimed to reanalyse the trial under different approaches and evaluate for heterogeneity in treatment effect (HTE).

Methods

All patients enrolled in SOAP II were eligible for reanalysis. We used a variety of methods including the win-ratio (WR), a Bayesian reanalysis stratified according to shock type, and both a risk-based and effect-based explorations for HTE. The methods were applied to different endpoints, including a hierarchy of death, new use of renal-replacement therapy (RRT), and new-onset arrhythmia; 28-day mortality; a composite endpoint (mortality, new use of RRT, and new-onset arrhythmia), and days alive and free of ICU at 28-days (DAFICU28).

Results

A total of 1679 patients were included (average age was 64.9 years, 57% male, 62% with septic and 17% with cardiogenic shock). All analysis favoured norepinephrine over dopamine. Under the WR approach, dopamine had fewer wins compared to norepinephrine (WR 0.79; 95% confidence intervals [CI] 0.68–0.92; p = 0.003), evident in both cardiogenic and septic shock subgroups. The Bayesian reanalysis for type of shock showed, for dopamine, a probability of harm of 0.95 for mortality, > 0.99 probability of harm for composite endpoint, and 0.91 probability of harm for DAFICU28. The fewer DAFICU28 with dopamine was more apparent in those with cardiogenic shock (0.92). Under the risk-based HTE, there was a high probability that dopamine resulted fewer DAFICU28 in the highest quartile of predicted mortality risk. The effect-based HTE assessment model did not recommended dopamine over norepinephrine for any combination of possible modifiers including age, type of shock, presence of cardiomyopathy, and SOFA score. Receiving dopamine when the effect-based model recommended norepinephrine was associated with an absolute increase in composite endpoint of 6%.

Conclusion

The harm associated with the use of dopamine for the management of shock appears to be present in both septic and cardiogenic shock patients. There was no suggestion of any subgroup in which dopamine was found to be favourable over norepinephrine.

Background

Choice of vasopressor is a crucial issue in patients with shock admitted to intensive care units (ICU). Dopamine was commonly used as first line vasopressor in patients with shock; however, norepinephrine has been favoured over the last couple decades after trials suggesting that norepinephrine could be associated with improved outcomes and fewer adverse effects [1, 2]. The landmark trial comparing norepinephrine and dopamine is the Comparison of Dopamine and Norepinephrine in the Treatment of Shock trial (SOAP II) which enrolled 1679 critically ill patients with shock to either norepinephrine or dopamine. SOAP II reported neutral results for the primary outcome of 28-day mortality (odds ratio [OR] with dopamine, 1.17; 95% confidence interval [CI] 0.97–1.42; p = 0.10) [3]. However, the trial found important differences, including more adverse events with dopamine [3].

Interpreting neutral clinical trials that include a heterogeneous population is challenging. On one hand, the effect size of SOAP II points towards possible harm with the use of dopamine; and a secondary subgroup analysis suggested that harm could be driven mostly by patients with cardiogenic shock, but the p-value for interaction was 0.12, making a strong statement not completely compatible under the frequentist framework. In addition, physicians continue to use dopamine as a vasopressor agent, suggesting that it may still be considered suitable in some conditions. There is growing interest for different statistical methods that can either provide hierarchical assessment of composite endpoints (such as the win ratio [WR]) [4], probabilistic statements on the results (mostly through Bayesian methods [5]) and exploring heterogeneity in treatment effects (HTE), which may provide additional information from randomized trials [6].

We therefore reanalyzed the SOAP II trial using different statistical methods, including a hierarchical assessment of composite endpoints using WR [4], Bayesian reanalysis accounting for interaction of shock type and intervention [5], probabilistic methods of HTE based on risk [6] and effect-based approaches [7]. We hypothesized that these methods would provide additional insights that would both serve to showcase novel methods of analysis for randomized trials and augment the findings from the original SOAP II results.

Methods

Patients

All patients included in the SOAP II trial (raw original dataset).

SOAP II overview

SOAP II was a multicenter, randomized, clinical trial, where patients were randomly allocated to receive either dopamine or norepinephrine as the first-line vasopressor therapy to manage shock. If the blood pressure could not be maintained with a dose of 20 μg per kilogram of body weight per minute for dopamine or a dose of 0.19 μg per kilogram per minute for norepinephrine, open-label norepinephrine, epinephrine, or vasopressin could be added. The primary endpoint was all-cause 28-day mortality. Fully anonymized data provided by SOAP II investigators was used for this analysis.

Missing data handling

Continuous or factor variables with missing values up to 1% of the sample size were imputed using median or most common level.

Endpoints and rationale

For the WR approach, two hierarchical approaches were used: (1) 28-day mortality, use of renal replacement therapy (RRT), and occurrence of new arrhythmia was used, and (2) 28-day mortality and ICU length-of-stay (LOS). Those endpoints encompass relevant hierarchies that could have been chosen in SOAP II; both use mortality as first endpoint (to avoid competing risks). The difference is how ties in mortality are dealt with. The first hierarchy adds two additional layers including use of RRT and arrhythmias, which could be affected by dopamine use. The second hierarchy, a simpler combination of mortality and LOS, is conceptually similar to days alive and free of ICU endpoints. There was no stopping due to early ties. WR above 1 would suggest dopamine is better (i.e., results in more wins). For the Bayesian reanalysis according to type of shock and HTE explorations, three different endpoints were used: 28-day mortality; a composite of 28-day mortality, new use of RRT and occurrence of new arrhythmia (hence using the first WR hierarchy in a single composite endpoint); and days alive and free of ICU at 28 days (DAFICU28; which conceptually approaches the second win ratio hierarchy used). Patients that died within 28 days received zero DAFICU28.

Statistical analysis

The analysis method depended on the approach used to explore SOAP II:

WR analysis

We calculated the WR for dopamine over norepinephrine using both hierarchical endpoints stratified according to shock type (grouped in cardiogenic, septic, and other; the “other” category used for analysis includes all types of shock that are not cardiogenic and septic; detailed explanation is provided in Table 1 legend) [8]. Results are presented as a WR with 95% confidence interval (CI), and p values. Values of WR below 1 suggest a harm for dopamine in the composite endpoint (i.e., fewer wins).

Table 1 Patient features according to randomization arm

Bayesian reanalysis of the main trial accounting for the interaction of the intervention with shock type

The effect of dopamine versus norepinephrine for 28-days mortality was assessed using a Bayesian regression model for randomization arm on outcome adjusted for type of shock and the interaction between type of shock and intervention. A neutral prior was used for the effect of the intervention (neutral prior for log odds ratio with mean zero and standard deviation of 0.355 [5]). Relevance of the interaction was estimated according for the Bayes factor (BF) for a model with and without interaction, with a BF above 10 being supportive of presence of interaction (i.e., the model with interaction is at least 10 times more likely to fit the data than the model without interaction). If the BF was above 10, we planned to explore the effects of the intervention on each shock type; otherwise, we planned to report only the result of the model without interaction, including the odds ratio and absolute risk differences for mortality with dopamine when compared to norepinephrine together with its probability of direction and its 95% high density intervals (HDI). A region of practical equivalence (ROPE) was defined as differences in absolute risks of 0.01. A probability of direction above 0.90 was suggestive of an effect; a probability of direction above 0.95 with less than 0.05 of the posterior mass of absolute risk difference within the ROPE (%ROPE < 0.05) was considered as strong evidence of effect. A similar approach was applied for the composite endpoint of death, new use of RRT or new arrhythmia. For DAFICU28, a Bayesian linear model was used with; cut off for probability of direction as well as BF criteria were similar; ROPE was defined as a difference of plus or minus one day for DAFICU28.

Risk-based heterogeneity in treatment effect

This approach followed the PATH statement [6]. We first performed a first level customization of APACHE II in the trial population and assessed its discrimination capability (through receiver-operating characteristics curve) and calibration (using calibration belts) [9]. Customization was necessary due to the lack of crude data on reason for admission allowing original APACHE II predictions to be calculated. We then generated quartiles of predicted risk of death based on customized APACHE II [10] (cAPACHE II) predicted probabilities. The prediction risks quartiles were then used for a Bayesian logistic regression with 28-day mortality as outcome together with study randomization arm and their interaction. Presence of risk-based HTE was defined as a BF above 10 for the model with versus a model without the interaction for cAPACHE II quartiles for the primary outcome. These risk categories from APACHE II were also used for risk-based HTE explorations for composite endpoint and DAFICU28. Results are presented as odds ratios, absolute risk differences, or mean differences, according to enrollment arm with their 95% HDI. Same rules for stating association and ROPE were used.

Effect-based heterogeneity in treatment effects

This analysis was designed to apply the S-learner technique using SOAP II for the composite endpoint [7]. In brief, the S-learner is a technique where a prediction model allowing interaction between intervention arm and key covariates is trained in part of the original dataset (“train set”; defined as a 0.60 of total SOAP II sample size) [7]. The model was adjusted for age, type of shock, presence of cardiomyopathy at baseline, and sequential organ failure assessment (SOFA [11]) score, all interacting with intervention arm. The model was then used to generate two counterfactual predictions in remaining sample (“test” set); counterfactuals were generating by first changing all patients to dopamine arm, obtaining predictions, then changing all patients to norepinephrine arm and obtaining a second set of predictions. The difference between predictions in the test set are an estimate of individualized treatment effects (ITE) for those patients. The model would then recommend one specific arm if the probability of benefit for that arm in that specific patient was above 0.90; if the differences in expected predicted outcomes did not reach the threshold, the model would not make a recommendation. If a recommendation could be made for at least 25% of the test set, the association between model recommendations and outcome would be assessed using a Bayesian logistic regression model with interaction between recommendation and study arm. We reported outcome differences according to recommended versus received treatment in an exploratory fashion following the same principles as the other Bayesian models.

Additional exploratory analyses

Additional exploratory analysis included assessment on whether heart rate would be an effect modifier for dopamine use; this included direct interaction analyses and an analysis also adjusted by type of shock.

All analysis were performed in R version 4.3.2 with packages {BuyseTest} [8], {brms} [12], and {marginaleffects} [13].

Results

A total of 1679 patients were included. Patient features are shown in Table 1. Average age was 64.9 years (standard deviation—SD—14.6) and most patients were male (57%). Septic shock was the most common etiology for shock (1044 patients—62%). Mechanical ventilation was used at enrollment in 1195 patients (71.2%). A summary of methods and results is shown in Table 2.

Table 2 Summary of analyses, results, and interpretations discussed in the manuscript

WR results

A cross tabulation of the components of the first WR hierarchy is shown in eTable 1. The overall WR for dopamine versus norepinephrine for the first hierarchical approach was 0.79 (95% CI 0.68–0.92; p = 0.003), suggesting that dopamine resulted in less wins than norepinephrine. For the secondary hierarchy, WR was 0.96 (95% CI 0.86–1.08; p = 0.51). Results for each hierarchical level according to strata for both approaches are shown in Fig. 1 and eFigure 1.

Fig. 1
figure 1

Results for the win ration (WR) analysis for the first hierarchical approach (A and B). A Stratified win ratio results for the first hierarchical approach (mortality, use of renal replacement therapy—RRT, and occurrence of arrhythmia). B Stratified win ratio results for the second hierarchical approach (mortality and ICU length-of-stay—LOS)

Bayesian reanalysis results for the interaction with shock type

For the standard Bayesian reanalysis, the BF for the primary endpoint for a model with interaction versus without interaction was 0.8, suggesting that it is unlikely that a significant interaction between enrolment arm and type of shock for this endpoint. There was a suggestion that dopamine was worse than norepinephrine (probability of harm of 0.95), with 0.09 of the posterior within the region of practical equivalence. The median increase in mortality with dopamine was 0.04; 95% HDI from − 0.01 to 0.09 (odds ratio of 1.18; 95% HDI from 0.96 to 1.39). Model syntax and diagnostics are shown in ESM (eFigure 2).

For the composite endpoint result, BF for model with interaction was 2.84, not supportive of relevance of the interaction between study arm and type of shock. There was strong association between dopamine and more composite endpoints than norepinephrine in this analysis (probability of harm > 0.99; %ROPE < 0.01; absolute risk difference of 0.08, 95% HDI 0.04–0.12; odds ratio of 1.44, 95% HDI from 1.18 to 1.73). Model syntax and diagnostics are shown in ESM (eFigure 3).

For the DAFICU28 analysis, results were supportive of an interaction between enrollment arm and outcome (BF of 20). For all types of shock, there was an association between dopamine use and less DAFICU28 (mean difference of − 0.70 days; 95% HDI − 1.58 to 0.16 days; probability of harm of 0.91); these results varied according to shock type: probability of harm of dopamine was 0.50, 0.92, and 0.85 for patients with other types of shock, cardiogenic shock, and septic shock, respectively. Differences in DAFICU28 also varied according to shock type (0.00, 95% HDI from − 0.68 to 0.68; − 1.79, 95% HDI − 4.31 to 0.62; and − 0.65, 95% HDI − 1.92 to 0.56 for other types of shock, cardiogenic shock, and septic shock, respectively). Results for these analyses are shown in Fig. 2 (28-day mortality and composite endpoint) and Fig. 3 (DAFICU28). Model syntax and diagnostic plots are shown in the ESM (eFigure 4).

Fig. 2
figure 2

Bayesian reanalysis for 28-day mortality and composite endpoint. Blue areas represent benefit for dopamine, red areas represent harm for dopamine (benefit for norepinephrine), and dark red areas represent the region of practical equivalence. A Posterior distribution of adjusted risk difference for mortality. B Posterior distribution of adjusted risk difference for composite endpoint of mortality, use of renal replacement therapy, and mortality

Fig. 3
figure 3

Bayesian reanalysis for day-alive and free of ICU at 28 days (DAFICU28). A Overall difference in DAFICU28 in a model adjusted for shock type. B, C, and D show the posterior distribution according to shock type. Areas in blue mark benefit for dopamine; areas in red mark harm for dopamine; region of practical equivalence—ROPE—is highlighted in dark red

Risk-based THE

Discrimination and calibration of customized APACHE II is shown in eFigure 5. Model syntaxes and diagnostic plots are provided in the ESM (eFigure 6–8). There was no clear sign of presence of risk-based HTE based on quartiles of APACHE II for 28-day mortality or the composite endpoint (BF of 2.5 and 1.40, respectively). Since there was no strong evidence of interaction, no summaries according to quartiles were calculated for those endpoints. There was a strong suggestion for the interaction between DAFICU28 and quartiles of cAPACHE II predicted mortality (BF of 75), with a suggestion of harm for dopamine concentrated in the higher quartile of predicted mortality for cAPACHE II (mean difference of − 1.64 days; 95% HDI − 3.62 to 0.48; probability of harm of 0.95). Results are visually shown in Fig. 4 with numeric values in Table 3.

Fig. 4
figure 4

Risk-based HTE results for days alive and free of ICU at 28 days—DAFICU28—according to customized APACHE II predicted mortality and randomization arm. Note that the differences in DAFICU28 are more pronounced for more severe patients (last quartile). Numeric values are shown in Table 3. Parenthesis are used to denote open intervals and brackets are used to denote the interval is closed at that end

Table 3 Results for risk-based HTE analysis

Effect-based THE

Details of model syntax and technical aspects are shown in ESM (including eFigure 9 for model diagnostics). The results for individual prediction of ITE (differences in composite endpoint based on S-Learner) is shown in Fig. 5A. The S-learner could make a recommendation in 73% of the test set (489 out of 669 patients). The model recommended norepinephrine for 489 patients, could not make a recommendation in 180, and never recommended dopamine. A comparison between patients recommended norepinephrine versus those without recommendation is shown in eTable 2; patients recommended norepinephrine more frequently had septic shock (78%) and had higher overall SOFA score. Receiving dopamine when the model recommended norepinephrine was associated with a 0.06 increase in composite endpoint (95% HDI − 0.01 to 0.14, 0.95 probability of harm, %ROPE 0.06), visual representation in Fig. 5B. Receiving dopamine when no recommendation could be made was not clearly associated with changes in composite endpoint (median difference of 0.95% HDI from − 0.14 to 0.15).

Fig. 5
figure 5

A “Individualized” (conditional) treatment effects (ITE) for patients in the test set according to the model adjusted for adjusted for age, type of shock, presence of cardiomyopathy at baseline, and SOFA score, all interacting with intervention arm. The differences between counterfactuals for each patient are displayed in (A), with patient identification sorted according to increasing ITE and y-axis showing the median difference and 95% HDI for the counterfactuals; a high value represents an increase in the counterfactuals between dopamine and norepinephrine, hence the higher the conceivable benefit for norepinephrine versus dopamine the higher the ITE. Note that the model could not make a recommendation in 27% of all patients. B Differences in outcome according to enrollment arm and suggestion obtained from the S-Learner. Receiving norepinephrine when the model recommended norepinephrine was associated with a decrease in composite endpoint; for patients for whom no recommendation was made, no clear differences between norepinephrine and dopamine were found for patients for whom no recommendation was made

Additional analyses

Results were neutral for both the continuous interaction of heart rate and intervention arm (eFigure 10) and heart rate (divided in quartiles) study arm and type of shock (eFigure 11).

Discussion

Overall results

In this exploratory post-hoc analysis of the SOAP II trial using additional analytic strategies, we have expanded interpretation of the main trial findings. SOAP II reported neutral results for mortality but suggested that dopamine was associated with more adverse events, including arrhythmias [1]. The present analysis brings some additional important conclusions, some which are complementary to the original report and some that are confirmatory. From a confirmatory perspective, our results suggest that if mortality and arrhythmia are bundled in a composite endpoint, the harm associated with dopamine is clearly presented (as shown in the WR analysis using the first hierarchy and the Bayesian reanalysis using composite endpoint). From a complementary perspective, our results highlight that a simple Bayesian reanalysis adjusted for type of shock already suggested a high probability of harm for dopamine for 28-day mortality (0.95). If more granular endpoints were used, such as days alive and free of ICU at day 28, then an interaction between shock type and outcome became more apparent, with a suggestion of increased harm with use of dopamine in cardiogenic shock. The heterogeneity in treatment effect for dopamine versus norepinephrine was also evident in risk-based HTE assessment based on APACHE II scores, making it uncertain whether type of shock or simply overall baseline risk are the mediators of the effect. A simple effect-based model based on S-learner confirmed the results with the model suggesting use of norepinephrine only for over 70% of the test set and never recommending dopamine.

Why is this analysis necessary now?

Revisiting trials through secondary exploratory analysis is a well-known phenomenon in critical care [14]. Traditionally, secondary analysis of a trial is done using specific subgroups of interest, which are prone to low power and magnification errors for positive results [15]. More recently, trials have started being reanalyzed using different endpoints and more advanced analytic methods, including Bayesian statistics or hierarchical endpoints [16, 17], often providing added insights beyond the original trial [16], but also generating concerns that excessive reanalysis and explorations will lead to spurious results or could be used to flip the primary conclusions of a clinical trial, specifically regarding the case of Bayesian analyses [18]. Dopamine is no longer recommended as first-line vasopressor in adults with septic shock [19] mostly due to SOAP II and the Patel et al. trials [1, 3, 20]. However, several guidelines for management of acute heart failure still mention dopamine use as a conceivable alternative mostly due to lack of supporting data against its use [21, 22]. In the recent DanGer trial, close to 25% of all patients with cardiogenic shock used dopamine [23]. Reappraising the SOAP II trial therefore represents an opportunity to both confirm the results for septic shock and provide further evidence for dopamine use in other types of shock, which are still lacking. This analysis therefore was needed to both confirm SOAP II results in the light of newer methods and to provide further evidence for other types of shock.

What those analyses add?

This reanalysis provides an in-depth exploration of SOAP II considering methods that were not disseminated at the time of the original trial’s publication. The results are consistent. All analyses favoured, to some extent, norepinephrine versus dopamine for shock. The findings were mostly inconclusive for patient without septic or cardiogenic shock (“other” types of shock), probably due to the low number of events in this less-severely ill population. Our results highlight that dopamine should not be a drug of choice for any shock type. This is particularly true for patients at higher risk of death (highest predicted APACHE II mortality) where dopamine was associated with fewer days alive and free of ICU. Some other important results include the high probability of harm for dopamine, also considering DAFICU28 as an endpoint, in septic and cardiogenic shock populations. The fact that a simple effect-based model using S-learner was unable to recommend dopamine regardless of the baseline combination of possible modifiers is also a strong argument for not using dopamine as first-line under most circumstances. Considering that it is very unlikely that new trials will test dopamine for shock given the overall context, this reanalysis provide further confirmatory evidence that extends from septic shock.

In addition, we hope that this reanalysis also reinforces that, when done properly, reanalysis can confirm and extend the main results of a trial in a concise way. The number of analyses is purposely high to showcase different approaches that could have been considered had the SOAP II trial been designed more recently.

Limitations

There are several limitations for this analysis. First, this is a post-hoc analysis based on previously published data, and hence should be viewed as exploratory. Second, the SOAP II trial was published more than 10 years ago and therefore, these findings may be limited due to the inevitable changes in ICU practices. Third, the number of analyses included in this analysis is extensive. We aimed to present the results of these analyses as streamlined as possible and provided a summary table of all our findings. Forth, all endpoints are anchored to mortality which was numerically higher in dopamine group in the trial; therefore, analyses using composite endpoints are essentially increasing power due to high number of events; this, however, is common practice in many trials in critical care and other fields. Fifth, definitions of ROPE were arbitrary in the lack of clear consensus on how to define equivalence for endpoints in critical care.

Conclusion

Dopamine may be associated with increased harm when used in patients with shock. This effect seems to be present in patients with either septic or cardiogenic shock. There was no suggestion of any subgroup in which dopamine should be favourable over norepinephrine.

References

  1. De Backer D, Aldecoa C, Njimi H, Vincent JL. Dopamine versus norepinephrine in the treatment of septic shock: a meta-analysis*. Crit Care Med. 2012;40(3):725–30. https://doi.org/10.1097/CCM.0b013e31823778ee.

    Article  PubMed  CAS  Google Scholar 

  2. Rui Q, Jiang Y, Chen M, Zhang N, Yang H, Zhou Y. Dopamine versus norepinephrine in the treatment of cardiogenic shock: a PRISMA-compliant meta-analysis. Medicine. 2017;96(43):e8402. https://doi.org/10.1097/MD.0000000000008402.PMID:29069037;PMCID:PMC5671870.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  3. De Backer D, Biston P, Devriendt J, Madl C, Chochrad D, Aldecoa C, Brasseur A, Defrance P, Gottignies P, Vincent JL, SOAP II Investigators. Comparison of dopamine and norepinephrine in the treatment of shock. N Engl J Med. 2010;362(9):779–89. https://doi.org/10.1056/NEJMoa0907118.

    Article  PubMed  Google Scholar 

  4. Pocock SJ, Ariti CA, Collier TJ, Wang D. The win ratio: a new approach to the analysis of composite endpoints in clinical trials based on clinical priorities. Eur Heart J. 2012;33(2):176–82. https://doi.org/10.1093/eurheartj/ehr352.

    Article  PubMed  Google Scholar 

  5. Zampieri FG, Casey JD, Shankar-Hari M, Harrell FE Jr, Harhay MO. Using bayesian methods to augment the interpretation of critical care trials. An overview of theory and example reanalysis of the alveolar recruitment for acute respiratory distress syndrome trial. Am J Respir Crit Care Med. 2021;203(5):543–52. https://doi.org/10.1164/rccm.202006-2381CP.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Kent DM, Paulus JK, van Klaveren D, D’Agostino R, Goodman S, Hayward R, Ioannidis JPA, Patrick-Lake B, Morton S, Pencina M, Raman G, Ross JS, Selker HP, Varadhan R, Vickers A, Wong JB, Steyerberg EW. The predictive approaches to treatment effect heterogeneity (PATH) statement. Ann Intern Med. 2020;172(1):35–45. https://doi.org/10.7326/M18-3667.

    Article  PubMed  Google Scholar 

  7. Zampieri FG, Damiani LP, Bagshaw SM, Semler MW, Churpek M, Azevedo LCP, Figueiredo RC, Veiga VC, Biondi R, Freitas FR, Machado FR, Cavalcanti AB, BRICNet. Conditional treatment effect analysis of two infusion rates for fluid challenges in critically ill patients: a secondary analysis of balanced solution versus saline in intensive care study (BaSICS) trial. Ann Am Thorac Soc. 2023;20(6):872–9. https://doi.org/10.1513/AnnalsATS.202211-946OC.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Buyse M. Generalized pairwise comparisons of prioritized outcomes in the two-sample problem. Stat Med. 2010;29:3245–57. https://doi.org/10.1002/sim.3923.

    Article  PubMed  Google Scholar 

  9. Nattino G, Finazzi S, Bertolini G, Rossi C, Carrara G. givitiR: the GiViTI calibration test and belt. R package version 1.3, 2017. https://CRAN.R-project.org/package=givitiR. Accessed 30 Mar 2024.

  10. Knaus WA, Draper EA, Wagner DP, Zimmerman JE. APACHE II: a severity of disease classification system. Crit Care Med. 1985;13(10):818–29.

    Article  PubMed  CAS  Google Scholar 

  11. Vincent JL, Moreno R, Takala J, Willatts S, De Mendonça A, Bruining H, Reinhart CK, Suter PM, Thijs LG. The SOFA (sepsis-related organ failure assessment) score to describe organ dysfunction/failure. On behalf of the working group on sepsis-related problems of the European society of intensive care medicine. Intensive Care Med. 1996;22(7):707–10. https://doi.org/10.1007/BF01709751.

    Article  PubMed  CAS  Google Scholar 

  12. Bürkner PC. brms: an R package for bayesian multilevel models using stan. J Stat Softw. 2017;80(1):1–28. https://doi.org/10.18637/jss.v080.i01.

    Article  Google Scholar 

  13. Arel-Bundock V. marginaleffects: predictions, comparisons, slopes, marginal means, and hypothesis tests. R package version 0.18.0, 2024. https://CRAN.R-project.org/package=marginaleffects. Accessed 30 Mar 2024.

  14. Ramanan M, Myburgh J, Finfer S, Bellomo R, Venkatesh B. Publication of secondary analyses from randomized trials in critical care. N Engl J Med. 2016;375(21):2105–6. https://doi.org/10.1056/NEJMc1610367.

    Article  PubMed  Google Scholar 

  15. Gelman A, Carlin J. Beyond power calculations: assessing type S (sign) and type M (magnitude) errors. Perspect Psychol Sci. 2014;9(6):641–51. https://doi.org/10.1177/1745691614551642.

    Article  PubMed  Google Scholar 

  16. Andersen-Ranberg NC, Poulsen LM, Perner A, Hästbacka J, Morgan M, Citerio G, Collet MO, Weber SO, Andreasen AS, Bestle M, Uslu B, Pedersen HS, Nielsen LG, Damgaard K, Jensen TB, Sommer T, Dey N, Mathiesen O, Granholm A. Haloperidol versus placebo for the treatment of delirium in ICU patients: a pre-planned, secondary Bayesian analysis of the AID-ICU trial. Intensive Care Med. 2023;49(4):411–20. https://doi.org/10.1007/s00134-023-07024-9.

    Article  PubMed  CAS  Google Scholar 

  17. Zampieri FG, Serpa-Neto A, Wald R, Bellomo R, Bagshaw SM, STARRT-AKI and RENAL Investigators. Hierarchical endpoints in critical care: a post-hoc exploratory analysis of the standard versus accelerated initiation of renal-replacement therapy in acute kidney injury and the intensity of continuous renal-replacement therapy in critically ill patients trials. J Crit Care. 2024;9(82):154767. https://doi.org/10.1016/j.jcrc.2024.154767.

    Article  CAS  Google Scholar 

  18. de Grooth HJ, Cremer OL. Bayes and the evidence base: reanalyzing trials using many priors does not contribute to consensus. Am J Respir Crit Care Med. 2024;209(5):483–4. https://doi.org/10.1164/rccm.202308-1455VP.

    Article  PubMed  Google Scholar 

  19. Evans L, Rhodes A, Alhazzani W, Antonelli M, Coopersmith CM, French C, Machado FR, Mcintyre L, Ostermann M, Prescott HC, Schorr C, Simpson S, Wiersinga WJ, Alshamsi F, Angus DC, Arabi Y, Azevedo L, Beale R, Beilman G, Belley-Cote E, Burry L, Cecconi M, Centofanti J, Coz Yataco A, De Waele J, Dellinger RP, Doi K, Du B, Estenssoro E, Ferrer R, Gomersall C, Hodgson C, Hylander Møller M, Iwashyna T, Jacob S, Kleinpell R, Klompas M, Koh Y, Kumar A, Kwizera A, Lobo S, Masur H, McGloughlin S, Mehta S, Mehta Y, Mer M, Nunnally M, Oczkowski S, Osborn T, Papathanassoglou E, Perner A, Puskarich M, Roberts J, Schweickert W, Seckel M, Sevransky J, Sprung CL, Welte T, Zimmerman J, Levy M. Surviving sepsis campaign: international guidelines for management of sepsis and septic shock 2021. Crit Care Med. 2021;49(11):e1063–143. https://doi.org/10.1097/CCM.0000000000005337.

    Article  PubMed  Google Scholar 

  20. Patel GP, Grahe JS, Sperry M, Singla S, Elpern E, Lateef O, Balk RA. Efficacy and safety of dopamine versus norepinephrine in the management of septic shock. Shock. 2010;33(4):375–80. https://doi.org/10.1097/SHK.0b013e3181c6ba6f.

    Article  PubMed  CAS  Google Scholar 

  21. McDonagh TA, Metra M, Adamo M, Gardner RS, Baumbach A, Böhm M, Burri H, Butler J, Čelutkienė J, Chioncel O, Cleland JGF, Coats AJS, Crespo-Leiro MG, Farmakis D, Gilard M, Heymans S, Hoes AW, Jaarsma T, Jankowska EA, Lainscak M, Lam CSP, Lyon AR, McMurray JJV, Mebazaa A, Mindham R, Muneretto C, Francesco Piepoli M, Price S, Rosano GMC, Ruschitzka F, Kathrine SA. Corrigendum to: 2021 ESC guidelines for the diagnosis and treatment of acute and chronic heart failure: developed by the task force for the diagnosis and treatment of acute and chronic heart failure of the European Society of Cardiology (ESC) With the special contribution of the Heart Failure Association (HFA) of the ESC. Eur Heart J. 2021;42(48):4901. https://doi.org/10.1093/eurheartj/ehab670. (Erratum for: Eur Heart J. 2021 Sep 21;42(36):3599–3726).

    Article  PubMed  Google Scholar 

  22. Heidenreich PA, Bozkurt B, Aguilar D, Allen LA, Byun JJ, Colvin MM, Deswal A, Drazner MH, Dunlay SM, Evers LR, Fang JC, Fedson SE, Fonarow GC, Hayek SS, Hernandez AF, Khazanie P, Kittleson MM, Lee CS, Link MS, Milano CA, Nnacheta LC, Sandhu AT, Stevenson LW, Vardeny O, Vest AR, Yancy CW. 2022 AHA/ACC/HFSA guideline for the management of heart failure: a report of the American College of Cardiology/American Heart Association Joint Committee on Clinical Practice Guidelines. Circulation. 2022;145(18):e895–1032. https://doi.org/10.1161/CIR.0000000000001063. (Erratum in: Circulation. 2022 May 3;145(18):e1033. Erratum in: Circulation. 2022 Sep 27;146(13):e185. Erratum in: Circulation. 2023 Apr 4;147(14):e674).

    Article  PubMed  Google Scholar 

  23. Møller JE, Engstrøm T, Jensen LO, Eiskjær H, Mangner N, Polzin A, Schulze PC, Skurk C, Nordbeck P, Clemmensen P, Panoulas V, Zimmer S, Schäfer A, Werner N, Frydland M, Holmvang L, Kjærgaard J, Sørensen R, Lønborg J, Lindholm MG, Udesen NLJ, Junker A, Schmidt H, Terkelsen CJ, Christensen S, Christiansen EH, Linke A, Woitek FJ, Westenfeld R, Möbius-Winkler S, Wachtell K, Ravn HB, Lassen JF, Boesgaard S, Gerke O, Hassager C. Microaxial flow pump or standard care in infarct-related cardiogenic shock. N Engl J Med. 2024;390(15):1382–93. https://doi.org/10.1056/NEJMoa2312572.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

SOAP II investigators: P. Biston, J. Devriendt, C. Madl, D. Chochrad, C. Aldecoa, A. Brasseur, P. Defrance, P. Gottignies, R. Kitzberger, U. Holzinger, A. Roman, D. De Bels, S. Anane, S. Brimioulle, M. Van Nuffelen, M. VanCutsem, J. Rico, J.I. Gomez Herreras, C. Mélot.

Funding

There was no specific funding for this analysis.

Author information

Authors and Affiliations

Authors

Consortia

Contributions

Conception: FGZ, SMB, JLV, DDD. Statistical analysis: FGZ, HJ. Data collection: JLV, DDD Manuscript draft: FGZ. Revision for relevant intelectual content: SMB, HJ, JLV, DDD. Data collection: SOAP II investigator.

Corresponding author

Correspondence to Fernando G. Zampieri.

Ethics declarations

Ethics approval and consent to participate

This analysis was made based on fully anonymized secondary analysis of SOAP II trial provided by the investigators; no new data except from the data already available from the trial was used. SOAP II study was approved by local ethics committees at enrolling sites.

Competing interests

Prof Jean-Louis Vincent is Critical Care´s journal editor. FGZ has received personal consulting fees from Baxter and Bactiguard (Sweden), unrelated to the scope of this article.

Additional information

Publisher's Note

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

The members of the “SOAP II investigators” group are listed in the acknowledgements section.

Supplementary Information

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zampieri, F.G., Bagshaw, S.M., Njimi, H. et al. Exploration of different statistical approaches in the comparison of dopamine and norepinephrine in the treatment of shock: SOAP II. Crit Care 28, 299 (2024). https://doi.org/10.1186/s13054-024-05016-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13054-024-05016-9

Keywords