pcnetmeta (version 2.8)

rank.prob: Plotting Treatment Rank Probabilities

Description

rank.prob generates a plot for treatment rank probabilities.

Usage

rank.prob(nma.obj, cex.axis = 1, cex.lab = 1, median = FALSE)

Value

None.

Arguments

nma.obj

a list object obtained by function nma.ab.bin, nma.ab.cont, nma.ab.py, or nma.ab.followup.

cex.axis

a numeric value specifying the size of the tick label numbers/text.

cex.lab

a numeric value specifying the size of the axis label text.

median

a logical value indicating whether the rank probabilities are based on the marginal effect measures (FALSE, the default) or the overall median odds ratios. This argument is only used when link = "logit" in nma.ab.bin.

Details

A plot of treatment rank probabilities is generated. Each vertical bar represents probabilities of being different ranks for a specific treatment. A darker area indicates the probability of being a higher rank. The black area indicates the probability of being the best treatment.

Examples

data("smoke")
# increase n.iter to reach convergence
set.seed(1234)
nma.out <- nma.ab.bin(s.id, t.id, r, n, data = smoke,
  trtname = c("NC", "SH", "IC", "GC"), param = "rank.prob", model = "het_cor",
  higher.better = TRUE, n.adapt = 1000, n.iter = 100, n.chains = 1)
rank.prob(nma.out)