Self-Guided Interventions for Depression
A living database of randomized controlled trials — filter the evidence and run a full meta-analysis online.
Additional Information
This online meta-analysis tool allows you to explore a simplified version of the “Self-Guided Interventions for Depression” dataset. It is part of the meta-analytic research domain (MARD) on psychotherapy for depression. The code used to simplify the data for this application can be downloaded here.
Select outcome
Choose which outcome domain to analyze. Only studies reporting that outcome will appear on the next step.
In this step you choose which of two outcome variables to analyze: depression or quality of life.
Not all studies in the database report both outcome variables. Once you select an outcome, only studies reporting that outcome appear in the data table. All filters on the next step apply only to those studies. You can return here at any time to select a different outcome.
Select your data
Apply filters to choose trials for your meta-analysis. All trials are selected by default — refine on the right, the table updates automatically, then Run meta-analysis.
Meta-analysis results
Results for the trials you selected. Use the tabs below for pooled effects, plots, outliers, publication bias, risk of bias and moderators.
This tab shows you the results of the meta-analysis conducted for your selected data. In the window below, you can switch from one panel to another to see various results, including the main effect size and heterogeneity, the forest plot, results of outlier analyses, publication bias analyses and moderator analyses.
By default, the web application automatically makes a few decisions concerning the analysis settings. These can be seen below under "Analysis Settings". By default, the application uses a Random-Effects Pooling Model for your meta-analysis. This is nearly always the correct choice when analyzing data in mental health research. In addition, the Paule-Mandel estimator is used by default to calculate the between-study heterogeneity. However, when clicking on the dropdown menu "Pooling Model", you can see that there are various estimators available when assuming a Random-Effects (RE) model. Sometimes, the choice of the between-study heterogeneity estimator can have an impact on the pooled results, and not every estimator is optimal under all circumstances. If you want to learn more about this, you can consult Harrer, Cuijpers, Furukawa & Ebert, 2020. It is also possible to calculate effects using a Fixed-Effect Model; but this model should only be used if you have good reasons for applying it.
Under "Analyzed Moderators", the moderating variables which are currently inspected in meta-regression/subgroup analyses in the "Moderator Analysis" tab are shown. By default, the publication year is used for a meta-regression on your selected data, and the country/region of a study is used as a subgroup analysis. You can easily add more moderator analyses by clicking on the white box. A dropdown will then appear, showing more variables which are available. You can then add those variables by clicking on them. It is also possible to remove variables: again click on the white box and use the backspace ← key.
By default, the application also uses Knapp-Hartung adjustments to calculate the confidence interval for your meta-analysis result. It is possible to disable this method by unchecking the "Use Knapp-Hartung Adjustments" checkbox. You can learn more about this method in Harrer, Cuijpers, Furukawa & Ebert, 2020.
Finally, after you have completed the reconfiguration of your meta-analysis settings, simply click the "Re-run Meta-Analysis" button. This will recalculate the results of your meta-analysis using the new settings you have provided. Please note that the reanalysis may take some time before it is finished.
A random-effects meta-analysis does not assume that every study shares one identical true effect. Instead, the true effects are assumed to be normally distributed. This chart shows that assumed distribution: a bell curve centred on the pooled estimate, whose width is the between-study standard deviation (τ). A wider curve means more heterogeneity between studies.
The thick vertical line marks the pooled overall effect (Hedges’ g). The shaded area under the curve is its 95% confidence interval — the range of plausible values for the average true effect.
The two short solid lines in the tails mark the 95% prediction interval: the range in which the true effect of a future, comparable study would be expected to fall.
The vertical line at 0 is the point of no effect. If the effect and its intervals lie clearly to one side of it, the intervention has a consistent direction of effect.
Individual studies with extreme effects can distort a pooled result and inflate heterogeneity. This is a sensitivity analysis: studies flagged as statistical outliers — those whose 95% confidence interval does not overlap the confidence interval of the pooled effect — are removed, and the meta-analysis is re-run on the remaining studies.
The numbers show the pooled effect and heterogeneity after removing outliers, and the list shows which studies were removed (with their effect sizes).
The chart compares the true-effect distribution with outliers removed (coloured) against the original analysis (grey). If the two are similar, the result is robust to outliers; if they differ markedly, the pooled effect is sensitive to a few extreme studies and should be interpreted with caution.
The Baujat plot helps identify studies that have a disproportionate impact on the meta-analysis.
The horizontal axis shows each study’s contribution to the overall heterogeneity; the vertical axis shows its influence on the pooled effect (how much the result changes when the study is removed).
Studies in the top-right corner are both major sources of heterogeneity and highly influential — these are the most important to inspect when judging how robust your findings are.
A funnel plot plots each study's effect against its precision (standard error). Without bias, studies scatter symmetrically around the pooled effect, forming an inverted funnel — smaller, less precise studies scatter more widely toward the bottom.
Asymmetry (a gap in one bottom corner) can signal publication bias — e.g. small studies with non-significant results going unpublished — though it can also reflect genuine heterogeneity or chance.
The shaded contours mark conventional significance regions (p < .05, .025, .01), which help judge whether any missing studies would fall in non-significant areas.
Egger's regression test quantifies the asymmetry by testing whether the regression intercept differs from zero. A significant intercept (p < .05) indicates asymmetry. Interpret with caution when few studies are included.
A p-curve examines the distribution of the statistically significant p-values (p < .05) across your studies to judge whether they hold evidential value — a genuine effect — or mainly reflect selective reporting / p-hacking.
If a true effect exists, significant p-values cluster near .01 (the curve is right-skewed). With no effect, they are uniform (flat). If results were pushed just past significance, the curve is left-skewed (piling up near .05).
Your observed curve is compared against two references: the null of no effect (flat, 20% per bin) and the null of 33% power.
Evidential value present means the observed curve is significantly right-skewed. Evidential value inadequate means it is flatter than expected even from studies with 33% power. The power estimate approximates the average statistical power of the significant studies.
Subgroup Analysis
About this web app
This application was developed by researchers at the Vrije Universiteit Amsterdam, The Netherlands, the Friedrich-Alexander-Universität Erlangen-Nürnberg, Germany, and the Technical University Munich, Germany. The data are based on a meta-analytic database that was developed by researchers at the Vrije Universiteit Amsterdam in the past 14 years.
This project, led by Prof. Pim Cuijpers and Dr. Eirini Karyotaki, has resulted in a long series of published studies in peer-reviewed journals. Researchers from the Friedrich-Alexander-Universität in Germany, led by Mathias Harrer and Dr. David Ebert (Technical University Munich), have developed the web application and the automated analyses that are conducted online.
About the Metapsy Web Application
This web application serves as an interface to the Metapsy meta-analytic database. The application uses a Shiny server instance to access the database, enable data downloads, and handle computations using state-of-the-art meta-analytic techniques.
The meta (Balduzzi, Rücker & Schwarzer, 2019) and dmetar (Harrer, Cuijpers, Furukawa & Ebert, 2020) R packages are used to perform the meta-analyses.
The metagen function is used internally to perform the meta-analytic pooling. By default, the Restricted Maximum-Likelihood estimator (Viechtbauer, 2005) is used to estimate the between-study heterogeneity/variance τ2 in a random-effects model; REML is also used for meta-regressions. To calculate the Number Needed to Treat (NNT), the method by Furukawa and Leucht (Furukawa & Leucht, 2011) is used, assuming Control Event Rate (CER) of 0.19, which is derived from Cuijpers et al., 2014. For subgroup analyses, the subgroup.analysis.mixed.effects function is used, implementing a mixed-effect model for which results within subgroups are pooled using a random-effects model (inheriting the τ2-estimator specified for the overall analysis), and results between subgroups are compared assuming a fixed-effects model.
To generate forest plots, the forest function in meta is used.
Outlier selection is conducted using the find.outliers function in dmetar. The function implements a simple outlier removal algorithm where all results for which the 95% confidence interval is outside the 95% confidence interval of the pooled effect are removed as outliers. Due to high computational costs, sensitivity analyses based on the “Leave-One-Out” paradigm are not conducted for large meta-analyses (k>50). For smaller meta-analyses, a Baujat plot (Baujat et al., 2002) is created using the baujat function in meta.
Publication bias analyses are conducted using the eggers.test function in dmetar (which is a wrapper for the metabias function in meta). A P-curve is created using the pcurve function in dmetar. Please note that several prerequisites should be considered before P-curves can be interpreted (Simonsohn et al., 2014).
For the risk of bias overview, the rob.summary function in dmetar is used.
The metapsyData R package
Access the Metapsy meta-analytic databases directly in R.
The metapsyData package lets you access the Metapsy meta-analytic psychotherapy databases — including the “Self-Guided Interventions for Depression” database powering this app — directly in your R environment.
Installation
if (!require("devtools"))
install.packages("devtools")
devtools::install_github("metapsy-project/metapsyData")
Usage
Once installed, run getData() to load a database locally. Databases are referenced by their identifier; this app uses "depression-selfguided-psyctr":
library(metapsyData)
library(dplyr)
d <- getData("depression-selfguided-psyctr")
glimpse(d$data)
#> Rows: 361
#> Columns: 90
#> $ study <chr> "Addington, 2019", "Addington, 2019", "Addington, …
#> $ condition_arm1 <chr> "other psy", "other psy", "other psy", "other psy"…
#> $ condition_arm2 <chr> "wl", "other ctr", "wl", "other ctr", "wl", "other…
#> $ descr_arm1 <chr> "positive emtion skills", "positive emtion skills"…
#> $ descr_arm2 <chr> "wl", "daily emotional report", "wl", "daily emoti…
#> $ multi_arm1 <chr> "other psy", "other psy", "other psy", "other psy"…
#> $ multi_arm2 <chr> "wl", "other ctr", "wl", "other ctr", "wl", "other…
#> $ .g <dbl> 0.22279049, 0.27338834, 0.32063561, 0.21349628, 0.…
#> $ .g_se <dbl> 0.3252040, 0.3133481, 0.3262391, 0.3127861, 0.3257…
#> $ primary_calc <dbl> 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,…
#> $ outcome_type <chr> "msd", "msd", "msd", "msd", "msd", "msd", "msd", "…
#> $ outcome_domain <chr> "depression", "depression", "depression", "depress…
#> $ time <chr> "post", "post", "post", "post", "FU1", "FU1", "FU1…
#> $ time_weeks <dbl> 7, 7, 7, 7, 11, 11, 11, 11, 19, 19, 19, 19, 4, 6, …
#> $ instrument <chr> "PHQ-8", "PHQ-8", "CES-D", "CES-D", "PHQ-8", "PHQ-…
#> $ rating <chr> "self-report", "self-report", "self-report", "self…
#> $ mean_arm1 <dbl> 14.78, 14.78, 24.63, 24.63, 14.75, 14.75, 25.01, 2…
#> $ mean_arm2 <dbl> 15.92, 16.18, 27.93, 26.83, 16.17, 15.70, 29.51, 2…
#> $ sd_arm1 <dbl> 5.05, 5.05, 10.15, 10.15, 5.05, 5.05, 10.20, 10.20…
#> $ sd_arm2 <dbl> 4.96, 4.99, 9.99, 10.06, 4.96, 5.15, 9.99, 10.35, …
#> $ n_arm1 <dbl> 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, NA…
#> $ n_arm2 <dbl> 15, 17, 15, 17, 15, 17, 15, 17, 15, 17, 15, 17, NA…
#> $ bl_mean_arm1 <dbl> 1546, 1546, 2403, 2403, 1546, 1546, 2403, 2403, 15…
#> $ bl_sd_arm1 <dbl> 739, 739, 1112, 1112, 739, 739, 1112, 1112, 739, 7…
#> $ bl_n_arm1 <dbl> 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 21…
#> $ bl_mean_arm2 <dbl> 1546, 1705, 2542, 2641, 1546, 1705, 2542, 2641, 15…
#> $ bl_sd_arm2 <dbl> 562, 565, 1099, 1100, 562, 565, 1099, 1100, 562, 5…
#> $ bl_n_arm2 <dbl> 15, 17, 15, 17, 15, 17, 15, 17, 15, 17, 15, 17, 20…
#> $ dich <chr> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, "r…
#> $ dich_type <chr> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, "r…
#> $ event_arm1 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 36…
#> $ event_arm2 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 34…
#> $ totaln_arm1 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 79…
#> $ totaln_arm2 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 10…
#> $ change <chr> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA…
#> $ mean_change_arm1 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA…
#> $ mean_change_arm2 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA…
#> $ sd_change_arm1 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA…
#> $ sd_change_arm2 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA…
#> $ n_change_arm1 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA…
#> $ n_change_arm2 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA…
#> $ other_statistic <chr> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA…
#> $ precalc_g <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA…
#> $ precalc_g_se <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA…
#> $ rand_arm1 <dbl> 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 26, 21…
#> $ rand_arm2 <dbl> 15, 17, 15, 17, 15, 17, 15, 17, 15, 17, 15, 17, 20…
#> $ n_drop_study1 <dbl> 0, 0, 0, 0, 4, 4, 4, 4, 4, 4, 4, 4, 132, 5, 9, 62,…
#> $ n_drop_study2 <dbl> 1, 1, 1, 1, 1, 2, 1, 2, 0, 2, 0, 2, 106, 0, 6, 10,…
#> $ n_drop_treat1 <dbl> NA, NA, NA, NA, 2, NA, NA, NA, NA, NA, NA, NA, NA,…
#> $ n_drop_treat2 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA…
#> $ def_drop_treat <chr> NA, NA, NA, NA, "did not complete at least 4 out o…
#> $ year <dbl> 2019, 2019, 2019, 2019, 2019, 2019, 2019, 2019, 20…
#> $ comorbid_mental <chr> "n", "n", "n", "n", "n", "n", "n", "n", "n", "n", …
#> $ comorbid_physical <chr> "n", "n", "n", "n", "n", "n", "n", "n", "n", "n", …
#> $ format_unguided <chr> "internet", "internet", "internet", "internet", "i…
#> $ commericial <chr> "NI", "NI", "NI", "NI", "NI", "NI", "NI", "NI", "N…
#> $ program_name <chr> "MARIGOLD", "MARIGOLD", "MARIGOLD", "MARIGOLD", "M…
#> $ support_unguided <chr> "pure", "pure", "pure", "pure", "pure", "pure", "p…
#> $ initial_screening <chr> "y", "y", "y", "y", "y", "y", "y", "y", "y", "y", …
#> $ forum_unguided <chr> "n", "n", "n", "n", "n", "n", "n", "n", "n", "n", …
#> $ n_sessions_planned <dbl> 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, NA, 6, 6, 10, …
#> $ n_sessions_received <dbl> 7.5, 7.5, 7.5, 7.5, NA, NA, NA, NA, NA, NA, NA, NA…
#> $ country <chr> "us", "us", "us", "us", "us", "us", "us", "us", "u…
#> $ income_level <chr> "high", "high", "high", "high", "high", "high", "h…
#> $ age_group <chr> "adul", "adul", "adul", "adul", "adul", "adul", "a…
#> $ mean_age <dbl> 42.16, 42.16, 42.16, 42.16, 42.16, 42.16, 42.16, 4…
#> $ percent_women <dbl> 0.655, 0.655, 0.655, 0.655, 0.655, 0.655, 0.655, 0…
#> $ recruitment <chr> "com", "com", "com", "com", "com", "com", "com", "…
#> $ diagnosis <chr> "cut", "cut", "cut", "cut", "cut", "cut", "cut", "…
#> $ target_group <chr> "adul", "adul", "adul", "adul", "adul", "adul", "a…
#> $ sg <chr> "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", …
#> $ ac <chr> "?", "?", "?", "?", "?", "?", "?", "?", "?", "?", …
#> $ ba <chr> "sr", "sr", "sr", "sr", "sr", "sr", "sr", "sr", "s…
#> $ itt <dbl> 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,…
#> $ rob <dbl> 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 4, 4, 3, 2,…
#> $ domain.1 <chr> "Some concerns", "Some concerns", "Some concerns",…
#> $ domain.2 <chr> "Low risk", "Low risk", "Low risk", "Low risk", "L…
#> $ domain.3 <chr> "Low risk", "Low risk", "Low risk", "Low risk", "L…
#> $ domain.4 <chr> "Low risk", "Low risk", "Low risk", "Low risk", "L…
#> $ domain.5 <chr> "Some concerns", "Some concerns", "Some concerns",…
#> $ overall <chr> "Some concerns", "Some concerns", "Some concerns",…
#> $ .id <chr> "Addington_2019_msd_PHQ-8_post_7_self-report_other…
#> $ .g_negative <dbl> -0.22279049, -0.27338834, -0.32063561, -0.21349628…
#> $ .log_rr <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 0.…
#> $ .log_rr_se <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 0.…
#> $ .event_arm1 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 36…
#> $ .event_arm2 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 34…
#> $ .totaln_arm1 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 79…
#> $ .totaln_arm2 <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 10…
#> $ full_reference <chr> "Addington, E. L., Cheung, E. O., Bassett, S. M., …
The raw data files can also be browsed in the GitHub repository ↗ under data, and the full package documentation is hosted on rdrr.io ↗.
Citation
Run citation("metapsyData") to retrieve the current citation:
Harrer, M., Ebert, D. D., Karyotaki, E., & Cuijpers, P. (2022).
metapsyData: Access the Meta-Analytic Psychotherapy Database.
R package version 0.1.0. DOI: 10.5281/zenodo.5171880.