data(Senn2013)
#
# Fixed effect model (default)
#
net1 <- netmeta(TE, seTE, treat1, treat2, studlab,
data=Senn2013, sm="MD")
forest(net1, ref="plac")
forest(net1, xlim=c(-1.5,1), ref="plac",
xlab="HbA1c difference",
leftcols="studlab", rightcols=FALSE,
leftlabs="Contrast to placebo")
#
# Random effects effect model
#
net2 <- netmeta(TE, seTE, treat1, treat2, studlab,
data=Senn2013, sm="MD", comb.random=TRUE)
forest(net2, xlim=c(-1.5,1), ref="plac",
xlab="HbA1c difference",
leftcols="studlab",
leftlabs="Contrast to placebo")
#
# Add column with P-Scores on right side of forest plot
#
forest(net2, xlim=c(-1.5,1), ref="plac",
xlab="HbA1c difference",
leftcols="studlab",
leftlabs="Contrast to placebo",
rightcols=c("effect", "ci", "Pscore"),
rightlabs="P-Score",
just.addcols="right")
#
# Add column with P-Scores on left side of forest plot
#
forest(net2, xlim=c(-1.5,1), ref="plac",
xlab="HbA1c difference",
leftcols=c("studlab", "Pscore"),
leftlabs=c("Contrast to placebo", "P-Score"),
just.addcols="right")
#
# Sort forest plot by descending P-Score
#
forest(net2, xlim=c(-1.5,1), ref="plac",
xlab="HbA1c difference",
leftcols="studlab",
leftlabs="Contrast to placebo",
rightcols=c("effect", "ci", "Pscore"),
rightlabs="P-Score",
just.addcols="right",
sortvar=-Pscore)
Run the code above in your browser using DataLab