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 refdist.

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 observed = FALSE.

...

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)


[Package poth version 0.3-0 Index]