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.

E…y† ˆa ‡ bx), or whether, as <strong>in</strong> (c), the l<strong>in</strong>es co<strong>in</strong>cide. In practice, the fitted<br />

regression l<strong>in</strong>es would rarely have precisely the same slope or position, and the<br />

question is to what extent differences between the l<strong>in</strong>es can be attributed to<br />

random variation. Differences <strong>in</strong> position between parallel l<strong>in</strong>es are discussed <strong>in</strong><br />

§11.5. In this section we concentrate on the question of differences between<br />

slopes.<br />

Suppose that there are k groups, with ni pairs of observations <strong>in</strong> the ith<br />

group. Denote the mean values of x and y <strong>in</strong> the ith group by xi and yi, and the<br />

regression l<strong>in</strong>e, calculated as <strong>in</strong> §7.2, by<br />

Yi ˆ y ‡ bi…x xi†: …11:12†<br />

If all the ni are reasonably large, a satisfactory approach is to estimate the<br />

variance of each bi by (7.16) and to ignore the imprecision <strong>in</strong> these estimates of<br />

variance. Chang<strong>in</strong>g the notation of (7.16) somewhat, we shall denote the residual<br />

mean square for the ith group by s 2 i and the sum of squares of x about xi by<br />

P<br />

…i† …x xi† 2 . Note that the parenthesized suffix i attached to the summation<br />

sign <strong>in</strong>dicates summation only over the specified group i; that is,<br />

P<br />

…x xi†<br />

…i†<br />

2 ˆ Pni<br />

…xij xi†<br />

jˆ1<br />

2 :<br />

To simplify the notation, denote this sum of squares about the mean of x <strong>in</strong><br />

the ith group by<br />

Sxxi, …11:13†<br />

the sum of products of deviations by Sxyi, and so on. Then follow<strong>in</strong>g the method<br />

of §8.2, we write<br />

and calculate<br />

wi ˆ 1<br />

var…bi†<br />

G ˆ P wib 2 i<br />

Sxxi<br />

ˆ<br />

s2 , …11:14†<br />

i<br />

P<br />

… wibi†<br />

2 = P wi: …11:15†<br />

On the null hypothesis that the true slopes bi are all equal, G follows approximately<br />

a x2 …k 1† distribution. High values of G <strong>in</strong>dicate departures from the null<br />

hypothesis, i.e. real differences between the bi.IfG is non-significant, and the<br />

null hypothesis is tentatively accepted, the common value b of the bi is best<br />

estimated by the weighted mean<br />

with an estimated variance<br />

11.4 Regression <strong>in</strong> groups 323<br />

b ˆ P wibi= P wi, …11:16†

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

Saved successfully!

Ooh no, something went wrong!