Exploring the efficacy of psychotherapies for depression: a multiverse meta-analysis
====================================================================================

* Constantin Yves Plessen
* Eirini Karyotaki
* Clara Miguel
* Marketa Ciharova
* Pim Cuijpers

## Abstract

**Background** Hundreds of randomised controlled trials and dozens of meta-analyses have examined psychotherapies for depression—yet not all points in the same direction. Are these discrepancies a result of specific meta-analytical decisions or do most analytical strategies reaching the same conclusion?

**Objective** We aim to solve these discrepancies by conducting a multiverse meta-analysis containing all possible meta-analyses, using all statistical methods.

**Study selection and analysis** We searched four bibliographical databases (PubMed, EMBASE, PsycINFO and Cochrane Register of Controlled Trials), including studies published until 1 January 2022. We included all randomised controlled trials comparing psychotherapies with control conditions without restricting the type of psychotherapy, target group, intervention format, control condition and diagnosis. We defined all possible meta-analyses emerging from combinations of these inclusion criteria and estimated the resulting pooled effect sizes with fixed-effect, random-effects, 3-level, robust variance estimation, *p*-uniform and PET-PEESE (precision-effect test and precision-effect estimate with SE) meta-analysis models. This study was preregistered ([https://doi.org/10.1136/bmjopen-2021-050197](https://doi.org/10.1136/bmjopen-2021-050197)).

**Findings** A total of 21 563 records were screened, and 3584 full texts were retrieved; 415 studies met our inclusion criteria containing 1206 effect sizes and 71 454 participants. Based on all possible combinations between inclusion criteria and meta-analytical methods, we calculated 4281 meta-analyses. The average summary effect size for these meta-analyses was Hedges’ *g* mean=0.56, a medium effect size, and ranged from *g*=−0.66 to 2.51. In total, 90% of these meta-analyses reached a clinically relevant magnitude.

**Conclusions and Clinical Implications** The multiverse meta-analysis revealed the overall robustness of the effectiveness of psychotherapies for depression. Notably, meta-analyses that included studies with a high risk of bias, compared the intervention with wait-list control groups, and not correcting for publication bias produced larger effect sizes.

*   Depression & mood disorders

#### WHAT IS ALREADY KNOWN ON THIS TOPIC

*   To this day, more than 80 meta-analyses have examined the efficacy of psychotherapies for depression. However, recent meta-research projects questioned whether the effect sizes from those published meta-analyses were inflated and discrepancies between these meta-analyses emerged.

#### WHAT THIS STUDY ADDS

*   In a so-called multiverse meta-analysis, we calculated over 4000 meta-analyses—most of them (90%) produced small but clinically relevant effect sizes. Our findings suggest that psychotherapies for depression are generally effective, but the specific type, format and other factors can affect the magnitude of the treatment effect slightly. However, meta-analyses that (1) restrict their control group to wait-list control groups, (2) do not exclude high risk of bias studies and (3) do not correct for publication bias are likely to produce larger effect sizes.

#### HOW THIS STUDY MIGHT AFFECT RESEARCH, PRACTICE OR POLICY

*   The effect size reported in any given meta-analysis on treatment efficacy for depression depends less on the type of psychotherapy, treatment format, diagnosis or target group but rather on the comparison with wait-list control groups, not excluding the high risk of bias studies, or correcting for publication bias. In general, future meta-analyses that diverge from small-to-medium summary effect size estimates may be indicative of extreme data analytical decisions and may therefore not contribute additional substantive knowledge. However, there are certain circumstances where higher treatment effects might be observed, such as when a future meta-analysis examines a specific subgroup of patients or a new type of therapy not previously studied.

## Background

Over the last four decades, more than 80 meta-analyses have examined the efficacy of psychotherapies for depression.1 In these meta-analyses, evidence from more than 700 randomised controlled trials (RCTs) is included, yet not all of these studies are pointing in the same direction.2 3 Contested evidence exists on efficacy claims between different psychotherapies for depression (eg, therapies based on cognitive–behavioural therapy (CBT) or other types of psychotherapy),2–12 target groups (eg, adults or general medical populations) and delivery formats (eg, individual or group therapy).

Some of the discrepancies in findings may be the result of publication bias leading to an overestimation of the effectiveness of psychotherapy or may be due to variations in inclusion criteria, such as the inclusion of low-quality studies or studies comparing interventions with wait-list control groups only.13–18 It is crucial to evaluate the influence such meta-analytical decisions have. For example, does it make a substantial difference when we correct for publication bias or not? Does the evidence depend on whether we include only the best evidence or all evidence? Are the results robust to slightly different inclusion criteria? This exploration is especially important when multiple meta-analyses with overlapping research questions reach different conclusions.19 20 To increase trust in the existing evidence, we need to ensure that the published results do not depend on these specific decisions in selecting and analysing the data but rather that most analytical strategies reach the same conclusion.

Although conventional meta-analyses exist on some of these specific aspects, a comprehensive bird’s-eye view of all meta-analyses for depression treatment research is missing. Ideally, this birds-eye view does not only include all published meta-analyses but also all possible meta-analyses based on defensible and reasonable analytical choices. To provide such an overview and fill substantive knowledge gaps, we conducted a so-called multiverse meta-analysis and calculated all possible meta-analyses on the efficacy of psychotherapies for depression in a single analysis. It can (1) integrate multiple meta-analyses like an umbrella review, (2) enable us to identify knowledge gaps and (3) investigate how flexibility in data selection and analysis might affect the overall interpretation of results. In doing so, this new approach can help solve diverging claims on the efficacy of psychotherapies once and for all. We replicated most meta-analyses that have ever been conducted in research on psychotherapy for depression and created additional evidence by conducting thousands of meta-analyses that were still missing in the literature.

Due to the sheer number of published meta-analyses and primary studies on these differences between psychotherapies, we aimed to summarise, integrate and visualise the entire evidence. As a result, we can highlight robustness—or lack thereof—by inspecting all possible meta-analyses to help resolve conflicting meta-analyses and contested evidence, alleviate the associated adverse effects of these phenomena on research progress and provide a birds-eye perspective of the entire field.21 

## Study selection and analysis

### Search strategy and selection criteria

We searched four major bibliographic sources (PubMed, PsycINFO, EMBASE and Cochrane Library; see online supplemental eMethods 1 for all search strings) for RCTs of psychotherapies for depression published until 1 January 2022.1 22 After title and abstract screening, two independent researchers conducted full-text screening of all records. We included all RCTs comparing a psychological intervention with any control condition written in English, German, Spanish or Dutch. We excluded maintenance and relapse prevention trials, dissertations and interventions not aimed at depression. Eligible were both self-reported and clinician-rated instruments measuring depression. Therapies could be delivered by any person trained to deliver the therapy. Two independent researchers extracted information on target groups, intervention formats, psychotherapy types, control conditions and countries. Inconsistencies were resolved by discussion. This study was preregistered ([https://doi.org/10.1136/bmjopen-2021-050197](https://doi.org/10.1136/bmjopen-2021-050197)).

### Supplementary data

[[bmjment-2022-300626supp001.pdf]](pending:yes) 

### Data analysis

Our study protocol outlines the analyses in more detail (see online supplemental eMethods 3 for deviations from our protocol).23 The R code and data to reproduce all analyses can be found at the Open Science Framework ([https://osf.io/mtx8a/](https://osf.io/mtx8a/)).24 All analyses were carried out using R (V.4.1.2)25 and the metafor package (V3.4.0).26 We calculated standardised mean differences (Hedges’ *g*) for postintervention comparisons between psychotherapy and control conditions based on continuous outcome data provided in the primary studies. If only change score or dichotomous outcome data were reported, we converted these data into Hedges’ *g* with the metapsyTools R package.27 We assessed the risk of bias of included studies using four criteria of the risk of bias assessment tool, developed by the Cochrane Collaboration28: adequate generation of allocation sequence, concealment of allocation to conditions, masking of assessors and dealing with incomplete outcome data (this criterion was met when intention-to-treat analyses were conducted). All items were rated as positive (the criterion was met) or negative (the criterion was not met or unclear). The total risk of bias score for each study was calculated as the sum of all positive scores (ranging from 0 to 4, with 4 indicating no risk of bias). We rated a study having overall ‘some concern’ for risk of bias when the study had a rating of 1 or higher. Two researchers conducted the risk of bias independently, and disagreements were resolved through discussion.

#### Multiverse meta-analysis

Our multiverse meta-analysis contains every single meta-analysis based on a defensible combination of subgroups (eg, target group of the intervention, type and format of intervention) and statistical methods investigating psychotherapies targeting depression that could possibly be conducted.

##### Descriptive specification curve

We specified seven the so-called *Which* factors—asking which data to meta-analyse—and one *How* factors—asking how to meta-analyse the data. Based on these *Which* factors, meta-analyses could include studies investigating different *target groups* that received different *types of psychotherapies* in different *formats* and *assessments/diagnoses* of depression. These studies could have different *risks of biases’* ratings and could be compared with different *control groups* (see online supplemental eMethods 2 for a detailed description of all *Which* factors). Based on our *How* factor, we used six different meta-analytical methods to analyse the data: random-effects, fixed-effect, 3-level, robust variance estimation (RVE),29 precision-effect test and precision-effect estimate with SEs (PET-PEESE)30 and *p*-uniform* meta-analysis models.31 

Some primary studies reported multiple effect sizes per study, that is, when multiple instruments were used to measure depression or when a study consisted of multiple interventions or control groups. These nested effect sizes are not independent as they are correlated (we assumed a correlation of *r*=0.5) and introduce a unit-of-analysis problem.32 33 Such effect size dependencies were handled by either averaging the effect sizes included in each study or modelling the dependency directly. For this hierarchical modelling of dependencies among effect sizes within studies, we included 3-level and RVE models. RVE methods allow for the inclusion of all relevant effect sizes in a single meta-regression model, regardless of the specific nature of the dependencies between them.

Additionally, we used two novel methods for addressing ‘small-study effects’, which can arise due to publication bias (eg, when statistically non-significant findings are not published), reporting bias (eg, when statistically significant results are selectively reported) or clinical heterogeneity (eg, when smaller studies include more severely ill patients than larger studies). This phenomenon, which is common in many scientific fields, can cause inflated effect sizes and therefore an overestimation of treatment effectiveness.34 

PET-PEESE is a regression-based method for addressing the issue of *small-study effects* in meta-analyses.30 34 It is based on the relationship between effect sizes and SEs and is part of a broader class of funnel plot–based methods, such as the trim-and-fill method or Egger’s regression test.35 36 In the absence of publication bias and reporting bias, the relationship between effect sizes and SEs should be unrelated. However, publication bias often results in a disproportionate number of larger studies being published, while smaller studies are only published if they show statistically significant results. This can lead to an over-representation of imprecise studies with inflated effect size estimates in the published literature.

PET-PEESE aims to correct for this bias and provides more accurate estimates of effect sizes.

*P*-uniform* is a selection model approach that uses a random-effects model as its effect size model.31 This method assumes that the probability of publishing a statistically significant or non-significant effect size is constant, but these probabilities may be different from each other. *P*-uniform* works by treating the primary studies’ effect sizes differently based on whether they are statistically significant or not. This method can be considered a selection model approach with a single cut-off value that determines whether an effect size is considered statistically significant.

The results of the multiverse meta-analysis were depicted in a descriptive specification curve plot, which is a graph that shows the results of all conducted meta-analyses, represented by points on the graph with 95% CIs. The summary effects are plotted in order of magnitude, from lower to higher, and connected by a solid line—the specification curve. To evaluate the number of meta-analyses that exceed relevant magnitudes, the plot focuses on two thresholds: meta-analyses that do not include 0 in their 95% CI (ie, meta-analyses that exceed a null effect) and meta-analyses that exceed a clinically relevant effect size of Hedges’ *g*=0.24, which was suggested as the minimal important difference for interventions targeting major depressive disorder.37 The plot includes vertical columns that represent different combinations of factors that may influence the meta-analyses, such as target group, therapy type and control group.

##### Inferential specification curve

In a second step, we evaluated if the findings of the descriptive specification curve (the magnitude sorted summary effect sizes from all meta-analyses of the multiverse) are likely to be true or if they could be due to chance. This so-called inferential specification curve analysis involves simulating new random-effect sizes for each primary study under the assumption that the null hypothesis, or the assumption that there is no psychological treatment effect on depression, is true.

These new data sets are created by drawing random values for the effect sizes from a normal distribution with mean zero and an SD that takes into account both the variance of the original effect size estimate for each study and a measure of between-study heterogeneity (τ values obtained from fitting a random-effects model with REML estimator and handling effect size dependency by averaging studies with multiple effect sizes on the entire data set).

We then applied a new descriptive specification curve analysis under three scenarios: a fixed-effect scenario with no heterogeneity, a scenario with heterogeneity equal to the random-effects model of the 415 studies and a scenario with the upper 95% CI estimate of τ from the random-effects model. By repeating this process 1000 times, we were able to identify the lower and upper limits of the resulting 1000 bootstrapped specification curves or their 2.5% and 97.5% quantiles. These limits, or quantiles, represent the range within which we would expect 95% of all meta-analyses to fall if the true underlying effect were a null effect. If the observed treatment effect falls outside of these limits, it is considered to be a deviation from the null hypothesis and is considered to be a likely true effect.

Additionally, we conducted a GOSH plot (Graphical Display of Study Heterogeneity) to identify the overall range of possible summary effect sizes for meta-analyses and visualise the heterogeneity of effect sizes.38 This plot is a visual tool that can be thought of as a brute force sensitivity check because it calculates all possible meta-analyses from all possible subsets of included studies. It is not restricted to the more theoretically guided comparisons defined by the descriptive specification curve analysis. The GOSH plot shows the relationship between the effect size of each meta-analysis and its heterogeneity, which is a measure of the variability of effect sizes across studies. We used a reduced set of 100 000 randomly drawn samples as this sensitivity check is computationally infeasible with 415 primary studies, as 2415=8.46×10124.

##### Conventional meta-analysis

To create a reference point for exploring heterogeneity in the data, we additionally fitted a 3-level model to the entire data set. This resulting summary effect size represented *one* possible meta-analysis—including the broadest inclusion of all *Which* factors and the 3-level *How* factor—out of the entire multiverse of meta-analyses. The amount of heterogeneity was estimated using the restricted maximum likelihood estimator. In addition to the estimate of ![Formula][1]</img> , both the Q-test for heterogeneity and the ![Formula][2]</img>  statistic were reported.

## Findings

We screened 21 563 titles and examined 3584 full-text papers, of which 415 were included. The PRISMA (Preferred Reporting Items for Systematic Reviews and Meta-Analyses) flow chart depicts the study inclusion process in figure 1. Overall, *kes *=1206 effect sizes from *kstudies *=415 studies were included. The sample sizes of the included primary studies ranged from n=4 to 1156, *Nmean *=103, *Nmedian *=67. The total sample size of all included samples from all primary studies was *Ntotal *=71 454. See online supplemental eTable 1 for a detailed description of all included primary studies.

![Figure 1](http://mentalhealth.bmj.com//https://mentalhealth.bmj.com/content/ebmental/26/1/e300626/F1.medium.gif)

[Figure 1](http://mentalhealth.bmj.com//content/26/1/e300626/F1)

Figure 1 
PRISMA (Preferred Reporting Items for Systematic Reviews and Meta-Analyses) flow chart for inclusion of studies.

Most included primary studies investigated the efficacy of psychotherapies for depression in adults (37%) or general medical populations (23%). Most studies were either conducted in Europe (37%) or in Northern America (36%). Most studies used CBT-based interventions (69%) and were primarily delivered in individual therapy format (35%). Care-as-usual was the most common type of control condition (49%), and depressive disorder was diagnosed by a clinician in 49% of studies. Only 35% of studies were rated with a low risk of bias. See table 1 for the study characteristics and online supplemental eTable 2 for the effect size characteristics.

View this table:
[Table 1](http://mentalhealth.bmj.com//content/26/1/e300626/T1)

Table 1 
Summary characteristics of included primary studies

Our multiverse meta-analysis produced 4281 unique meta-analyses, with effect sizes ranging from Hedges’ *g*=−0.66 to 2.51. Half of those effect sizes were in the IQR of Hedges’ *g*=0.42 to 0.71, representing small-to-medium effect sizes. In total, 97% of the effect sizes were greater than 0, and 84% of these had 95% CIs that did not include 0 (ie, estimated *g* was greater than 0 which would have returned a two-tailed p-value of less than 0.05). In total, 90% reached a clinically relevant magnitude of Hedges’ *g* >0.24, and 68% of the summary effect sizes had 95% CIs above the clinically relevant cut-off.

The overall pattern of the descriptive specification curve indicates that larger meta-analyses, including more primary studies, had medium-to-large effect sizes and were close to the median-estimated effect size of the multiverse (see figure 2). More extreme meta-analytical effect sizes were associated with few included studies and therefore broader CIs.

![Figure 2](http://mentalhealth.bmj.com//https://mentalhealth.bmj.com/content/ebmental/26/1/e300626/F2.medium.gif)

[Figure 2](http://mentalhealth.bmj.com//content/26/1/e300626/F2)

Figure 2 
Descriptive specification curve: psychotherapies for depression. The top panel shows the outcome of all 4281 meta-analyses (Hedges’ *g*) with their 95% CIs. The summary effects are sorted by magnitude, from lower to higher. Connecting the different summary effects results in the solid line, which is the specification curve. A horizontal dashed line of no effect is shown at *g*=0, and a red dotted line indicates a clinically relevant effect size at *g*=0.24. The vertical columns in the bottom panel represent factor combinations of *How* factors (different target groups, formats, therapy types, control groups and diagnoses) and *Which* factors (3-LVL, 3-level model; FEM, fixed-effect model; PET-PEESE, *p*-uniform*; REML, random-effects model; RVE, robust variance estimation) constitute a given specification. The location of each *Which* factors’ largest meta-analysis—containing information from all 415 studies—is depicted on the specification curve. Each vertical column is colour-coded, signifying the number of samples included in a specification (hot spectral colours for more included samples vs cool spectral colours for less included samples). CBT, cognitive–behavioural therapy; PET-PEESE, precision-effect test and precision-effect estimate with SE; ROB, risk of bias.

Several *Which* and *How* factors produced—on average—systematically different summary effect size estimates compared with others. In the following, we descriptively summarise the most important results. It is important to note that this examination of these differences is based on descriptive analysis rather than a formal statistical comparison. For a more detailed breakdown of each *Which* factor, see online supplemental eFigures 1–7 and online supplemental eTables 3–9.

### Target group

Meta-analyses including only student populations, mean *g*=0.82, 95% CI (0.51, 1.12) from *k*=54 included meta-analyses, produced larger effect size estimates than meta-analyses on adults, mean *g*=0.51, 95% CI (0.26, 0.75) with *k*=1076.

### Format

Meta-analyses including studies delivered in a group format, mean *g*=0.76, 95% CI (0.32, 1.19) with *k*=536, produced larger effect size estimates than meta-analyses delivered as guided self-help interventions, mean *g*=0.50, 95% CI (0.27, 0.72) with *k*=586.

### Type

Meta-analyses focusing only on CBT-based treatments, mean *g*=0.60, 95% CI (0.29, 0.91) with *k*=1620, produced larger effect size estimates than meta-analyses focusing on non-CBT treatments, mean *g*=0.48, 95% CI (0.19, 0.78) with *k*=733.

### Control group

Meta-analyses that included samples compared with a wait-list control group, mean *g*=0.66, 95% CI (0.35, 0.96) with *k*=836, produced larger effect size estimates than treatments compared with care-as-usual, mean *g*=0.52, 95% CI (0.22, 0.82) with *k*=1194.

### Risk of bias

Meta-analyses that excluded high risk of bias studies, mean *g*=0.61, 95% CI (0.27, 0.95) with *k*=2413 included samples, produced larger effect size estimates than meta-analyses including only low risk of bias studies, mean *g*=0.45, 95% CI (0.19, 0.72) with *k*=1034.

### Diagnosis

Meta-analyses that included studies in which depression was diagnosed by a clinician or was self-reported produced similar results, mean *g*=0.56 and 0.57, respectively.

### Meta-analytical method

Meta-analyses analysed with 3-level models, mean *g*=0.66, 95% CI (0.36, 0.96) with *k*=963, produced larger effect size estimates than meta-analyses analysed with PET-PEESE, mean *g*=0.18, 95% CI (−0.24, 0.59) with *k*=591.

We further investigated the studies in which there was no strong evidence that psychotherapies were effective indicated by including zero in the 95% CI. On closer inspection, we found that the observed null effects were largely caused by using different *How* factors rather than different *Which* factors. Of the 688 meta-analyses that included a zero in their 95% CI, the PET-PEESE method accounted for 408, the RVE estimation method for 143 (as it tends to produce larger CIs), *p*-uniform for 68 and 3-level modelling for 59. The fixed-effect (only one meta-analysis) and REML models (nine meta-analyses) barely produced any statistically non-significant meta-analyses. No other systematic differences were observed among those meta-analyses.

The results of the inferential specification curve analysis shown in figure 3 indicate that in most cases, treatments for depression had a substantial effect, as indicated by the deviation from the scenario of no effect (*g*=0). This was true for a range of scenarios, from simulating no heterogeneity (*τ*=0) to the identified heterogeneity in the random-effects model of all 415 studies (*τ*=0.53) and the respective 95% CI upper limit of *τ*=0.71. In fact, 98% of the meta-analyses were outside the expected range under the scenario of no effect, suggesting that the meta-analyses found in the multiverse meta-analysis were significantly different from a null effect. In addition, we conducted a similar analysis for a scenario in which the simulated studies had a clinically relevant effect size of Hedges’ *g*=0.24, rather than *g*=0. We found that, consistent with our previous analysis, most treatments were more effective than the simulated effect sizes in this scenario (see online supplemental eFigure 8).

![Figure 3](http://mentalhealth.bmj.com//https://mentalhealth.bmj.com/content/ebmental/26/1/e300626/F3.medium.gif)

[Figure 3](http://mentalhealth.bmj.com//content/26/1/e300626/F3)

Figure 3 
Inferential specification curve. Depicted are three inferential specification curves (red lines), each summarising the magnitude-sorted meta-analytic summary effects and correspond to the descriptive specification curve from figure 2 (for computational reasons only the REML, FE and PET-PEESE models were used). The grey area represents the corresponding 97.5% and 2.5% quantiles of 1000 specification curves that were simulated assuming no true effect. The left panel depicts a fixed-effect scenario of no heterogeneity (*τ*=0), the middle panel a scenario of heterogeneity equal to the random-effects model (*τ*=0.53) and the right panel the scenario of the upper 95% CI estimate of *τ* from the random-effects model (*τ*=0.71). Each is simulated under the null hypothesis for a given specification number using a parametric bootstrap procedure, but they differ in underlying heterogeneity assumptions. If the specification curve exceeds the limits of the 95% CI (as is the case in this plot), there is evidence against the null hypothesis (*g*=0), indicating that there is a substantial effect for the effectiveness of psychotherapies for depression. FE, fixed effect; PET-PEESE, precision-effect test and precision-effect estimate with SE; REM, random-effects model.

The GOSH plot (see figure 4) from 100 000 random samples from all possible subset combinations of 1206 included effect sizes revealed a similar picture as the descriptive specification curve: most meta-analyses fell in the IQR 0.43–0.7 of the multiverse meta-analysis, while heterogeneity is substantial.

![Figure 4](http://mentalhealth.bmj.com//https://mentalhealth.bmj.com/content/ebmental/26/1/e300626/F4.medium.gif)

[Figure 4](http://mentalhealth.bmj.com//content/26/1/e300626/F4)

Figure 4 
GOSH (Graphical Display of Study Heterogeneity) plot. This GOSH plot visualises the heterogeneity of a random sample of 100 000 subsets for the combinatorial meta-analysis. The y-axis depicts Higgins *I* 2 statistics for heterogeneity, and the summary effect size is visualised on the x-axis. Density distributions are visualised next to the respective axes.

We chose to include a conventional 3-level meta-analysis, including all 1206 effect sizes from all 415 primary studies, as an exemplary meta-analytical specification out of the 4281 meta-analyses. The pooled Hedges’ *g* based on this 3-level meta-analytic model was *g*=0.72, 95% CI (0.66, 0.78), p<0.001). The estimated variance components were *τ* 2 Level 3=0.364 and *τ* 2 Level 2=0.004. Overall, *I* 2 Level 3=86% of the total variation can be attributed to between-study and *I* 2 Level 2=1% to within-study heterogeneity. We found that the 3-level model provided a significantly better fit compared with a 2-level model with level 3 heterogeneity constrained to zero (*χ* 2 * 1 *=654.03, p<0.001). The overall *I* 2 value, indicating how much of the total variance can be attributed to the total amount of heterogeneity, is very large, with approximately 87% of the total variance attributable to heterogeneity.

## Conclusions and clinical implications

We investigated the efficacy of psychotherapies for depression by simultaneously analysing 4281 meta-analyses based on all reasonable combinations of inclusion criteria and meta-analytical methods. We investigated the influence of different treatment groups, types of psychotherapy, treatment formats, control groups, diagnoses for depression, risk of bias assessments and different meta-analytical methods. We found that most meta-analyses produced small-to-medium, but clinically relevant effect sizes, suggesting the overall robustness of psychotherapies for depression. These results are supported by our descriptive and inferential multiverse meta-analysis and the combinatorial meta-analytical approach. Notably, specific patterns emerged in the descriptive specification curve analysis.

Meta-analyses that compared interventions with wait-list control groups had larger effect sizes compared with those including care-as-usual or other control conditions, such as attention placebo. This aligns with previous research suggesting that effect sizes obtained in RCTs should be interpreted differently depending on the used control condition and warning especially against lumping control conditions into one comparison group for network meta-analyses.14 It is likely that the larger effect sizes found in meta-analyses that compare interventions with wait-list control groups are due to the fact that wait-list control groups are considered to be less effective treatments than other types of control conditions, such as treatment-as-usual. This means that when interventions are compared with wait-list control groups, the difference in outcome is likely to be greater, resulting in a larger effect size estimate. In other words, the comparison group in these meta-analyses are not as good as other control conditions and therefore the difference with the intervention is greater, which results in a larger effect size.12 39 Therefore, it is important to keep in mind that this increase in effect size estimates should be considered when interpreting studies that used wait-list control groups as a comparison.

Our findings indicate that the effect sizes found in depression research might be inflated, even when high risk of bias studies are removed. To support this conclusion, we compared the results of meta-analyses that excluded high risk of bias studies with meta-analyses that did not exclude them. We found that excluding high risk of bias studies did not result in a reduction in effect size estimates. In fact, the effect sizes were very similar to those observed when all studies, regardless of study quality, were included in the meta-analysis. However, when we only included studies with a low risk of bias in our meta-analyses, we did observe a reduction in effect size estimates. This finding is in line with earlier research that has also suggested that effect sizes in depression research may be inflated, as including only low risk of bias studies yields substantially smaller effect size estimates.13 

Our analysis of publication bias (small-study effects) in psychotherapy research on depression showed that one of the methods used to correct publication bias resulted in substantially smaller effect size estimates compared with conventional methods that did not correct for publication bias. Even though the PET-PEESE estimator may have overcorrected for biases in our data, resulting in negative effect sizes, these negative effect sizes indicate that the meta-analyses under study had no effect when correcting for small-sample effects.40 This finding suggests that the previously reported effect sizes for psychotherapy in the treatment of depression may be overestimated and potentially inflated. This is consistent with previous research, which has found significant publication bias for psychotherapies of major depressive disorder,41 for digital psychological interventions for depression42 and for most evidence-supported therapies for adult depression.2 In other words, the true effect size of psychotherapy in treating depression may be smaller than what has been previously reported in the literature.

In addition, we found that meta-analyses focused on student populations and interventions delivered in a group format had higher effect size estimates. This finding is consistent with previous meta-analyses that found interventions to be more effective in young adults compared with middle-aged adults,43 although individual formats tend to be at least as effective as group formats.44 

Our multiverse meta-analysis was able to demonstrate that the inclusion of high and medium risk of bias studies, the comparison with wait-list control groups and models not accounting for publication bias yielded larger—and potentially inflated—effect size estimates for the efficacy of psychotherapies for depression. At the same time, we were able to demonstrate that even after considering this inflation, the effectiveness of psychotherapies remains clinically significant. Despite these strengths of our study, several limitations need to be acknowledged. First, our database contains several studies with unreasonably large effect sizes (Hedges’ *g*>3). These enormous effect sizes might distort meta-analyses with only a few included studies towards more extreme summary effect sizes. For this reason, we included only meta-analyses with more than 10 studies in our multiverse meta-analysis to avoid such extreme meta-analyses. We assessed the influence of several cut-offs in online supplemental eTable 10 and online supplemental eFigure 9. Our sensitivity analyses indicate that the overall mean summary effect size does not change for multiverse meta-analyses limited to at least 10, 25 or 50 primary studies, yet the spread of possible summary effect sizes changes substantially. Second, we had to slightly deviate from our preregistered protocol as we merged several *Which* factor categories to ensure that both computations and visualisations remained feasible. For instance, we initially planned to investigate each therapy type individually, yet we had to create broader categories and combined CBT-based and non-CBT-based therapy approaches together. These merges decreased the level of detail of our analyses, but at the same time ensured the interpretability, as a visualisation of over 40 000 meta-analyses (resulting from not merging different therapy types) was simply not possible. Third, the methods we used for correcting for publication bias have some limitations in performing in environments with low sample size and high heterogeneity—as was the case with the body of evidence in this multiverse meta-analysis. This might have caused the PET-PEESE method to underestimate and *p*-uniform to overestimate the effect size. Overall, the presence of publication bias and its inflating influence on effect size estimates remain highly likely and have been reported in other publications as well.12 18 45 It is important to note that the investigated small-study effects can arise from different biases (publication bias, reporting bias), but also might be indicative of a genuine effect. This can, for example, be the case when smaller studies are performed in different populations (ie, difficult to research, high disease burden) or different clinical settings, where the effect is genuinely larger than in larger trials.

Most studies in the present literature only consider one or two of the presented *Which* factors when evaluating the efficacy of psychological treatments for depression. To ensure the comparability of our results with such effect sizes, we also assessed the effect of each ‘Which’ factor separately in our analyses. Presumably, because not all *Which* factors were considered simultaneously as predictors of depressive symptoms, as it is done, for example, in multiple linear regression analyses, effect sizes for each separate factor might have been overestimated. This is true for both the included primary studies and meta-analyses that took such an approach, as well as for the derived effect sizes of our multiverse meta-analysis. This fact, however, again underpins one of the main claims of our study, namely those future meta-analyses that diverge from small-to-medium summary effect size estimates are likely indicative of extreme data analytical decisions and should be evaluated with great care.

In summary, our multiverse meta-analysis completes the evidence that psychotherapies are effective for treating depression in a wide range of patient populations. Because we evaluated the entire multiverse of defensible combinations of inclusion criteria and statistical methods on depression research, our results suggest that this line of research and the debate about whether treatment is effective can now end once and for all. Future research can and should be less concerned with *whether* therapies work but rather investigate *how* they work and who benefits most from which type of intervention. New approaches like individual-patient data meta-analyses and (component) network meta-analyses, as well as longitudinal approaches, are needed to investigate these more relevant and critical issues for the individual patient. Finally, this study provides the most exhaustive overview of psychological depression research that is available and possible today. It can guide future research as knowledge gaps were closed and is a valuable source for policy-makers to inform evidence-based decision-making.

## Data availability statement

Data are available in a public, open access repository. The R code and data to reproduce all analyses can be found at the Open Science Framework.

## Footnotes

*   Twitter @CYPlessen

*   Contributors CYP conceived of the presented idea and is responsible for the overall content as guarantor. PC, CM, MC and EK prepared the data curation. CYP analysed the data and wrote the initial draft. The draft was reviewed and edited by CM, MC, PC and EK. PC and EK supervised the findings of this work.

*   Funding The authors have not declared a specific grant for this research from any funding agency in the public, commercial or not-for-profit sectors.

*   Competing interests None declared.

*   Provenance and peer review Not commissioned; externally peer reviewed.

*   Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.

[https://creativecommons.org/licenses/by/4.0/](https://creativecommons.org/licenses/by/4.0/)
This is an open access article distributed in accordance with the Creative Commons Attribution 4.0 Unported (CC BY 4.0) license, which permits others to copy, redistribute, remix, transform and build upon this work for any purpose, provided the original work is properly cited, a link to the licence is given, and indication of whether changes were made. See: [https://creativecommons.org/licenses/by/4.0/](https://creativecommons.org/licenses/by/4.0/).

## References

1.   Cuijpers P , Karyotaki E . Psychological treatment of depression: an updated meta-analytic database of randomized trials. 2022. [doi:10.17605/OSF.IO/825C6](http://dx.doi.org/10.17605/OSF.IO/825C6) Available: [https://doi.org/10.17605/OSF.IO/825C6](https://doi.org/10.17605/OSF.IO/825C6) 
    
    

2.   Cuijpers P , Karyotaki E , de Wit L , et al . The effects of fifteen evidence-supported therapies for adult depression: a meta-analytic review. Psychother Res 2020;30:279–93. [doi:10.1080/10503307.2019.1649732](http://dx.doi.org/10.1080/10503307.2019.1649732) 
    
    [PubMed](http://mentalhealth.bmj.com//lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Febmental%2F26%2F1%2Fe300626.atom) 

3.   Cuijpers P . Four decades of outcome research on psychotherapies for adult depression: an overview of a series of meta-analyses. Can Psychol 2017;58:7–19. [doi:10.1037/cap0000096](http://dx.doi.org/10.1037/cap0000096) 
    
    

4.   Ekers D , Webster L , Van Straten A , et al . Behavioural activation for depression; an update of meta-analysis of effectiveness and sub group analysis. PLoS One 2014;9:e100100. [doi:10.1371/journal.pone.0100100](http://dx.doi.org/10.1371/journal.pone.0100100) 
    
    

5.   Leichsenring F , Rabung S . Effectiveness of long-term psychodynamic psychotherapy: a meta-analysis. JAMA 2008;300:1551–65. [doi:10.1001/jama.300.13.1551](http://dx.doi.org/10.1001/jama.300.13.1551) 
    
    [CrossRef](http://mentalhealth.bmj.com//lookup/external-ref?access_num=10.1001/jama.300.13.1551&link_type=DOI) 
    
    [PubMed](http://mentalhealth.bmj.com//lookup/external-ref?access_num=18827212&link_type=MED&atom=%2Febmental%2F26%2F1%2Fe300626.atom) 
    
    [Web of Science](http://mentalhealth.bmj.com//lookup/external-ref?access_num=000259617700019&link_type=ISI) 

6.   Cuijpers P , Cristea IA , Karyotaki E , et al . How effective are cognitive behavior therapies for major depression and anxiety disorders? A meta-analytic update of the evidence. World Psychiatry 2016;15:245–58. [doi:10.1002/wps.20346](http://dx.doi.org/10.1002/wps.20346) 
    
    [CrossRef](http://mentalhealth.bmj.com//lookup/external-ref?access_num=10.1002/wps.20346&link_type=DOI) 
    
    [PubMed](http://mentalhealth.bmj.com//lookup/external-ref?access_num=27717254&link_type=MED&atom=%2Febmental%2F26%2F1%2Fe300626.atom) 

7.   Cuijpers P , Donker T , Weissman MM , et al . Interpersonal psychotherapy for mental health problems: a comprehensive meta-analysis. Am J Psychiatry 2016;173:680–7. [doi:10.1176/appi.ajp.2015.15091141](http://dx.doi.org/10.1176/appi.ajp.2015.15091141) 
    
    [CrossRef](http://mentalhealth.bmj.com//lookup/external-ref?access_num=10.1176/appi.ajp.2015.15091141&link_type=DOI) 
    
    [PubMed](http://mentalhealth.bmj.com//lookup/external-ref?access_num=27032627&link_type=MED&atom=%2Febmental%2F26%2F1%2Fe300626.atom) 

8.   Cuijpers P , de Wit L , Kleiboer A , et al . Problem-solving therapy for adult depression: an updated meta-analysis. Eur Psychiatry 2018;48:27–37. [doi:10.1016/j.eurpsy.2017.11.006](http://dx.doi.org/10.1016/j.eurpsy.2017.11.006) 
    
    [CrossRef](http://mentalhealth.bmj.com//lookup/external-ref?access_num=10.1016/j.eurpsy.2017.11.006&link_type=DOI) 
    
    [PubMed](http://mentalhealth.bmj.com//lookup/external-ref?access_num=29331596&link_type=MED&atom=%2Febmental%2F26%2F1%2Fe300626.atom) 

9.   Cuijpers P , Driessen E , Hollon SD , et al . The efficacy of non-directive supportive therapy for adult depression: a meta-analysis. Clin Psychol Rev 2012;32:280–91. [doi:10.1016/j.cpr.2012.01.003](http://dx.doi.org/10.1016/j.cpr.2012.01.003) 
    
    [CrossRef](http://mentalhealth.bmj.com//lookup/external-ref?access_num=10.1016/j.cpr.2012.01.003&link_type=DOI) 
    
    [PubMed](http://mentalhealth.bmj.com//lookup/external-ref?access_num=22466509&link_type=MED&atom=%2Febmental%2F26%2F1%2Fe300626.atom) 

10.  Driessen E , Hegelmaier LM , Abbass AA , et al . The efficacy of short-term psychodynamic psychotherapy for depression: a meta-analysis update. Clin Psychol Rev 2015;42:1–15. [doi:10.1016/j.cpr.2015.07.004](http://dx.doi.org/10.1016/j.cpr.2015.07.004) 
    
    [CrossRef](http://mentalhealth.bmj.com//lookup/external-ref?access_num=10.1016/j.cpr.2015.07.004&link_type=DOI) 
    
    [PubMed](http://mentalhealth.bmj.com//lookup/external-ref?access_num=26281018&link_type=MED&atom=%2Febmental%2F26%2F1%2Fe300626.atom) 

11.  Churchill R , Moore THM , Furukawa TA , et al . “Third wave” cognitive and behavioural therapies versus treatment as usual for depression. Cochrane Database Syst Rev 2013:CD008705. [doi:10.1002/14651858.CD008705.pub2](http://dx.doi.org/10.1002/14651858.CD008705.pub2) 
    
    

12.  Leichsenring F , Steinert C , Rabung S , et al . The efficacy of psychotherapies and pharmacotherapies for mental disorders in adults: an umbrella review and meta-analytic evaluation of recent meta-analyses. World Psychiatry 2022;21:133–45. [doi:10.1002/wps.20941](http://dx.doi.org/10.1002/wps.20941) 
    
    

13.  Cuijpers P , van Straten A , Bohlmeijer E , et al . The effects of psychotherapy for adult depression are overestimated: a meta-analysis of study quality and effect size. Psychol Med 2010;40:211–23. [doi:10.1017/S0033291709006114](http://dx.doi.org/10.1017/S0033291709006114) 
    
    [CrossRef](http://mentalhealth.bmj.com//lookup/external-ref?access_num=10.1017/S0033291709006114&link_type=DOI) 
    
    [PubMed](http://mentalhealth.bmj.com//lookup/external-ref?access_num=19490745&link_type=MED&atom=%2Febmental%2F26%2F1%2Fe300626.atom) 

14.  Michopoulos I , Furukawa TA , Noma H , et al . Different control conditions can produce different effect estimates in psychotherapy trials for depression. J Clin Epidemiol 2021;132:59–70. [doi:10.1016/j.jclinepi.2020.12.012](http://dx.doi.org/10.1016/j.jclinepi.2020.12.012) 
    
    [PubMed](http://mentalhealth.bmj.com//lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Febmental%2F26%2F1%2Fe300626.atom) 

15.  Miguel C , Karyotaki E , Cuijpers P , et al . Selective outcome reporting and the effectiveness of psychotherapies for depression. World Psychiatry 2021;20:444–5. [doi:10.1002/wps.20900](http://dx.doi.org/10.1002/wps.20900) 
    
    

16.  Lakens D , Hilgard J , Staaks J . On the reproducibility of meta-analyses: six practical recommendations. BMC Psychol 2016;4:24.:24. [doi:10.1186/s40359-016-0126-3](http://dx.doi.org/10.1186/s40359-016-0126-3) 
    
    [CrossRef](http://mentalhealth.bmj.com//lookup/external-ref?access_num=10.1186/s40359-016-0126-3&link_type=DOI) 
    
    [PubMed](http://mentalhealth.bmj.com//lookup/external-ref?access_num=27241618&link_type=MED&atom=%2Febmental%2F26%2F1%2Fe300626.atom) 

17.  López-Nicolás R , López-López JA , Rubio-Aparicio M , et al . A meta-review of transparency and reproducibility-related reporting practices in published meta-analyses on clinical psychological interventions (2000-2020). Behav Res Methods 2022;54:334–49. [doi:10.3758/s13428-021-01644-z](http://dx.doi.org/10.3758/s13428-021-01644-z) 
    
    [CrossRef](http://mentalhealth.bmj.com//lookup/external-ref?access_num=10.3758/s13428-021-01644-z&link_type=DOI) 
    
    [PubMed](http://mentalhealth.bmj.com//lookup/external-ref?access_num=34173943&link_type=MED&atom=%2Febmental%2F26%2F1%2Fe300626.atom) 

18.  Cuijpers P , Karyotaki E , Reijnders M , et al . Was eysenck right after all? A reassessment of the effects of psychotherapy for adult depression [10]. Epidemiol Psychiatr Sci 2019;28:21–30. [doi:10.1017/S2045796018000057](http://dx.doi.org/10.1017/S2045796018000057) 
    
    [PubMed](http://mentalhealth.bmj.com//lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Febmental%2F26%2F1%2Fe300626.atom) 

19.  Cristea IA , Stefan S , Karyotaki E , et al . The effects of cognitive behavioral therapy are not systematically falling: a revision of johnsen and friborg (2015). Psychol Bull 2017;143:326–40. [doi:10.1037/bul0000062](http://dx.doi.org/10.1037/bul0000062) 
    
    [PubMed](http://mentalhealth.bmj.com//lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Febmental%2F26%2F1%2Fe300626.atom) 

20.  Johnsen TJ , Friborg O . The effects of cognitive behavioral therapy as an anti-depressive treatment is falling: a meta-analysis. Psychol Bull 2015;141:747–68. [doi:10.1037/bul0000015](http://dx.doi.org/10.1037/bul0000015) 
    
    [CrossRef](http://mentalhealth.bmj.com//lookup/external-ref?access_num=10.1037/bul0000015&link_type=DOI) 
    
    [PubMed](http://mentalhealth.bmj.com//lookup/external-ref?access_num=25961373&link_type=MED&atom=%2Febmental%2F26%2F1%2Fe300626.atom) 

21.  Voracek M , Kossmeier M , Tran US . Which data to meta-analyze, and how?: A specification-curve and multiverse-analysis approach to meta-analysis. Z Psychol 2019;227:64–82. [doi:10.1027/2151-2604/a000357](http://dx.doi.org/10.1027/2151-2604/a000357) 
    
    

22.  Cuijpers P , van Straten A , Warmerdam L , et al . Psychological treatment of depression: a meta-analytic database of randomized studies. BMC Psychiatry 2008;8:36. [doi:10.1186/1471-244X-8-36](http://dx.doi.org/10.1186/1471-244X-8-36) 
    
    

23.  Plessen CY , Karyotaki E , Cuijpers P . Exploring the efficacy of psychological treatments for depression: a multiverse meta-analysis protocol. BMJ Open 2022;12:e050197. [doi:10.1136/bmjopen-2021-050197](http://dx.doi.org/10.1136/bmjopen-2021-050197) 
    
    

24.  Plessen C . Data from: exploring the efficacy of psychological treatments for depression: A multiverse meta-analysis. Open Science Framework November 8, 2022. [doi:10.17605/OSF.IO/MTX8A](http://dx.doi.org/10.17605/OSF.IO/MTX8A) Available: [https://doi.org/10.17605/OSF.IO/MTX8A](https://doi.org/10.17605/OSF.IO/MTX8A) 
    
    

25.  Team RC . R: A language and environment for statistical computing. R foundation for statistical computing, vienna. 2020. Available: [www.R-project.org](http://www.R-project.org) 
    
    

26.  Viechtbauer W . Conducting meta-analyses in R with the metafor package. J Stat Soft 2010;36. [doi:10.18637/jss.v036.i03](http://dx.doi.org/10.18637/jss.v036.i03) 
    
    

27.  Harrer M , Kuper P , Sprenger A , et al . MetapsyTools: several R helper functions for the ‘metapsy’ database. 2022. Available: [https://tools.metapsy.org](https://tools.metapsy.org) 
    
    

28.  Higgins JPT , Altman DG , Gøtzsche PC . The cochrane collaboration’s tool for assessing risk of bias in randomised trials. BMJ 2011;343:d5928. [doi:10.1136/bmj.d5928](http://dx.doi.org/10.1136/bmj.d5928) 
    
    

29.  Pustejovsky JE , Tipton E . Meta-Analysis with robust variance estimation: expanding the range of working models. Prev Sci 2022;23:425–38. [doi:10.1007/s11121-021-01246-3](http://dx.doi.org/10.1007/s11121-021-01246-3) 
    
    

30.  Bartoš F , Maier M , Quintana DS , et al . Adjusting for publication bias in JASP and R: selection models, PET-PEESE, and robust bayesian meta-analysis. Advances in Methods and Practices in Psychological Science 2022;5:251524592211092. [doi:10.1177/25152459221109259](http://dx.doi.org/10.1177/25152459221109259) 
    
    

31.  van der Meulen E , van der Velden PG , van Aert RCM , et al . Longitudinal associations of psychological resilience with mental health and functioning among military personnel: a meta-analysis of prospective studies. Soc Sci Med 2020;255:112814. [doi:10.1016/j.socscimed.2020.112814](http://dx.doi.org/10.1016/j.socscimed.2020.112814) 
    
    

32.  Borenstein M , Hedges LV , Higgins JPT , et al . Introduction to meta-analysis. Chichester, U.K: John Wiley & Sons, 2011.
    
    

33.  Higgins JPT , Thomas J , Chandler J , et al . Cochrane handbook for systematic reviews of interventions. John Wiley & Sons, 23 September 2019. [doi:10.1002/9781119536604](http://dx.doi.org/10.1002/9781119536604) 
    
    

34.  Stanley TD , Doucouliagos H . Meta-Regression approximations to reduce publication selection bias. Res Synth Methods 2014;5:60–78. [doi:10.1002/jrsm.1095](http://dx.doi.org/10.1002/jrsm.1095) 
    
    

35.  Egger M , Davey Smith G , Schneider M , et al . Bias in meta-analysis detected by a simple, graphical test. BMJ 1997;315:629–34. [doi:10.1136/bmj.315.7109.629](http://dx.doi.org/10.1136/bmj.315.7109.629) 
    
    [Abstract/FREE Full Text](http://mentalhealth.bmj.com//lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjEyOiIzMTUvNzEwOS82MjkiO3M6NDoiYXRvbSI7czoyNzoiL2VibWVudGFsLzI2LzEvZTMwMDYyNi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

36.  Duval S , Tweedie R . Trim and fill: a simple funnel-plot-based method of testing and adjusting for publication bias in meta-analysis. Biometrics 2000;56:455–63. [doi:10.1111/j.0006-341x.2000.00455.x](http://dx.doi.org/10.1111/j.0006-341x.2000.00455.x) 
    
    [CrossRef](http://mentalhealth.bmj.com//lookup/external-ref?access_num=10.1111/j.0006-341X.2000.00455.x&link_type=DOI) 
    
    [PubMed](http://mentalhealth.bmj.com//lookup/external-ref?access_num=10877304&link_type=MED&atom=%2Febmental%2F26%2F1%2Fe300626.atom) 
    
    [Web of Science](http://mentalhealth.bmj.com//lookup/external-ref?access_num=000087677500019&link_type=ISI) 

37.  Cuijpers P , Turner EH , Koole SL , et al . What is the threshold for a clinically relevant effect? the case of major depressive disorders. Depress Anxiety 2014;31:374–8. [doi:10.1002/da.22249](http://dx.doi.org/10.1002/da.22249) 
    
    [CrossRef](http://mentalhealth.bmj.com//lookup/external-ref?access_num=10.1002/da.22249&link_type=DOI) 
    
    [PubMed](http://mentalhealth.bmj.com//lookup/external-ref?access_num=24677535&link_type=MED&atom=%2Febmental%2F26%2F1%2Fe300626.atom) 
    
    [Web of Science](http://mentalhealth.bmj.com//lookup/external-ref?access_num=000334803100002&link_type=ISI) 

38.  Olkin I , Dahabreh IJ , Trikalinos TA . GOSH - a graphical display of study heterogeneity: all-subsets meta-analysis for heterogeneity. Res Synth Methods 2012;3:214–23. [doi:10.1002/jrsm.1053](http://dx.doi.org/10.1002/jrsm.1053) 
    
    

39.  Furukawa TA , Noma H , Caldwell DM , et al . Waiting list may be a nocebo condition in psychotherapy trials: a contribution from network meta-analysis. Acta Psychiatr Scand 2014;130:181–92. [doi:10.1111/acps.12275](http://dx.doi.org/10.1111/acps.12275) 
    
    [CrossRef](http://mentalhealth.bmj.com//lookup/external-ref?access_num=10.1111/acps.12275&link_type=DOI) 
    
    [PubMed](http://mentalhealth.bmj.com//lookup/external-ref?access_num=24697518&link_type=MED&atom=%2Febmental%2F26%2F1%2Fe300626.atom) 

40.  Harrer M , Cuijpers P , Furukawa TA , et al . Doing meta-analysis with R. 1st ed. Boca Raton: Chapman & Hall/CRC Press, 2021. [doi:10.1201/9781003107347](http://dx.doi.org/10.1201/9781003107347) 
    
    

41.  Driessen E , Hollon SD , Bockting CLH , et al . Does publication bias inflate the apparent efficacy of psychological treatment for major depressive disorder? A systematic review and meta-analysis of US National Institutes of health-funded trials. PLOS ONE 2015;10:e0137864. [doi:10.1371/journal.pone.0137864](http://dx.doi.org/10.1371/journal.pone.0137864) 
    
    

42.  Moshe I , Terhorst Y , Philippi P , et al . Digital interventions for the treatment of depression: a meta-analytic review. Psychol Bull 2021;147:749–86. [doi:10.1037/bul0000334](http://dx.doi.org/10.1037/bul0000334) 
    
    [PubMed](http://mentalhealth.bmj.com//lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Febmental%2F26%2F1%2Fe300626.atom) 

43.  Cuijpers P , Karyotaki E , Eckshtain D , et al . Psychotherapy for depression across different age groups: a systematic review and meta-analysis. JAMA Psychiatry 2020;77:694–702. [doi:10.1001/jamapsychiatry.2020.0164](http://dx.doi.org/10.1001/jamapsychiatry.2020.0164) 
    
    

44.  Cuijpers P , Noma H , Karyotaki E , et al . Effectiveness and acceptability of cognitive behavior therapy delivery formats in adults with depression: a network meta-analysis. JAMA Psychiatry 2019;76:700–7. [doi:10.1001/jamapsychiatry.2019.0268](http://dx.doi.org/10.1001/jamapsychiatry.2019.0268) 
    
    

45.  Cuijpers P , Karyotaki E , Reijnders M , et al . Is psychotherapy effective? pretending everything is fine will not help the field forward. Epidemiol Psychiatr Sci 2019;28:356–7. [doi:10.1017/S204579601800080X](http://dx.doi.org/10.1017/S204579601800080X)

 [1]: /embed/mml-math-1.gif
 [2]: /embed/mml-math-2.gif