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.

170 Bayesian methods<br />

<strong>in</strong> this way. For example, <strong>in</strong> the present problem, the user might believe that the<br />

mean lies <strong>in</strong> the neighbourhood of either of two values, u0 or u1. It might then be<br />

appropriate to use a bimodal prior distribution with peaks at these two values. In<br />

that case, the simplicity afforded by the conjugate family would be lost, and the<br />

posterior distribution would no longer take the normal form (6.1).<br />

Secondly, if either n is very large (when the evidence from the data overwhelms<br />

the prior <strong>in</strong>formation) or if s2 0 is very large (when the prior evidence is<br />

very weak and the prior distribution is non-<strong>in</strong>formative), the posterior distribution<br />

(6.1) tends towards the likelihood N(x, s2 =n).<br />

In pr<strong>in</strong>ciple, once the formulations for the prior and likelihood have been<br />

accepted as appropriate, the posterior distribution provides all we need for<br />

<strong>in</strong>ference about m. In practice, as <strong>in</strong> frequentist <strong>in</strong>ference, it will be useful to<br />

consider ways of answer<strong>in</strong>g specific questions about the possible value of m. In<br />

particular, what are the Bayesian analogues of the two pr<strong>in</strong>cipal modes of<br />

<strong>in</strong>ference discussed <strong>in</strong> §4.1: significance tests and confidence <strong>in</strong>tervals?<br />

Bayesian significance tests<br />

Suppose that, <strong>in</strong> the formulation lead<strong>in</strong>g up to (6.1), we wanted to ask whether<br />

there was strong evidence that m < 0orm > 0. In frequentist <strong>in</strong>ference we should<br />

test the hypothesis that m ˆ 0, and see whether it was strongly contradicted by a<br />

significant result <strong>in</strong> either direction. In the present Bayesian formulation there is<br />

no po<strong>in</strong>t <strong>in</strong> consider<strong>in</strong>g the probability that m is exactly 0, s<strong>in</strong>ce that probability<br />

is zero (although m ˆ 0 has a non-zero density). However, we can state directly<br />

the probability that, say m < 0 by calculat<strong>in</strong>g the tail area to the left of zero <strong>in</strong> the<br />

normal distribution (6.1).<br />

It is <strong>in</strong>structive to note what happens <strong>in</strong> the limit<strong>in</strong>g case considered above,<br />

when the sample size is large or the prior is non-<strong>in</strong>formative and the posterior<br />

distribution is N(x, s2 =n). The posterior probability that m < 0 is the probability<br />

of a standardized normal deviate less than<br />

0 x x n<br />

p ˆ<br />

s= n<br />

p<br />

s ,<br />

and this is precisely the same as the one-sided P value obta<strong>in</strong>ed <strong>in</strong> a frequentist<br />

test of the null hypothesis that m ˆ 0. The posterior tail area and the one-sided P<br />

value are thus numerically the same, although of course their strict <strong>in</strong>terpretations<br />

are quite different.<br />

Example 6.1<br />

Example 4.1 described a frequentist significance test based on a sample of n ˆ 100 survival<br />

times of patients with a form of cancer. The observed mean was x ˆ 46 9 months, and the

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

Saved successfully!

Ooh no, something went wrong!