---
title: "Zplot Publication-Bias Diagnostics"
author: "František Bartoš"
date: "8th of September 2025 (updated: 29th of April 2026)"
output:
  rmarkdown::html_vignette:
    self_contained: yes
bibliography: ../inst/REFERENCES.bib
csl: ../inst/apa.csl
link-citations: true
vignette: >
  %\VignetteIndexEntry{Zplot Publication-Bias Diagnostics}
  %\VignetteEngine{knitr::rmarkdown}
  %\VignetteEncoding{UTF-8}
  %\VignetteEngine{knitr::rmarkdown_notangle}
---

```{r child = "_vignette-nowrap.md", echo = FALSE, eval = TRUE}
```

```{r setup, include = FALSE}
source("vignette-cache.R", local = knitr::knit_global())

cached_fits <- vignette_cache(
  name     = "v36-zplot",
  objects  = c(
    "fit_RE_Weingarten2018",
    "fit_RoBMA_Weingarten2018",
    "fit_RE_Hoppen2025",
    "fit_RoBMA_Hoppen2025",
    "fit_RE_Wang2025",
    "fit_RoBMA_Wang2025",
    "fit_RE_ManyLabs16",
    "fit_RoBMA_ManyLabs16"
  ),
  packages = c("metafor")
)

knitr::opts_chunk$set(
  collapse   = TRUE,
  comment    = "#>",
  eval       = vignette_cache_eval(cached_fits),
  message    = FALSE,
  warning    = FALSE,
  fig.width  = 7,
  fig.height = 5,
  dev        = "png",
  fig.retina = 3
)
if (.Platform$OS.type == "windows") {
  knitr::opts_chunk$set(dev.args = list(type = "cairo"))
}
```

```{r load-models, include = FALSE}
library("RoBMA")
vignette_cache_load(cached_fits)
```

```{r precompute-models, include = FALSE, eval = FALSE}
library("RoBMA")

# Ease of retrieval - Weingarten2018
data("Weingarten2018", package = "RoBMA")
Weingarten2018 <- Weingarten2018[Weingarten2018$standard_paradigm & Weingarten2018$proximal_dataset, ]
Weingarten2018_z <- metafor::escalc(
  ri = r_xy, ni = round(N), measure = "ZCOR", data = Weingarten2018
)
fit_RE_Weingarten2018 <- brma(
  yi = yi, vi = vi, measure = "ZCOR",
  cluster = paper_id, data = Weingarten2018_z,
  sample = 10000, burnin = 10000, adapt = 10000,
  chains = 5, parallel = TRUE, seed = 1
)
fit_RoBMA_Weingarten2018 <- RoBMA(
  yi = yi, vi = vi, measure = "ZCOR",
  cluster = paper_id, data = Weingarten2018_z,
  sample = 10000, burnin = 10000, adapt = 10000,
  chains = 5, parallel = TRUE, seed = 1
)

# Social comparison - Hoppen2025
data("Hoppen2025", package = "RoBMA")
fit_RE_Hoppen2025 <- brma(
  yi = d, sei = sqrt(v), measure = "SMD",
  data = Hoppen2025,
  sample = 10000, burnin = 5000, adapt = 5000,
  chains = 5, parallel = TRUE, seed = 1
)
fit_RoBMA_Hoppen2025 <- RoBMA(
  yi = d, sei = sqrt(v), measure = "SMD",
  data = Hoppen2025,
  sample = 10000, burnin = 5000, adapt = 5000,
  chains = 5, parallel = TRUE, seed = 1
)

# ChatGPT - Wang2025
data("Wang2025", package = "RoBMA")
Wang2025 <- Wang2025[Wang2025$Learning_effect == "Learning performance", ]
fit_RE_Wang2025 <- brma(
  yi = g, sei = se, measure = "SMD",
  data = Wang2025,
  sample = 10000, burnin = 5000, adapt = 5000,
  chains = 5, parallel = TRUE, seed = 1
)
fit_RoBMA_Wang2025 <- RoBMA(
  yi = g, sei = se, measure = "SMD",
  data = Wang2025,
  sample = 10000, burnin = 5000, adapt = 5000,
  chains = 5, parallel = TRUE, seed = 1
)

# Many Labs 2 - ManyLabs16
data("ManyLabs16", package = "RoBMA")
fit_RE_ManyLabs16 <- brma(
  yi = y, sei = se, measure = "SMD",
  data = ManyLabs16,
  sample = 10000, burnin = 5000, adapt = 5000,
  chains = 5, parallel = TRUE, seed = 1
)
fit_RoBMA_ManyLabs16 <- RoBMA(
  yi = y, sei = se, measure = "SMD",
  data = ManyLabs16,
  sample = 10000, burnin = 5000, adapt = 5000,
  chains = 5, parallel = TRUE, seed = 1
)

vignette_cache_save(cached_fits)
```

**This vignette accompanies the manuscript "Zplot: A Visual Diagnostic for Publication Bias in Meta-Analysis" [@bartos2025zcurve].**

The manuscript provides the theoretical foundation and detailed methodology for zplot diagnostics, while this vignette demonstrates their practical implementation using the `RoBMA` R package [@RoBMA].

Zplots visually assess meta-analytic model fit, focusing on misfit due to publication bias [@bartos2025zcurve].
The methodology overlays model-implied posterior predictive distributions of z-statistics on the observed distribution of z-statistics [@gabry2019visualization], allowing researchers to assess how well different meta-analytic models capture the patterns in their data.
The approach builds on earlier work by @brunner2020estimating for assessing the quality of research.

The zplot approach complements statistical tests of publication bias (such as inclusion Bayes factors) with intuitive visualizations that can highlight model misfit due to publication bias.
The method also allows extrapolation to the pre-publication bias state, providing estimates of key metrics such as the expected discovery rate and the number of missing studies.

We illustrate the zplot diagnostics using four empirical examples from recent meta-analyses that exhibit different degrees of publication bias.
The examples demonstrate how to interpret zplots and use them to guide model selection in meta-analytic practice.
For details refer to the accompanying manuscript [@bartos2025zcurve].

## Getting Started

Before we start, we load the `RoBMA` R package:

```{r, message = FALSE}
library("RoBMA")
```

The zplot workflow consists of the following steps:

1. Fit meta-analytic models to your data using `brma()` for random-effects models or `RoBMA()` for publication-bias-adjusted models [@maier2020robust; @bartos2021no].
2. Call `zplot()` directly on a fitted model for a one-shot diagnostic plot, or convert via `as_zplot()` when you want to overlay multiple models.
3. Generate histograms of observed z-statistics using `hist()` on a `zplot_brma` object.
4. Overlay model fits using `lines()` to compare different models.
5. Use the visualization to supplement statistical model comparison.

`brma()` fits a single random-effects meta-analytic model, while `RoBMA()` fits a Bayesian model-averaged ensemble that includes publication-bias adjustment.
The zplot diagnostics work with fitted `brma`/`RoBMA` objects.

## Applied Examples

We demonstrate the zplot diagnostics on four empirical meta-analyses that represent different publication-bias scenarios:

1. **Ease-of-retrieval effect in the few/many paradigm** [@weingarten2018does]: an example with extreme publication bias.
2. **Social comparison as behavior change technique** [@hoppen2025meta]: an example with strong publication bias.
3. **ChatGPT effects on learning performance** [@wang2025effect]: an example with moderate publication bias.
4. **Framing effects from Many Labs 2** [@klein2018many]: an example with no publication bias (registered replication reports).

For each example, we fit both a simple random-effects model using `brma()` and a publication-bias-adjusted ensemble using `RoBMA()` [@maier2020robust; @bartos2021no].
We then use zplots to compare how well each model captures the observed distribution of z-statistics.

### Example 1: Ease-of-Retrieval Effect in the Few/Many Paradigm

This example examines the ease-of-retrieval effect, a well-established phenomenon in cognitive psychology where the subjective difficulty of recalling information influences judgments more than the actual number of items recalled [@schwarz1991ease].

We analyze 298 estimates from 111 studies examining the ease-of-retrieval effect in the few/many standard paradigm in the proximal dataset collected by Weingarten and colleagues [@weingarten2018does].
The original analysis reported a pooled effect size r = 0.25, 95% CI [0.22, 0.28].
When adjusted for publication bias using PET-PEESE, the effect was reduced but remained substantial at r = 0.19, 95% CI [0.15, 0.23].

**Data and model fitting.**
We begin by loading the ease-of-retrieval dataset and examining its structure, focusing on studies using the standard paradigm in proximal dataset conditions:

```{r}
data("Weingarten2018", package = "RoBMA")
Weingarten2018 <- Weingarten2018[Weingarten2018$standard_paradigm & Weingarten2018$proximal_dataset, ]
head(Weingarten2018)
```

The filtered dataset contains 298 effect size estimates (correlation coefficients) and sample sizes.
Since the effect size estimates are nested within studies (indicated by the `paper_id` variable), we specify a multilevel random-effects model using `brma()` and a multilevel publication-bias-adjusted model using `RoBMA()` [@bartos2025robust].
We use `metafor::escalc()` to compute Fisher's z effect sizes from the correlation coefficients and sample sizes prior to fitting.

```{r, eval = FALSE}
Weingarten2018_z <- metafor::escalc(
  ri = r_xy, ni = round(N), measure = "ZCOR", data = Weingarten2018
)

fit_RE_Weingarten2018 <- brma(
  yi = yi, vi = vi, measure = "ZCOR",
  cluster = paper_id, data = Weingarten2018_z,
  sample = 10000, burnin = 10000, adapt = 10000,
  chains = 5, parallel = TRUE, seed = 1)

fit_RoBMA_Weingarten2018 <- RoBMA(
  yi = yi, vi = vi, measure = "ZCOR",
  cluster = paper_id, data = Weingarten2018_z,
  sample = 10000, burnin = 10000, adapt = 10000,
  chains = 5, parallel = TRUE, seed = 1)
```

**Model results.**
We examine the results from both models using the `summary()` function.

```{r}
summary(fit_RE_Weingarten2018)

summary(fit_RoBMA_Weingarten2018)
```

The Bayesian multilevel random-effects model finds results similar to those reported in the original publication-bias-unadjusted analysis.

The RoBMA analysis reveals extreme evidence for publication bias.
However, RoBMA still finds extreme evidence for the presence of an effect, with a substantially reduced but non-zero model-averaged effect size estimate.

**Zplot diagnostics.**
We now generate zplots to assess how well each model captures the observed distribution of z-statistics.
The `hist()` function creates a histogram of the observed z-statistics, while `lines()` overlays the model-implied posterior predictive distributions [@gabry2019visualization].
The `as_zplot()` function transforms a fitted `brma`/`RoBMA` object into a `zplot_brma` object that `hist()`, `lines()`, and `plot()` can consume:

```{r, fig.cap="Ease-of-Retrieval Effect: Model Fit Assessment"}
hist(as_zplot(fit_RoBMA_Weingarten2018), from = -3, to = 6, by = 0.25)

lines(as_zplot(fit_RE_Weingarten2018),    from = -3, to = 6, col = "black", lty = 2, lwd = 2)
lines(as_zplot(fit_RoBMA_Weingarten2018), from = -3, to = 6, col = "blue",  lty = 2, lwd = 2)

legend(
  "topright",
  legend = c("Random-Effects", "RoBMA"),
  col = c("black", "blue"),
  lty = 2, lwd = 2)
```

The zplot reveals clear evidence of extreme publication bias in the ease-of-retrieval literature.
Two extreme discontinuities are visible in the observed distribution of z-statistics (gray bars):

1. **Marginal significance threshold (z ≈ 1.64).** There is a sharp increase in the frequency of test statistics just above the threshold for marginal significance (α = 0.10).
2. **Zero threshold (z = 0).** A weaker discontinuity occurs at zero, with additional suppression of studies reporting negative effects (z < 0).

The random-effects model (black dashed line) fails to capture these patterns.
It systematically overestimates the number of negative results and non-significant positive results.

RoBMA (blue dashed line) captures both discontinuities and approximates the observed data much better.
These results provide extreme evidence for the presence of publication bias and highlight the need to interpret the publication-bias-adjusted model.

**Extrapolation to pre-publication bias.**
The package also allows us to extrapolate what the distribution of z-statistics might have looked like in the absence of publication bias.
This is achieved by calling the `zplot()` function (with the default `plot_extrapolation = TRUE` argument), which builds the diagnostic plot directly from a fitted model.

```{r, fig.cap="Ease-of-Retrieval Effect: Extrapolation Analysis"}
zplot(fit_RoBMA_Weingarten2018, from = -3, to = 6, by.hist = 0.25)
```

The extrapolated distribution (blue line) shows what we would expect to observe if studies were published regardless of their results.
Comparing the fitted distribution (accounting for publication bias) with the extrapolated distribution reveals the extent of the bias.
The large discrepancy between these distributions quantifies the substantial impact of publication bias in this literature.

**Zplot summary metrics.**
This discrepancy can be summarized with the additional statistics provided by the `summary()` function applied to a `zplot_brma` object.

```{r}
summary(as_zplot(fit_RoBMA_Weingarten2018))
```

The summary provides several key results.
The observed discovery rate (ODR) substantially exceeds the expected discovery rate (EDR), indicating that many more significant results appear in the published literature than we would expect.
The estimated number of missing studies suggests that a substantial number of non-significant or negative results may be absent from the published literature.
The false discovery risk (FDR) provides an upper bound on the proportion of statistically significant results that may be false positives, though this risk remains moderate due to evidence for a genuine underlying effect despite the extreme publication bias.

### Example 2: Social Comparison and Behavior Change

This example examines a meta-analysis of randomized controlled trials evaluating the efficacy of social comparison as a behavior change technique [@hoppen2025meta].
The analysis includes 37 trials comparing social comparison interventions to passive controls across domains including climate-change mitigation, health, performance, and service outcomes.

**Data and model fitting.**
Again, we begin by loading the social comparison dataset and examining its structure:

```{r}
data("Hoppen2025", package = "RoBMA")
head(Hoppen2025)
```

The dataset contains effect sizes (`d`) and sampling variances (`v`) from individual studies.
We fit both a random-effects model using `brma()` and a publication-bias-adjusted model using `RoBMA()`:

```{r, eval = FALSE}
fit_RE_Hoppen2025 <- brma(
  yi = d, sei = sqrt(v), measure = "SMD",
  data = Hoppen2025,
  sample = 10000, burnin = 5000, adapt = 5000,
  chains = 5, parallel = TRUE, seed = 1)

fit_RoBMA_Hoppen2025 <- RoBMA(
  yi = d, sei = sqrt(v), measure = "SMD",
  data = Hoppen2025,
  sample = 10000, burnin = 5000, adapt = 5000,
  chains = 5, parallel = TRUE, seed = 1)
```

**Model results.**
We examine the key results from both models:

```{r}
summary(fit_RE_Hoppen2025)
```

The random-effects model estimates a positive effect size, suggesting a positive effect of social comparison interventions.
However, this analysis does not account for potential publication bias.

```{r}
summary(fit_RoBMA_Hoppen2025)
```

When accounting for the possibility of publication bias with RoBMA, we find extreme evidence for the presence of publication bias.
Consequently, the effect size shrinks to essentially zero, with moderate evidence against the presence of an effect.

**Zplot diagnostics.**

```{r, fig.cap="Social Comparison: Model Fit Assessment"}
hist(as_zplot(fit_RoBMA_Hoppen2025))

lines(as_zplot(fit_RE_Hoppen2025),    col = "black", lty = 2, lwd = 2)
lines(as_zplot(fit_RoBMA_Hoppen2025), col = "blue",  lty = 2, lwd = 2)

legend(
  "topright",
  legend = c("Random-Effects", "RoBMA"),
  col = c("black", "blue"),
  lty = 2, lwd = 2)
```

The zplot highlights the clear evidence of publication bias in this dataset.
We can observe pronounced discontinuities in the observed distribution (gray bars) at critical thresholds; particularly at the transition to marginal significance (z ≈ 1.64) and at zero (indicating selection against negative results).

The random-effects model (black dashed line) fails to capture these patterns, systematically overestimating the number of negative and non-significant results.
In contrast, RoBMA (blue dashed line) successfully models both discontinuities, providing a markedly better fit to the observed data.
This visual assessment aligns with the statistical evidence: RoBMA yields extreme evidence for publication bias and suggests that the unadjusted pooled effect is misleading.

**Extrapolation to pre-publication bias.**
The package also allows us to extrapolate what the distribution of z-statistics might have looked like in the absence of publication bias.
This is achieved by either setting `extrapolate = TRUE` in the `lines()` function, or calling the `zplot()` function (with the default `plot_extrapolation = TRUE` argument).

```{r, fig.cap="Social Comparison: Extrapolation Analysis"}
zplot(fit_RoBMA_Hoppen2025)
```

The extrapolated distribution (blue line) shows what we would expect to observe if studies were published regardless of their results.
Comparing the fitted distribution (accounting for publication bias) with the extrapolated distribution reveals the extent of the bias.
The large discrepancy between these distributions quantifies the substantial impact of publication bias in this literature, with implications for the estimated effect size and number of missing studies.

**Zplot summary metrics.**
This discrepancy can be summarized with the additional statistics provided by the `summary()` function.

```{r}
summary(as_zplot(fit_RoBMA_Hoppen2025))
```

The summary provides several results.
The observed discovery rate (ODR, i.e., the observed proportion of statistically significant results) in the dataset matches the expected discovery rate (EDR) due to the one-sided selection for marginally significant results (instead of statistically significant results).
The estimated number of missing studies suggests that a considerable number of non-significant or negative results may be absent from the published literature.
The false discovery risk (FDR), which corresponds to the upper bound on the false positive rate, is not extremely inflated due to the possible small positive and negative effects under the moderate heterogeneity.


### Example 3: ChatGPT and Learning Performance

This example examines the effectiveness of ChatGPT-based interventions on students' learning performance [@wang2025effect].
This meta-analysis includes 42 randomized controlled trials comparing experimental groups (using ChatGPT for tutoring or learning support) with control groups (without ChatGPT) on learning outcomes such as exam scores and final grades.

**Data and model fitting.**
We follow the same procedure as in the previous example:

```{r}
data("Wang2025", package = "RoBMA")
Wang2025 <- Wang2025[Wang2025$Learning_effect == "Learning performance", ]
head(Wang2025)
```

```{r, eval = FALSE}
fit_RE_Wang2025 <- brma(
  yi = g, sei = se, measure = "SMD",
  data = Wang2025,
  sample = 10000, burnin = 5000, adapt = 5000,
  chains = 5, parallel = TRUE, seed = 1)
fit_RoBMA_Wang2025 <- RoBMA(
  yi = g, sei = se, measure = "SMD",
  data = Wang2025,
  sample = 10000, burnin = 5000, adapt = 5000,
  chains = 5, parallel = TRUE, seed = 1)
```

**Zplot diagnostics.**

```{r, fig.cap="ChatGPT: Model Fit Assessment"}
hist(as_zplot(fit_RoBMA_Wang2025), from = -2, to = 8)
lines(as_zplot(fit_RE_Wang2025),    col = "black", lty = 2, lwd = 2, from = -2, to = 8)
lines(as_zplot(fit_RoBMA_Wang2025), col = "blue",  lty = 2, lwd = 2, from = -2, to = 8)
legend(
  "topright",
  legend = c("Random-Effects", "RoBMA"),
  col = c("black", "blue"),
  lty = 2, lwd = 2)
```

The zplot for the ChatGPT data shows a different pattern than the extreme publication bias observed in the social comparison example.
While we do not see strong selection at conventional significance thresholds, there is a moderate discontinuity at the transition to non-conforming results (z = 0), suggesting some degree of selection against negative findings.

The random-effects model (black dashed line) provides a better fit to the data than in the previous example; however, RoBMA (blue dashed line) captures the discontinuity at zero slightly better.
This visual pattern corresponds to moderate statistical evidence for publication bias, highlighting a case where both models might be considered, though the RoBMA model incorporates the uncertainty about the best model and provides a more complete account of the data patterns.

**Extrapolation to pre-publication bias.**
We can examine the extrapolation to assess the impact of publication bias:

```{r, fig.cap="ChatGPT: Extrapolation Analysis"}
zplot(fit_RoBMA_Wang2025, from = -2, to = 8)
```

The extrapolated distribution (blue line) shows a more modest difference between the fitted and extrapolated distributions compared to the extreme bias example, reflecting the moderate degree of publication bias in this literature.

**Model results.**
To quantify these visual patterns, we examine the model summaries:

```{r}
summary(fit_RE_Wang2025)
summary(fit_RoBMA_Wang2025)
```

The random-effects model yields a substantial positive effect size estimate, while the RoBMA model accounting for selection produces a more conservative estimate with a wider credible interval.
Importantly, the results show an extreme degree of between-study heterogeneity that greatly complicates feasible implications and recommendations.
This demonstrates the moderate nature of the publication bias, where the adjusted estimate is meaningfully smaller but not completely reduced.
The Bayes factor for publication bias provides moderate evidence for publication bias, and the evidence for the effect becomes weak.

**Zplot summary metrics.**
This moderate publication-bias pattern is reflected in the summary statistics,

```{r}
summary(as_zplot(fit_RoBMA_Wang2025))
```

which show a moderate-to-high EDR and a small number of missing estimates.

### Example 4: Framing Effects from Many Labs 2

Our final example analyzes registered replication reports of the classic framing effect on decision making [@tversky1981framing] conducted as part of the Many Labs 2 project [@klein2018many].
This dataset provides an ideal test case for zplot diagnostics because the pre-registered nature of these studies does not allow for publication bias.
The analysis includes 55 effect size estimates that examine how framing influences decision-making preferences.

**Data and model fitting.**

```{r}
data("ManyLabs16", package = "RoBMA")
head(ManyLabs16)
```

```{r, eval = FALSE}
fit_RE_ManyLabs16 <- brma(
  yi = y, sei = se, measure = "SMD",
  data = ManyLabs16,
  sample = 10000, burnin = 5000, adapt = 5000,
  chains = 5, parallel = TRUE, seed = 1)
fit_RoBMA_ManyLabs16 <- RoBMA(
  yi = y, sei = se, measure = "SMD",
  data = ManyLabs16,
  sample = 10000, burnin = 5000, adapt = 5000,
  chains = 5, parallel = TRUE, seed = 1)
```

**Zplot diagnostics.**

```{r, fig.cap="Framing Effects: Model Fit Assessment"}
hist(as_zplot(fit_RoBMA_ManyLabs16))
lines(as_zplot(fit_RE_ManyLabs16),    col = "black", lty = 2, lwd = 2)
lines(as_zplot(fit_RoBMA_ManyLabs16), col = "blue",  lty = 2, lwd = 2)
legend(
  "topleft",
  legend = c("Random-Effects", "RoBMA"),
  col = c("black", "blue"),
  lty = 2, lwd = 2)
```

The zplot for the Many Labs 2 framing effects demonstrates what we expect to see in the absence of publication bias.
The observed distribution of z-statistics (gray bars) appears smooth without sharp discontinuities at significance thresholds or at zero.
Both the random-effects model (black dashed line) and RoBMA (blue dashed line) provide essentially identical fits to the data, with their posterior predictive distributions overlapping almost perfectly.

This close agreement between models indicates that either approach would be appropriate for these data.
The absence of publication bias is further confirmed by the statistical evidence: RoBMA provides moderate evidence against publication bias, demonstrating how the method appropriately penalizes unnecessary model complexity when simpler models explain the data equally well.

**Extrapolation to pre-publication bias.**
We can examine whether there would be any difference in the absence of publication bias:

```{r, fig.cap="Framing Effects: Extrapolation Analysis"}
zplot(fit_RoBMA_ManyLabs16)
```

The extrapolated distribution (blue line) shows virtually no difference between the fitted and extrapolated distributions, confirming that publication bias has minimal impact in this well-designed replication project.
This example illustrates the ideal scenario where traditional meta-analytic approaches are fully justified.

**Model results.**
The quantitative results confirm the visual impression:

```{r, eval = FALSE}
summary(fit_RE_ManyLabs16)
summary(fit_RoBMA_ManyLabs16)
```

Both models yield virtually identical effect size estimates.
The Bayes factor for publication bias provides moderate evidence against the presence of publication bias, appropriately penalizing the more complex model when it offers no advantage.
This demonstrates the method's ability to distinguish between necessary and unnecessary model complexity.

**Zplot summary metrics.**
The absence of publication bias is reflected in the publication-bias assessment statistics: a moderate EDR matching the ODR and no missing studies.

```{r}
summary(as_zplot(fit_RoBMA_ManyLabs16))
```


## Conclusions

Zplots are an intuitive diagnostic tool for assessing publication bias and model fit [@bartos2025zcurve].
By visualizing the distribution of test statistics and comparing observed patterns with model-implied expectations, researchers can make more informed decisions about their analytic approach using the `RoBMA` R package [@RoBMA].

Zplot diagnostics are particularly informative when applied to moderate to large meta-analyses (typically >20-30 studies), where histogram patterns become interpretable.
Publication bias and questionable research practices (QRPs) might produce similar patterns of results.
The zplot diagnostics cannot distinguish between them; however, they might help in assessing whether the model approximates the observed data well.
They are especially useful for model comparison scenarios when they provide a visual supplement to statistical tests like inclusion Bayes factors.

The following points are important for interpreting zplot diagnostics:

 - Sharp drops in the observed distribution at z = 0, z = ±1.64, or z = ±1.96 suggest publication bias.
 - Models whose posterior predictive distributions closely match the observed pattern should be preferred.
 - Large differences between fitted and extrapolated distributions indicate substantial publication bias.
 - Visual assessment should complement formal statistical tests.

## References
