Code
source("R/nma_functions.R")
dat <- load_pairwise_data()
nm <- run_model(dat, model_type = "simple", effect_model = "random")

P-score ranking

Code
ranking <- tryCatch(extract_ranking(nm, effect_model = "random"), error = identity)
if (inherits(ranking, "error")) {
  cat("Ranking failed:", conditionMessage(ranking))
} else {
  knitr::kable(ranking)
}
Treatment P_score
Clb 0.967
Clb + O 0.919
FCR 0.857
FCR + BR 0.788
BR 0.712
Ven + R 0.642
I + R 0.511
Ven + O 0.472
Ibr + Ven 0.384
Ibr + O 0.343
Acal 0.316
I 0.219
Ven + O + Ibr 0.216
Zanubrutinib 0.103
Acal + O 0.051