data("OW_below")
#Plot of positive likelihood ratios for inadequate GWG for overweight women
plot(OW_below$month, OW_below$lr.pos.est,
xlim=c(1,10),
ylim=c(0,20),
main="LR+ for inadequate GWG",
ylab="likelihood ratio",
xlab="month",
cex.main=2,
font.main=1,
cex.lab=1.9,
cex.axis=1.5)
lines(OW_below$month, OW_below$lr.pos.est, lty=1,lwd=4)
lines(OW_below$month, OW_below$lr.pos.lower, lty=2,lwd=4)
lines(OW_below$month, OW_below$lr.pos.upper, lty=2,lwd=4)
abline(h=2,lwd=1, col="grey")
abline(h=4,lwd=1, col="grey")
abline(h=6,lwd=1, col="grey")
abline(h=8,lwd=1, col="grey")
abline(h=10,lwd=1, col="grey")
abline(h=12,lwd=1, col="grey")
abline(h=14,lwd=1, col="grey")
abline(h=16,lwd=1, col="grey")
abline(h=18,lwd=1, col="grey")
axis(1, c(3,5,7,9), cex.axis=1.5)
legend("topleft", c("lr+ estimate", "lr+ confidence interval"),
col=c("black", "black"), lty=c(1,2),
bg="white", cex=1.5, lwd=c(3,3))
Run the code above in your browser using DataLab