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.

526 Empirical methods for categorical data<br />

the marg<strong>in</strong>al totals, the test <strong>in</strong>volves evaluation of the probability that P rixi is<br />

at least as large as the observed value (Agresti, 1990, §4.8.2). The exact method<br />

<strong>in</strong>volves a large amount of calculation and is only feasible with appropriate<br />

statistical software, such as StatXact.<br />

The exact test may be extended to cover the test<strong>in</strong>g of a trend when there is a<br />

stratify<strong>in</strong>g variable to take <strong>in</strong>to account (§15.7). The x 2 test for this situation is<br />

(15.21). An exact test is based on the probability of obta<strong>in</strong><strong>in</strong>g a value of P<br />

h Th at<br />

least as large as that observed and, from (15.4), this is equivalent to bas<strong>in</strong>g the<br />

test on the sum of P rixi over the h strata. Condition<strong>in</strong>g is on the marg<strong>in</strong>al totals<br />

<strong>in</strong> all the strata.<br />

Example 15.9<br />

In Table 15.8 there are several low counts and so there might be some doubt as to the<br />

accuracy of the analysis shown <strong>in</strong> Example 15.8. Us<strong>in</strong>g StatXact the test of trend is as<br />

follows.<br />

The sum of P rixi over the three strata is 35. The exact probability that this sum would<br />

be greater than or equal to 35 if there were no effect of length of exposure, and subject to<br />

the marg<strong>in</strong>al totals <strong>in</strong> each of the three age groups, is 0 0050 from StatXact. The<br />

correspond<strong>in</strong>g one-sided probability us<strong>in</strong>g the x2 …1† of 7 42 <strong>in</strong> Example 15.8 is 0 0032.<br />

However, this is based on a test without a cont<strong>in</strong>uity correction and is, therefore,<br />

equivalent to a mid-P test (see §4.4), whereas the exact value is the usual P value.<br />

If the test statistic <strong>in</strong> Example 15.8 had been corrected for cont<strong>in</strong>uity, then its value<br />

would have been 7 242 =8 45 ˆ 6 52, giv<strong>in</strong>g a one-sided P value of 0 0054very similar to<br />

the exact value of 0 0050.<br />

Conversely, s<strong>in</strong>ce from StatXact the probability that P rixi over the three strata is<br />

1<br />

exactly 35 is 0 0033, the exact one-sided mid-P value is 0 0050 2 …0 0033† ˆ0 0034, very<br />

near to the value from the uncorrected x2 of 0 0032.<br />

The tests have been compared on the basis of their one-sided values s<strong>in</strong>ce, with the<br />

exact test, there can be ambiguity on how to obta<strong>in</strong> a two-sided value (see discussion <strong>in</strong><br />

Example 4.13). The two-sided levels us<strong>in</strong>g the x2 are, of course, simply double the onesided<br />

levels, giv<strong>in</strong>g 0 0065 and 0 011 for the uncorrected and corrected values. One option<br />

of obta<strong>in</strong><strong>in</strong>g a two-sided exact level, which we advocated, is to double the one-sided level,<br />

giv<strong>in</strong>g 0 0067 and 0 010 for the mid-P and P values, respectively.<br />

So, <strong>in</strong> this example, even though the frequencies are not large, the tests based on x2 statistics proved very satisfactory approximations to the exact tests.<br />

Comb<strong>in</strong>ation of 2 2 tables<br />

This is an extension of the exact test <strong>in</strong> a s<strong>in</strong>gle 2 2 table (§4.5). The one-tailed<br />

significance level is the probability that P rAi is equal to or greater than its<br />

observed value, where for any value of P rAi the probability is calculated by<br />

consider<strong>in</strong>g all the possible comb<strong>in</strong>ations of tables, with the same marg<strong>in</strong>al<br />

totals, over the strata that produce this total (Mehta et al., 1985; Hirji et al.,<br />

1988; Agresti, 1990, §7.4.4).

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

Saved successfully!

Ooh no, something went wrong!