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.

368 Modell<strong>in</strong>g cont<strong>in</strong>uous data<br />

Table 11.7 Residual diagnostic statistics for Fig. 11.10.<br />

Fig. 11.10 Fitted slope s r10 t10 Leverage h10 Cook's D10<br />

(a) 1 00 1 00<br />

(b) 1 03 0 94 0 14 0 13 0 54 0 01<br />

(c) 0 64 1 11 1 51 1 68 0 54 1 37<br />

(d) 1 61 1 37 2 07 2 84 0 54 2 55<br />

(e) 1 21 2 15 2 55 5 49 0 11 0 41<br />

and 7 DF, respectively, the studentized and jackknife residuals are clearly not<br />

remarkable for either Fig. 11.10(b) or (c). The greater sensitivity of the<br />

jackknife residual is shown for Fig. 11.10(d and e); the lower values of the<br />

studentized residual, r10, <strong>in</strong> these two cases is a consequence of the <strong>in</strong>crease <strong>in</strong><br />

standard deviation due to the outlier under test. For Fig. 11.10(d and e), the<br />

extra po<strong>in</strong>t is suggested as a possible outlier by the values of the jackknife<br />

residual. We now consider the formulation of a significance test of the residuals.<br />

For Fig. 11.10(d), t10 ˆ 2 84 and if this is assessed as a t statistic with 7 DF<br />

then P ˆ 0 025. However, this takes no account of the fact that this residual has<br />

been chosen as the largest, <strong>in</strong> absolute magnitude, and the effect of this is quite<br />

considerable and <strong>in</strong>creases with <strong>in</strong>creas<strong>in</strong>g n. The required significance level is<br />

then the probability that the largest residual <strong>in</strong> a set of n residuals will be as large<br />

as or larger than the observed value. This probability is the probability that not<br />

all 10 residuals are <strong>in</strong> the central part of the distribution and is given approximately<br />

by P ˆ 1 …1 0 025† 10 ˆ 0 22, s<strong>in</strong>ce the jackknife residuals are almost<br />

<strong>in</strong>dependent. This is the Bonferroni correction (see §13.3). An equivalent means<br />

of assessment is to work with a lower critical significance probability, and the<br />

largest residual would be significant at level a only if the t statistic gave P < a=n.<br />

For Fig. 11.10(d), and a significance test at the 0 05 level, we would require<br />

P < 0 005 so that the largest residual with a nom<strong>in</strong>al P of 0 025 is not significant.<br />

In Fig. 11.10(e) the largest residual has a nom<strong>in</strong>al P of just less than 0 001 and so<br />

there is significant evidence …P < 0 01…ˆ 10 0 001†† that the 10th po<strong>in</strong>t is an<br />

outlier. Kle<strong>in</strong>baum et al. (1998) give a table (Table A±8A) to facilitate the<br />

assessment of the significance of jackknife residuals and from this table the<br />

significance level is slightly less than 0 01 (critical value 5 41).<br />

Check<strong>in</strong>g for <strong>in</strong>fluential po<strong>in</strong>ts<br />

A po<strong>in</strong>t may be regarded as <strong>in</strong>fluential if it exerts more than its fair share <strong>in</strong><br />

determ<strong>in</strong><strong>in</strong>g the values of the regression coefficients. Conceptually the <strong>in</strong>fluence<br />

of the ith po<strong>in</strong>t may be envisaged <strong>in</strong> terms of the changes <strong>in</strong> the estimates of the

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

Saved successfully!

Ooh no, something went wrong!