data(smoke)
# increase n.iter to reach convergence
nma.out <- nma.ab(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)
rank.prob(nma.out)Run the code above in your browser using DataLab