next up previous contents
Next: Fitting Alternating Logistic Regressions Up: Worked Examples Previous: Modelling over-dispersed Poisson counts

Modelling binary longitudinal data with ALR: the cerebro data

  As a final example, we turn to the cerebrovascular deficiency crossover trial discussed in Diggle et al. (1994), Section 8.3. In this trial, 67 patients were given an active drug (A) and a placebo (B) in the order A,B (the AB group) for 34 of the patients (allocated at random), and in the order B,A (the BA group) for the remaining 33. For each patient and each treatment there is a binary response: 0 for an abnormal ECG reading, and 1 for a normal reading. The data are available in the cerebro balanced data structure, and tabulated (after first combining the period 1 and 2 responses) below:


\begin{Example}
\gt combinations <- apply(tsy(cerebro),1,paste,collapse='','')
\...
 ...nations, groups(cerebro))
 AB BA 
0,0 6 9
0,1 0 4
1,0 6 2
1,1 22 18\end{Example}

We wish to determine whether the active drug increases the probability of a normal ECG reading.



 

David Smith
4/2/1998