gemtc (version 0.1-2)

rank.probability: Calculating rank-probabilities

Description

Rank probabilities indicate the probability for each treatment to be best, second best, etc.

Usage

rank.probability(result)

Arguments

result

Object of S3 class mtc.result to be used in creation of the rank probability table

Value

• A matrix with the treatments as rows and the ranks as columns

encoding

utf8

Details

For each MCMC iteration, the treatments are ranked by their effect relative to an arbitrary baseline. A frequency table is constructed from these rankings and normalized by the number of iterations to give the rank probabilities.

See Also

relative.effect

Examples

file <- system.file("extdata/luades-smoking.gemtc", package="gemtc")
network <- read.mtc.network(file)
model <- mtc.model(network)
# To save computation time we load the samples instead of running the model
results <- mtc.run(model)
file.samples <- system.file("extdata/luades-smoking.samples.gz", package="gemtc")
results <- read.mtc.result.samples(file.samples, model, "rjags")

ranks <- rank.probability(results)
print(ranks)
##       [,1]      [,2]     [,3]     [,4]
## A 0.0000500 0.0024250 0.104000 0.893525
## B 0.0571125 0.1802375 0.661125 0.101525
## C 0.2409625 0.5893875 0.169375 0.000275
## D 0.7020250 0.2279000 0.065450 0.004625
barplot(t(ranks), beside=TRUE)