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.

16.2 Prior and posterior distributions 537<br />

improper prior is used, there is no guarantee that a proper posterior distribution<br />

will result.<br />

Another way to model ignorance is to use a uniform distribution over a f<strong>in</strong>ite<br />

range. The f<strong>in</strong>ite range means that the value of p…u† can be chosen so that<br />

the distribution <strong>in</strong>tegrates to one. The improper prior, mentioned above, can<br />

be obta<strong>in</strong>ed as the limit as the range tends to <strong>in</strong>f<strong>in</strong>ity. The uniformity of the<br />

distribution means that no one value is preferred over any other, which might<br />

seem as good an implementation of ignorance as one might expect. However, the<br />

use of a uniform distribution to model ignorance is not as satisfactory as it might<br />

at first appear. Suppose u is a scalar. If there is ignorance about the value of u,<br />

then surely there is equal ignorance about any function of u, such as u2 or<br />

exp( u). However, a uniform distribution for u implies a non-uniform distribution<br />

for these (and most other) functions of u.<br />

A similar problem applies when the range of the parameter is <strong>in</strong>f<strong>in</strong>ite. The<br />

discussion of improper priors given above tacitly assumed that the range of the<br />

parameters is the whole real l<strong>in</strong>e. However, for many parameters, this is not true;<br />

for example, a standard deviation, s, is necessarily positive. However, log s can<br />

take any value and an improper prior for this quantity that is constant on the<br />

whole real l<strong>in</strong>e transforms to a prior p…s† ˆ1=s for the standard deviation. This<br />

form of improper prior is widely used for parameters that are necessarily<br />

positive.<br />

A way to def<strong>in</strong>e prior distributions that express ignorance and are <strong>in</strong>variant<br />

under transformations of the parameters was proposed by Jeffreys (1961, p. 181).<br />

For scalar u the Jeffreys prior is taken to be proportional to the square root of<br />

the Fisher <strong>in</strong>formation, i.e.<br />

…<br />

q<br />

p…u† /<br />

2<br />

" ( ) # 1=2<br />

2 log p…xju† p…xju† dx : …16:3†<br />

qu<br />

When u is a vector, the Jeffreys prior is, strictly speak<strong>in</strong>g, proportional to the<br />

square root of the determ<strong>in</strong>ant of the matrix whose (i, j)th element is<br />

…<br />

q 2<br />

( )<br />

log p…xju† p…xju† dx:<br />

quiquj<br />

The expression for vector u is often cumbersome and it is quite common to form<br />

the Jeffreys prior separately for each component of u us<strong>in</strong>g (16.3) and then use<br />

the product of these as the prior for u. The argument here is that, broadly<br />

speak<strong>in</strong>g, ignorance is consistent with <strong>in</strong>dependence, so form<strong>in</strong>g the jo<strong>in</strong>t prior<br />

<strong>in</strong> this way is reasonable.<br />

Jeffreys priors can lead to some slightly unusual results. If the data have a<br />

univariate normal distribution, with mean m and standard deviation s, then the<br />

follow<strong>in</strong>g can be obta<strong>in</strong>ed:

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

Saved successfully!

Ooh no, something went wrong!