bacteria                package:MASS                R Documentation

_P_r_e_s_e_n_c_e _o_f _B_a_c_t_e_r_i_a _a_f_t_e_r _D_r_u_g _T_r_e_a_t_m_e_n_t_s

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

     Tests of the presence of the bacteria H. influenzae in children
     with otitis media in the Northern Territory of Australia.

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

     data(bacteria)

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

     This data frame has 220 rows and the following columns:

     _y presence or absence: a factor with levels `n' and `y'.

     _a_p active/placebo: a factor with levels `a' and `p'.

     _h_i_l_o hi/low compliance: a factor with levels `hi' amd `lo'.

     _w_e_e_k numeric: week of test.

     _I_D subject ID: a factor.

     _t_r_t a factor with levels `placebo', `drug' and `drug+', a
          re-coding of `ap' and `hilo'.

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

     Dr A. Leach tested the effects of a drug on 50 children with a
     history of otitis media in the Northern Territory of Australia. 
     The children were randomized to the drug or the a placebo, and
     also to receive active encouragement to comply with taking the
     drug.

     The presence of H. influenzae was checked at weeks 0, 2, 4, 6 and
     11: 30 of the checks were missing and are not included in this
     data frame.

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

     Menzies School of Health Research 1999-2000 Annual Report pp.
     18-21 (<URL:
     http://www.menzies.edu.au/publications/anreps/MSHR00.pdf>).

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

     data(bacteria)
     contrasts(bacteria$trt) <- structure(contr.sdif(3),
          dimnames = list(NULL, c("drug", "encourage")))
     ## fixed effects analyses
     summary(glm(y ~ trt * week, binomial, data = bacteria))
     summary(glm(y ~ trt + week, binomial, data = bacteria))
     summary(glm(y ~ trt + I(week > 2), binomial, data = bacteria))

     # conditional random-effects analysis
     library(survival)
     bacteria$Time <- rep(1, nrow(bacteria))
     coxph(Surv(Time, unclass(y)) ~ week + strata(ID),
           data = bacteria, method = "exact")
     coxph(Surv(Time, unclass(y)) ~ factor(week) + strata(ID),
           data = bacteria, method = "exact")
     coxph(Surv(Time, unclass(y)) ~ I(week > 2) + strata(ID),
           data = bacteria, method = "exact")

     # PQL glmm analysis
     library(nlme)
     summary(glmmPQL(y ~ trt + I(week > 2), random = ~ 1 | ID,
                     family = binomial, data = bacteria))

