plot.refdist {poth} | R Documentation |
Plot reference distribution for given network
Description
Plot a histogram of simulated reference distribution for POTH
Usage
## S3 method for class 'refdist'
plot(x, observed = TRUE, probability = FALSE, ...)
Arguments
x |
R object of class |
observed |
A logical indicating whether a line with the observed POTH should be drawn in the plot. |
probability |
A logical indicating whether the probability of a POTH
greater than the observed POTH should be printed in the legend. This
argument is ignored if |
... |
Additional arguments (ignored). |
Details
Plots a histogram of the simulated POTH values from the reference distribution. (Wigle et al., 2024).
Value
A ggplot2 object.
Author(s)
Augustine Wigle amhwigle@uwaterloo.ca, Guido Schwarzer guido.schwarzer@uniklinik-freiburg.de
References
Wigle A, Béliveau A, Salanti G, Rücker G, Schwarzer G, Mavridis D, Nikolakopoulou A (2024): Precision of treatment hierarchy: A metric for quantifying uncertainty in treatment hierarchies in network meta-analysis. Preprint on arXiv, doi:10.48550/arXiv.2501.11596
Examples
library("netmeta")
data(smokingcessation)
p1 <- pairwise(list(treat1, treat2, treat3),
event = list(event1, event2, event3), n = list(n1, n2, n3),
data = smokingcessation, sm = "OR")
net1 <- netmeta(p1, random = FALSE)
# Reference distribution
dist1 <- refdist(net1, pooled = "common", nsim = 2)
plot(dist1, observed = FALSE)
data(Senn2013)
net2 <- netmeta(TE, seTE, treat1.long, treat2.long, studlab,
data = Senn2013, sm = "MD", random = FALSE)
# Reference distribution
dist2 <- refdist(net2, pooled = "common", nsim = 25)
plot(dist2, probability = TRUE)