tomogplot              package:MCMCpack              R Documentation

_T_o_m_o_g_r_a_p_h_y _P_l_o_t

_D_e_s_c_r_i_p_t_i_o_n:

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

_U_s_a_g_e:

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

_A_r_g_u_m_e_n_t_s:

      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

_D_e_t_a_i_l_s:

     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).

_R_e_f_e_r_e_n_c_e_s:

     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.

_S_e_e _A_l_s_o:

     'MCMChierEI', 'MCMCdynamicEI', 'dtomogplot'

_E_x_a_m_p_l_e_s:

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

