Calculate SUCRA and MNBT based on MCMC draws in NMA or RaCE-NMA models

Source: R/calculate_SUCRA_MNBT.R

This function creates a results table of SUCRA and Median Number of Better Treatments (MNBT) based on MCMC draws from a standard NMA model or a fitted RaCE-NMA model.

calculate_SUCRA_MNBT(data = NULL, mcmc = NULL, level = 0.5, names = NULL)

Arguments

data

A NxJ matrix of data, where N is the number of observations and J the number of treatments. This feature is designed to display results from a standard NMA study.

mcmc

MCMC draws from the RaCE NMA model, in the form of the model output of the mcmc_RCMVN function.

level

The level level for the interval estimate of MNBT. Defaults to 0.50.

names

A vector of intervention names (optional)

Value

A table containing SUCRA and MNBT values for each treatment, ordered by descending SUCRA values.

Examples

mcmc <- mcmc_raceNMA(mu_hat=c(0,0,1,1), s=c(.1,.1,.1,.1), seed=1)
#> [1] "Estimating chain 1 of 2."
#> [1] "Estimating chain 2 of 2."
calculate_SUCRA_MNBT(mcmc=mcmc)
#>   Treatment     SUCRA MNBT (50% CI)
#> 1         1 0.9876833      0 (0, 0)
#> 2         2 0.9874000      0 (0, 0)
#> 3         4 0.3202167      2 (2, 2)
#> 4         3 0.3200333      2 (2, 2)