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]