TenStudies <- matrix(
c(215, 229, 311-215, 306-229,
38, 33, 59-38, 51-33,
161, 174, 293-161, 293-174,
76, 88, 164-76, 163-88,
103, 105, 129-103, 133-105,
65, 67, 120-65, 125-67,
81, 75, 113-81, 110-75,
48, 63, 160-48, 159-63,
22, 21, 60-22, 62-21,
56, 51, 137-56, 140-51
), 10, 4, byrow=TRUE)
ElevenStudies <- rbind(TenStudies, c(468, 480, 229, 205))
# Figure 8-4 in Chapter 8 of Rothman's textbook.
pvpORMH(TenStudies)
pvpORMH(ElevenStudies, add=TRUE, lty=2)
segments(1, 0, 1, 1, lwd=2)
Run the code above in your browser using DataLab