tomogplot {MCMCpack}R Documentation

Tomography Plot

Description

tomogplot is used to produce a tomography plot (see King, 1997) for a series of partially observed 2 x 2 contingency tables.

Usage

tomogplot(r0, r1, c0, c1, xlab="fraction of r0 in c0 (p0)",
          ylab="fraction of r1 in c0 (p1)", bgcol="white", ...)
   

Arguments

r0 An (ntables * 1) vector of row sums from row 0.
r1 An (ntables * 1) vector of row sums from row 1.
c0 An (ntables * 1) vector of column sums from column 0.
c1 An (ntables * 1) vector of column sums from column 1.
xlab The x axis label for the plot.
ylab The y axis label for the plot.
bgcol The background color for the plot.
... further arguments to be passed

Details

Consider the following partially observed 2 by 2 contingency table:

| Y=0 | Y=1 |
- - - - - - - - - - - - - - - - - - - -
X=0 | Y0 | | r0
- - - - - - - - - - - - - - - - - - - -
X=1 | Y1 | | r1
- - - - - - - - - - - - - - - - - - - -
| c0 | c1 | N

where r0, r1, c0, c1, and N are non-negative integers that are observed. The interior cell entries are not observed. It is assumed that Y0|r0 ~ Binomial(r0, p0) and Y1|r1 ~ Binomial(r1,p1).

This function plots the bounds on the maximum likelihood estimatess for (p0, p1).

References

Gary King, 1997. A Solution to the Ecological Inference Problem. Princeton: Princeton University Press.

Jonathan Wakefield. 2001. ``Ecological Inference for 2 x 2 Tables,'' Center for Statistics and the Social Sciences Working Paper no. 12. University of Washington.

See Also

MCMChierEI, MCMCdynamicEI, dtomogplot

Examples

r0 <- rpois(100, 500) 
r1 <- rpois(100, 200) 
c0 <- rpois(100, 100) 
c1 <- (r0 + r1) - c0 
tomogplot(r0, r1, c0, c1) 

[Package MCMCpack version 0.5-2 Index]