23.10.2012 Views

View PDF Version - RePub - Erasmus Universiteit Rotterdam

View PDF Version - RePub - Erasmus Universiteit Rotterdam

View PDF Version - RePub - Erasmus Universiteit Rotterdam

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.

For the indirect approach we rewrite Pr(R =1| Y = y) as<br />

Pr(R =1| Y = y) =<br />

Dynamic prediction of response using longitudinal profi les 83<br />

pSE<br />

pSE+(1− p)(1− SP)<br />

using Bayes’ theorem. Here p = Pr(R = 1) is the overall fraction of respon-<br />

ders (the ’prior probability’) and Pr(Y =1| R = 1) is the behavior of Y in the<br />

group of responders, i.e. the sensitivity (SE) of Y and Pr(Y =1| R =0)isthe<br />

behavior of Y in the group of non-responders, corresponding to 1-specificity (1-<br />

SP). Bayes’ theorem then states that the left hand side is the ’posterior probability’<br />

of response given Y = y. In this simple situation the direct and the<br />

indirect methods are similar and α = log ((1 − SE)/SP ) + log (p/(1 − p)) and<br />

β = log (SE · SP/((1 − SE)(1 − SP))).<br />

Suppose Y is continuous, then the indirect approach is written<br />

Pr(R =1| Y = y) =<br />

p f1(y)<br />

p f1(y)+(1− p) f0(y)<br />

where f1(y) represents the probability density function of Y at the point y in the<br />

presence of response (R=1), and f0(y) the probability density function of Y at y in<br />

the absence of response (R=0), illustrated in Figure 1a. The direct method and<br />

the indirect method are the same when Y is normal distributed in both populations<br />

of responders and non-responders with mean μ1 and μ0, respectively, and common<br />

variance σ2 . In this case α = − μ2 1−μ2 0<br />

2σ2 and β = μ1−μ0<br />

σ2 . Both methods are easily<br />

generalized to handle a multivariate normal vector Y of covariates with common<br />

covariance matrix for the responders and non-responders. For the indirect approach<br />

the requirement of a common covariance matrix for responders and non-responders<br />

is ignored in the longitudinal prediction.<br />

In the following the two approaches are extended to handle repeated measurements<br />

of Y over time.<br />

Longitudinal prediction with the direct approach<br />

The observed makers as predictors<br />

Suppose a baseline prediction at t=0 exist and at visit j the markers Y i,j are observed.<br />

A simple update of the baseline prediction model is to extend the prediction at t=0,<br />

pi,0 with the new information observed at visit j, as suggested by Steyerberg(chapter<br />

20), 11 i.e.

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

Saved successfully!

Ooh no, something went wrong!