1 Introduction

The nmaplateplot package is a graphical display of results from network meta-analysis (NMA). It is suitable for outcomes like odds ratio (OR), risk ratio (RR), risk difference (RD) and standardized mean difference (SMD). It also has an option to visually display and compare the surface under the cumulative ranking (SUCRA) of different treatments. Moreover, nmaplateplot is good at details, including choosing colors, text size, circle size, layout and adding title.

1.1 Quick example: Efficacy and acceptability of 12 Antidepressants

We use a dataset to quickly demonstrate the plate plot.

The ad12.eff.acc dataset contains the point estimates (posterior medians), interval estimates (95% credible intervals), p-values of odds ratios (ORs) for both efficacy and acceptability of 12 antidepressants based on arm-based network meta-analysis.

library(nmaplateplot)
library(svglite)
## Warning: package 'svglite' was built under R version 4.3.3
data("ad12.eff.acc")
plateplot(ad12.eff.acc, design_method = c("circle", "circle"),
             upper_diagonal_name = "Efficacy", lower_diagonal_name = "Acceptability")

Let us explain the components of the plate plot:

The upper diagonal part graphically presents ORs for efficacy, while the lower diagonal part graphically presents ORs for acceptability.

Each cell in off-diagonal part contains 2 (for p > 0.05) or 3 (for p < 0.05) concentric circles with white circle (inner circle) representing lower bound of interval estimates if result is significant, grey circle (middle circle) representing point estimates, and outer circle with varying colors representing upper bound of interval estimates. The varying colors are associated with p-values shown in the color bar. The circles’ size is positively associated with (point/interval) estimates.

For example: 1) the plate with (blue/grey/white) color in the 1st row and 12th column cell means MIR is significantly more effective than REB in terms of response rate (efficacy); 2) the plate with (red/grey) color in the 2nd row and 1st column cell means MIR is less acceptable (higher dropout rate) than ESC, though this effect is not statistically significant.

The diagonal texts are names of 12 antidepressant drugs with 2 numbers in each cell standing for SUCRA(%) for efficacy and acceptability of these treatments. The background colors in diagonal part are associated with SUCRA values as shown in the color bar.

In addition to the plate plot which graphically and intuitively displays the NMA results, the package can directly print the point estimates with 95% credible intervals.

# option design_method controls whether you want upper and lower diagonal parts to display "text" or "circle" respectively.
plateplot(ad12.eff.acc, design_method = c("text", "text"),
             upper_diagonal_name = "Efficacy", lower_diagonal_name = "Acceptability")

The estimates with p-value < 0.05 are underscored. Symbols “a”, “b” and “c” are added to estimates with p-value < 0.05, < 0.01 and < 0.001 respectively. The varying colors of text are associated with p-values shown in the color bar.

For example: 1) 1.91 (1.41, 2.56) in the 1st row and 12th column cell prints odds of efficacy in the upper-left defined treatment (MIR) compared with odds of efficacy in the lower-right defined treatment (REB); 2) 1.27 (0.98, 1.66) in the 2nd row and 1st column cell prints odds of acceptability in the upper-left defined treatment (MIR) compared with odds of acceptability in the lower-right defined treatment (ESC).

1.2 Input dataset

The input dataset must be a list containing at least 5 data frames named: Point_estimates, Interval_estimates_LB, Interval_estimates_UB, Pvalues, and Treatment_specific_values.

nma_res_minimum <- ad12.eff.acc
nma_res_minimum$Treatment_specific_values$Value_Upper <- NULL
nma_res_minimum$Treatment_specific_values$Value_Lower <- NULL
# This is the minimum input for the plateplot function
nma_res_minimum$Treatment_specific_values
##    Trt_ID Trt_abbrv
## 1       1       BUP
## 2       2       CIT
## 3       3       DUL
## 4       4       ESC
## 5       5       FLU
## 6       6       FVX
## 7       7       MIL
## 8       8       MIR
## 9       9       PAR
## 10     10       REB
## 11     11       SER
## 12     12       VEN
round(nma_res_minimum$Point_estimates,2)
##      V1   V2   V3   V4   V5   V6   V7   V8   V9  V10  V11  V12
## 1  0.00 0.99 1.13 0.83 1.10 1.13 1.18 0.80 1.07 1.53 0.89 0.85
## 2  0.89 0.00 1.14 0.84 1.11 1.14 1.19 0.80 1.08 1.54 0.89 0.86
## 3  1.24 1.39 0.00 0.74 0.97 1.00 1.04 0.71 0.95 1.35 0.78 0.75
## 4  0.87 0.97 0.70 0.00 1.32 1.36 1.42 0.96 1.29 1.83 1.06 1.02
## 5  1.09 1.23 0.88 1.26 0.00 1.03 1.08 0.73 0.98 1.39 0.81 0.78
## 6  1.27 1.42 1.02 1.47 1.16 0.00 1.04 0.71 0.95 1.35 0.78 0.75
## 7  1.25 1.40 1.01 1.45 1.14 0.98 0.00 0.68 0.91 1.30 0.75 0.72
## 8  1.10 1.23 0.89 1.27 1.01 0.86 0.88 0.00 1.35 1.92 1.11 1.07
## 9  1.18 1.32 0.95 1.36 1.08 0.93 0.94 1.07 0.00 1.43 0.83 0.80
## 10 1.42 1.58 1.15 1.64 1.29 1.11 1.13 1.29 1.20 0.00 0.58 0.56
## 11 0.93 1.04 0.75 1.08 0.85 0.73 0.74 0.85 0.79 0.66 0.00 0.96
## 12 1.15 1.29 0.93 1.33 1.05 0.90 0.92 1.04 0.98 0.81 1.23 0.00
round(nma_res_minimum$Interval_estimates_LB,2)
##      V1   V2   V3   V4   V5   V6   V7   V8   V9  V10  V11  V12
## 1  0.00 0.79 0.87 0.68 0.92 0.86 0.85 0.63 0.88 1.15 0.73 0.71
## 2  0.70 0.00 0.88 0.70 0.93 0.88 0.87 0.64 0.90 1.18 0.73 0.71
## 3  0.91 1.02 0.00 0.59 0.77 0.74 0.73 0.54 0.76 0.98 0.61 0.59
## 4  0.70 0.79 0.54 0.00 1.13 1.05 1.04 0.76 1.09 1.40 0.89 0.86
## 5  0.91 1.00 0.68 1.06 0.00 0.82 0.81 0.60 0.85 1.09 0.69 0.68
## 6  0.94 1.05 0.72 1.10 0.89 0.00 0.74 0.54 0.74 0.98 0.61 0.58
## 7  0.90 0.99 0.69 1.05 0.86 0.69 0.00 0.49 0.68 0.89 0.55 0.53
## 8  0.85 0.94 0.64 0.98 0.80 0.65 0.61 0.00 1.10 1.42 0.90 0.87
## 9  0.96 1.06 0.73 1.12 0.93 0.71 0.68 0.85 0.00 1.09 0.70 0.67
## 10 1.03 1.18 0.79 1.20 0.99 0.78 0.75 0.92 0.90 0.00 0.44 0.42
## 11 0.74 0.83 0.55 0.87 0.71 0.55 0.53 0.66 0.65 0.48 0.00 0.81
## 12 0.94 1.02 0.69 1.09 0.90 0.69 0.66 0.83 0.81 0.60 1.01 0.00
round(nma_res_minimum$Interval_estimates_UB,2)
##      V1   V2   V3   V4   V5   V6   V7   V8   V9  V10  V11  V12
## 1  0.00 1.24 1.47 1.02 1.31 1.49 1.63 1.01 1.30 2.06 1.07 1.02
## 2  1.15 0.00 1.47 1.01 1.31 1.48 1.63 1.01 1.30 2.03 1.08 1.05
## 3  1.69 1.89 0.00 0.91 1.21 1.36 1.46 0.93 1.17 1.87 0.99 0.95
## 4  1.08 1.20 0.91 0.00 1.55 1.77 1.94 1.19 1.52 2.44 1.27 1.22
## 5  1.33 1.50 1.17 1.52 0.00 1.30 1.43 0.88 1.12 1.79 0.93 0.89
## 6  1.73 1.92 1.50 1.97 1.51 0.00 1.47 0.91 1.20 1.86 1.00 0.96
## 7  1.79 2.02 1.49 2.05 1.56 1.45 0.00 0.93 1.21 1.86 1.02 0.98
## 8  1.44 1.62 1.23 1.64 1.25 1.17 1.24 0.00 1.64 2.61 1.38 1.31
## 9  1.46 1.64 1.25 1.65 1.24 1.21 1.27 1.34 0.00 1.87 0.98 0.94
## 10 1.97 2.15 1.64 2.24 1.71 1.62 1.69 1.81 1.62 0.00 0.76 0.73
## 11 1.16 1.31 1.02 1.33 1.01 0.97 1.02 1.08 0.96 0.89 0.00 1.14
## 12 1.42 1.62 1.26 1.63 1.22 1.19 1.26 1.32 1.17 1.09 1.51 0.00
round(nma_res_minimum$Pvalues,3)
##       V1    V2    V3    V4    V5    V6    V7    V8    V9   V10   V11   V12
## 1  0.000 0.948 0.353 0.073 0.304 0.367 0.327 0.062 0.491 0.006 0.207 0.087
## 2  0.390 0.000 0.307 0.060 0.251 0.316 0.290 0.061 0.417 0.002 0.247 0.130
## 3  0.165 0.033 0.000 0.006 0.800 0.995 0.803 0.012 0.628 0.068 0.044 0.020
## 4  0.196 0.770 0.008 0.000 0.001 0.020 0.030 0.705 0.004 0.000 0.509 0.797
## 5  0.341 0.046 0.357 0.010 0.000 0.795 0.613 0.001 0.732 0.010 0.004 0.000
## 6  0.116 0.023 0.902 0.012 0.265 0.000 0.812 0.007 0.660 0.066 0.053 0.024
## 7  0.178 0.055 0.966 0.027 0.351 0.932 0.000 0.017 0.515 0.164 0.066 0.036
## 8  0.462 0.125 0.468 0.067 0.956 0.345 0.463 0.000 0.004 0.000 0.336 0.522
## 9  0.120 0.010 0.693 0.002 0.307 0.574 0.689 0.557 0.000 0.009 0.026 0.007
## 10 0.034 0.003 0.472 0.002 0.064 0.564 0.551 0.139 0.215 0.000 0.000 0.000
## 11 0.534 0.733 0.066 0.505 0.066 0.032 0.067 0.179 0.017 0.006 0.000 0.671
## 12 0.175 0.034 0.598 0.006 0.527 0.470 0.596 0.727 0.789 0.178 0.036 0.000
# The text output based on minimum input
plateplot(nma_res_minimum, design_method = c("text", "text"),
             upper_diagonal_name = "Efficacy", lower_diagonal_name = "Acceptability")