data(dietaryfat)
# Transform data from arm-based format to contrast-based format
p1 <- pairwise(list(treat1, treat2, treat3),
list(d1, d2, d3),
time=list(years1, years2, years3),
studlab=ID,
data=dietaryfat)
p1
# Conduct network meta-analysis using incidence rate ratios (sm="IRR").
# Note, the argument 'sm' is not necessary as this is the default in R
# function metainc called internally
net1 <- netmeta(TE, seTE, treat1, treat2, studlab, data=p1, sm="IRR")
summary(net1)
# Conduct network meta-analysis using incidence rate differences
# (sm="IRD").
net2 <- netmeta(TE, seTE, treat1, treat2, studlab, data=p1, sm="IRD")
summary(net2)
# Draw network graph
netgraph(net1, points=TRUE, cex.points=3, cex=1.25)
tname <- c("Control","Diet", "Diet 2")
netgraph(net1, points=TRUE, cex.points=3, cex=1.25, labels=tname)
Run the code above in your browser using DataLab