Draw a ‘comparison-adjusted’ radial plot to assess funnel plot asymmetry in network meta-analysis.
# S3 method for netmeta
radial(
x,
order,
pooled = ifelse(x$random, "random", "common"),
level = x$level,
pch,
col = "black",
legend = FALSE,
pos.legend = "topright",
pos.tests = "topleft",
lump.comparator = FALSE,
text.comparator = "comparator",
method.bias,
text.linreg = "(Egger)",
text.rank = "(Begg-Mazumdar)",
text.mm = "(Thompson-Sharp)",
sep.trts = x$sep.trts,
nchar.trts = x$nchar.trts,
digits.pval = gs("digits.pval"),
warn.deprecated = gs("warn.deprecated"),
linreg = FALSE,
rank = FALSE,
mm = FALSE,
...
)A data frame with the following columns:
Study label.
Label/Number for first treatment.
Label/Number for second treatment.
Treatment comparison.
Estimate of treatment effect, e.g., log odds ratio.
Pooled estimate from direct evidence.
‘Comparison-adjusted’ treatment effect (TE - TE.direct).
Standard error of treatment estimate.
Plotting symbol(s).
Colour of plotting symbol(s).
An object of class netmeta.
A mandatory character or numerical vector specifying the order of treatments or list of comparators (see Details).
A character string indicating whether results for the
common ("common") or random effects model ("random")
should be plotted. Can be abbreviated.
The confidence level utilised in the plot.
The plotting symbol(s) used for individual studies within direct comparisons.
The colour(s) used for individual studies within direct comparisons.
A logical indicating whether a legend with information on direct comparisons should be added to the plot.
The position of the legend describing plotting symbols and colours for direct comparisons.
The position of results for test(s) of funnel plot asymmetry.
A logical indicating whether comparators should be lumped, e.g., to specify inactive treatments. information on direct comparisons should be added to the plot.
A character string used in the plot to label
the comparator if lump.comparator is TRUE.
A character vector indicating which test(s) for
funnel plot asymmetry to use. Admissible values are
"Begg", "Egger", and "Thompson", can be
abbreviated. See function metabias.
A character string used in the plot to label the Egger test for funnel plot asymmetry.
A character string used in the plot to label the Begg test for funnel plot asymmetry.
A character string used in the plot to label the Thompson-Sharp test for funnel plot asymmetry.
A character used in comparison names as separator between treatment labels.
A numeric defining the minimum number of
characters used to create unique treatment names (see
netmeta).
Minimal number of significant digits for p-value of test(s) for funnel plot asymmetry.
A logical indicating whether warnings should be printed if deprecated arguments are used.
Deprecated argument (replaced by method.bias).
Deprecated argument (replaced by method.bias).
Deprecated argument (replaced by method.bias).
Additional graphical arguments passed as arguments to
radial.meta.
Guido Schwarzer guido.schwarzer@uniklinik-freiburg.de
A ‘comparison-adjusted’ radial plot is drawn in the active graphics window. The idea of this radial plot is similar to a ‘comparison-adjusted’ funnel plot (Chaimani & Salanti, 2012).
Argument order is mandatory to determine the order of
treatments (Chaimani et al., 2013):
“Before using this plot, investigators should order the treatments in a meaningful way and make assumptions about how small studies differ from large ones. For example, if they anticipate that newer treatments are favored in small trials, then they could name the treatments from oldest to newest so that all comparisons refer to ‘old versus new intervention’. Other possibilities include defining the comparisons so that all refer to an active treatment versus placebo or sponsored versus non-sponsored intervention.”
Alternatively, it is possible to only provide a single or few
treatment name(s) in argument order to define the
comparator(s). In this case only comparisons with this / these
treatment(s) will be considered. If argument lump.comparator
is TRUE, all comparators will be lumped into a single
group. The text for this group can be specified with argument
text.comparator.
Chaimani A & Salanti G (2012): Using network meta-analysis to evaluate the existence of small-study effects in a network of interventions. Research Synthesis Methods, 3, 161--76
Chaimani A, Higgins JP, Mavridis D, Spyridonos P, Salanti G (2013): Graphical tools for network meta-analysis in STATA. PLOS ONE, 8, e76654
netmeta, radial.meta,
funnel.meta, metabias
# \donttest{
data(Senn2013)
net1 <- netmeta(TE, seTE, treat1, treat2, studlab,
data = Senn2013, sm = "MD")
# 'Comparison-adjusted' radial plot not created as argument 'order'
# is missing
#
try(radial(net1))
# Only show comparisons with placebo
#
radial(net1, order = "pl")
# Add result for Egger test of funnel plot asymmetry
#
radial(net1, order = "pl", method.bias = "Egger",
digits.pval = 2)
# (Non-sensical) alphabetic order of treatments with placebo as
# last treatment
#
ord <- c("a", "b", "me", "mi", "pi", "r", "si", "su", "v", "pl")
radial(net1, order = ord)
# Add results for tests of funnel plot asymmetry and use different
# plotting symbols and colours
#
radial(net1, order = ord,
pch = rep(c(15:18, 1), 3), col = 1:3,
method.bias = c("Egger", "Begg", "Thompson"), digits.pval = 2)
# Same results for tests of funnel plot asymmetry using reversed
# order of treatments
#
radial(net1, order = rev(ord),
pch = rep(c(15:18, 1), 3), col = 1:3,
method.bias = c("Egger", "Begg", "Thompson"), digits.pval = 2)
# Calculate tests for funnel plot asymmetry
#
f1 <- radial(net1, order = ord)
#
metabias(metagen(TE.adj, seTE, data = f1))
metabias(metagen(TE.adj, seTE, data = f1), method = "Begg")
metabias(metagen(TE.adj, seTE, data = f1), method = "Thompson")
# }
Run the code above in your browser using DataLab