01.06.2013 Views

Statistical Methods in Medical Research 4ed

Statistical Methods in Medical Research 4ed

Statistical Methods in Medical Research 4ed

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

drop-out. In this it is supposed that complete data are observed on an <strong>in</strong>dividual<br />

up to and <strong>in</strong>clud<strong>in</strong>g measurement occasion D 1, where 2 D k ‡ 1, with all<br />

subsequent outcomes be<strong>in</strong>g miss<strong>in</strong>g: D ˆ k ‡ 1 <strong>in</strong>dicates that the <strong>in</strong>dividual did<br />

not drop out. This reduces the modell<strong>in</strong>g problem, as the miss<strong>in</strong>g value process<br />

has been greatly restricted. While this undoubtedly limits the practical utility of<br />

the method, it is an important case and also allows progress <strong>in</strong> a methodologically<br />

challeng<strong>in</strong>g area.<br />

Diggle and Kenward (1994) and Diggle (1998) suggest that the miss<strong>in</strong>g value<br />

process, which can now also be called a drop-out process, can be modelled by<br />

log<br />

Pr…D ˆ d jD d, y†<br />

1 Pr…D ˆ d jD d, y† ˆ ad ‡ Ps<br />

12.6 Longitud<strong>in</strong>al data 447<br />

jˆ0<br />

b jyd j, …12:57†<br />

where s is a prespecified <strong>in</strong>teger and for notational convenience we take yj ˆ 0<br />

for j 0. This model is not the only one that could be considered, but it is<br />

plausible; it assumes that the probability of drop-out at any given stage can<br />

depend on the outcomes up to and <strong>in</strong>clud<strong>in</strong>g that stage, but not on the future of<br />

the process. If all the bs are 0 then the process is MCAR, and if b 0 ˆ 0 the process<br />

is MAR. If neither of these simplifications obta<strong>in</strong>s then the log-likelihood for the<br />

observed data has an extra term over those <strong>in</strong> (12.56), which accommodates the<br />

way the parameters b and u jo<strong>in</strong>tly determ<strong>in</strong>e aspects of the data. This will<br />

<strong>in</strong>volve contributions from probabilities <strong>in</strong> (12.57) but this depends on the<br />

unobserved value yd, so it is necessary to compute an expression for<br />

P(D ˆ d jD d, y o) by <strong>in</strong>tegrat<strong>in</strong>g the probability <strong>in</strong> (12.57) over the conditional<br />

distribution of the unobserved value, f …yd jy o†. By this means valid <strong>in</strong>ferences<br />

can be drawn about u even when data are miss<strong>in</strong>g not at random.<br />

The MCAR, MAR classification of Little and Rub<strong>in</strong> (1987) is naturally<br />

expressed by factoriz<strong>in</strong>g f …y, r† as f …rjy†f …y† and then specify<strong>in</strong>g particular<br />

forms for the first factor; these are known as selection models. Another approach<br />

would be to use the alternative factorization, f …yjr†f …r†, which gives rise to<br />

pattern mixture models, so called because they view the jo<strong>in</strong>t distribution as a<br />

mixture of different distributions on the outcomes, one distribution for each<br />

pattern of miss<strong>in</strong>g data. This view, which is extensively discussed <strong>in</strong> Little (1993),<br />

Hogan and Laird (1997), Diggle (1998) and Kenward and Molenberghs (1999),<br />

gives an alternative and <strong>in</strong> many ways complementary approach to the problem.<br />

The reliance of the MAR, MCAR classification on the form of f …rjy† means that<br />

the extent to which this classification can be carried over <strong>in</strong>to pattern mixture<br />

models is not transparent. Molenberghs et al. (1998) have shown that for dropout,<br />

but not for more general patterns of miss<strong>in</strong>g data, MAR models can be seen<br />

to correspond to certa<strong>in</strong> classes of pattern mixture model.<br />

In order to make valid <strong>in</strong>ferences about f …y, r† from <strong>in</strong>complete data …y o, r†,<br />

various models have to be constructed or assumptions made. A fundamental

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!