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

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

728 Laboratory assays<br />

Proportion positive, P<br />

1<br />

0<br />

Distribution of tolerances<br />

Dose, X<br />

Fig. 20.4 Quantal-response curve with the correspond<strong>in</strong>g tolerance distribution.<br />

In an assay of T and S, the response curves, plotted aga<strong>in</strong>st x ˆ log X,willbe<br />

parallel <strong>in</strong> the sense that they differ by a constant horizontal distance, log r, and<br />

a full analysis requires that the data be fitted by two parallel curves of this sort. A<br />

natural approach is to l<strong>in</strong>earize the response curves by apply<strong>in</strong>g to the response<br />

one of the standard transformations for b<strong>in</strong>ary data (§§10.8 and 14.1). One can<br />

then suppose that the expectations of the transformed values are l<strong>in</strong>early related<br />

to the log dose, x, by the equations<br />

)<br />

Y ˆ aS ‡ bx for S<br />

and<br />

: …20:20†<br />

Y ˆ aT ‡ bx for T<br />

The slope parameter b is <strong>in</strong>versely proportional to the standard deviation of<br />

the log-tolerance distribution. For the logistic transformation (14.5) with<br />

Y ˆ ln‰P=…1 P†Š, for <strong>in</strong>stance, the proportionality factor is p= 3<br />

p ˆ 1 814.<br />

Although (20.20) has essentially the same form as (20.3) and (20.4), the<br />

methods described <strong>in</strong> §20.1 for parallel-l<strong>in</strong>e assays cannot be used, s<strong>in</strong>ce the<br />

random variation about the l<strong>in</strong>es is not normally distributed with constant<br />

variance. The model may be fitted iteratively by the maximum likelihood procedure<br />

for logistic regression described <strong>in</strong> §14.2, us<strong>in</strong>g a dummy variable to<br />

dist<strong>in</strong>guish between S and T, as <strong>in</strong> (11.55). Writ<strong>in</strong>g aT ˆ aS ‡ d, the log potency<br />

ratio is given by log r ˆ d=b. A computer program for logistic regression may be<br />

used to obta<strong>in</strong> the estimates aS, d and b of aS, d and b, respectively, and hence log<br />

r is estimated by M ˆ d=b. From (5.17), an approximate formula for var(M) is<br />

var…M† ' Vd<br />

b 2<br />

2dCdb<br />

b 3<br />

‡ d2 Vb<br />

b 4 ;

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

Saved successfully!

Ooh no, something went wrong!