CoalMiners                package:vcd                R Documentation

_B_r_e_a_t_h_l_e_s_s_n_e_s_s _a_n_d _W_h_e_e_z_e _i_n _C_o_a_l _M_i_n_e_r_s

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

     Data from Ashford & Snowden (1970) given by Agresti (1990) on the
     association between two pulmonary conditions, breathlessness and
     wheeze, in a large sample of coal miners.

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

     data(CoalMiners)

_F_o_r_m_a_t:

     A 3-dimensional array resulting from cross-tabulating variables
     for 16,330 coal miners.  The variables and their levels are as
     follows:

       No  Name            Levels
        1  Wheeze          W, NoW
        2  Breathlessness  B, NoB
        3  Age             25-29, 30-34, ..., 60-64

_S_o_u_r_c_e:

     Michael Friendly (2000), Visualizing Categorical Data, pages
     82-83, 319-322.

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

     A. Agresti (1990), _Categorical Data Analysis_.
     Wiley-Interscience, New York.

     J. R. Ashford & R. D. Snowdon (1970), Multivariate probit
     analysis, _Biometrics_, *26*, 535-546.

     M. Friendly (2000), _Visualizing Categorical Data_. SAS Institute,
     Cary, NC.

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

     data(CoalMiners)

     ## Fourfold display, both margins equated
     fourfold(CoalMiners, mfcol = c(2,4))

     ## Log Odds Ratio Plot
     summary(l <- oddsratio(CoalMiners))
     g <- seq(25, 60, by = 5)
     plot(l,
          xlab = "Age Group",
          main = "Breathlessness and Wheeze in Coal Miners")
     m <- lm(l ~ g + I(g^2))
     lines(fitted(m), col = "red")

     ## Fourfold display, strata equated
     fourfold(CoalMiners, std = "ind.max", mfcol = c(2,4))

