MCMCpack: tomogplot – R documentation

Become an expert in R — Interactive courses, Cheat Sheets, certificates and more!

tomogplot

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 \times 1) vector of row sums from row 0.
`r1`	An (ntables \times 1) vector of row sums from row 1.
`c0`	An (ntables \times 1) vector of column sums from column 0.
`c1`	An (ntables \times 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	\| Y_0	\|	\| r_0
---------	---------	---------	---------
X=1	\| Y_1	\|	\| r_1
---------	---------	---------	---------
	\| c_0	\| c_1	\| N

where r_0, r_1, c_0, c_1, and N are non-negative integers that are observed. The interior cell entries are not observed. It is assumed that Y_0|r_0 \sim \mathcal{B}inomial(r_0, p_0) and Y_1|r_1 \sim \mathcal{B}inomial(r_1, p_1).

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 C. Wakefield. 2004. “Ecological Inference for 2 x 2 Tables.” Journal of the Royal Statistical Society, Series A. 167(3): 385445.

Examples

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

MCMCpack

Markov Chain Monte Carlo (MCMC) Package

v1.5-0

GPL-3

Authors

Andrew D. Martin [aut], Kevin M. Quinn [aut], Jong Hee Park [aut,cre], Ghislain Vieilledent [ctb], Michael Malecki[ctb], Matthew Blackwell [ctb], Keith Poole [ctb], Craig Reed [ctb], Ben Goodrich [ctb], Ross Ihaka [cph], The R Development Core Team [cph], The R Foundation [cph], Pierre L'Ecuyer [cph], Makoto Matsumoto [cph], Takuji Nishimura [cph]

Initial release

2021-01-19