pvpORMH: Drawing p-value function plot by stratified or pooled cross...

In fmsb: Functions for Medical Statistics Book with some Demographic Data

Drawing p-value function plot by stratified or pooled cross tables

Description

Drawing the p-value function (a.k.a. nested confidence intervals) plot of pooled odds ratios (pORs) for several 2 by 2 crosstables, which are stratified by a confounding variable or pooled for several studies, with Mantel-Haenszel's method.

Usage

 pvpORMH(XTAB, xrange, add, ...)

Arguments

XTAB	A matrix with 4 columns. The first column is the number of exposed cases. The second column is the number of unexposed cases. The third column is the number of exposed controls. The forth column is the number of unexposed controls. Rows should be composed of different strata or studies.
xrange	A numeric vector includes 2 elements for minimum and maximum of x axis. Default is c(0.6, 1.2).
add	Logical. If TRUE, the line is overlayed on the existing pvalueplot, otherwise the graph is newly plotted. Default FALSE.
...	Other options handed down to plot() or lines(). pch, lty or col may be useful.

Value

A data.frame containing the set of p-values (ranging from 0.0005 to 1) and corresponding pORs are returned.

Author(s)

Minato Nakazawa [e-mail] https://minato.sip21c.org/

References

Rothman KJ (2012) Epidemiology: An Introduction. 2nd Ed., Oxford University Press, Oxford.

Examples

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)

fmsb documentation built on May 29, 2024, 2:02 a.m.