data(smokingcessation)
# Add variable with (fictious) risk of bias values
# with 1 = "low risk" and 2 = "high risk"
#
smokingcessation$rob <- rep(1:2, 12)
p1 <- pairwise(list(treat1, treat2, treat3),
event = list(event1, event2, event3), n = list(n1, n2, n3),
data = smokingcessation, sm = "OR")
net1 <- netmeta(p1, common = FALSE, ref = "A")
# Generate network graph with information on risk of bias
#
col.rob <- netmatrix(net1, rob, ties.method = "last",
levels = 1:2, labels = c("green", "yellow"))
#
netgraph(net1, plastic = FALSE, col = col.rob,
cex.points = 5, bg.points = "gray", adj = 0.5)
netgraph(net1, plastic = FALSE, col = col.rob,
cex.points = n.trts, bg.points = "blue",
labels = paste0(trts, " (n=", n.trts, ")"),
offset = c(0.05, 0.035, 0.05, 0.025))
Run the code above in your browser using DataLab