next up previous contents
Next: Prior Elicitation and Non-informative Up: Bayesian Approach Previous: Kernel Density Estimation

A Simple Example - tex2html_wrap_inline2155

Consider the case of drawing from a population of unknown mean, tex2html_wrap_inline2157 , but known variance tex2html_wrap_inline2159 . ( tex2html_wrap_inline1949 is termed precision, and is just the reciprocal of variance.)

The model is that the data, X, will be normally distributed with unknown mean but given variance. Thus, in terms of a single observation, x, we can write down the likelihood;

displaymath2145

The next step is to elicit a prior for tex2html_wrap_inline2157 . It may be reasonable to assume that the prior beliefs about tex2html_wrap_inline2157 can be expressed as a normal distribution, that is

displaymath2146

where both tex2html_wrap_inline2171 and tex2html_wrap_inline2173 are specified. Typically tex2html_wrap_inline2171 is the expected location of tex2html_wrap_inline2157 , and tex2html_wrap_inline2173 is an expression of how precise that estimate is. In general, tex2html_wrap_inline2173 will be small.

Thus, having collected data, it is possible to derive the posterior for tex2html_wrap_inline2157 according to Bayes theorem for random variables;

eqnarray243

where tex2html_wrap_inline2185 is independent of tex2html_wrap_inline2157 . Defining

displaymath2147

and multiplying by tex2html_wrap_inline2189 which is independent of tex2html_wrap_inline2157 , the above is

displaymath2148

which reduces to

displaymath2149

which is the form of the normal density with mean tex2html_wrap_inline2193 and precision tex2html_wrap_inline2195 . Thus, in the case of inference for the unknown mean, with normal prior, the posterior is normal. This simple form of the posterior depends on the choice of the prior, given the likelihood. The choice of prior that leads to the simple posterior, is called a conjugate prior; more formally, given a likelihood, tex2html_wrap_inline2073 , then a prior chosen from a family of densities, such that the posterior is also from that family, is said to be conjugate.

As can be seen from the above, in the case of conjugate densities, the problem of obtaining a posterior is simplified [5]. However, this is only appropriate where the chosen prior distribution, with suitable parameters can accurately represent the prior knowledge. The alternative is to use numerical techniques to obtain the properties of interest from the posterior distribution.

The question of prior elicitation is one that needs mentioning also. Apart from the philosophical difficulties that many have with prior probabilities, there are practical problems which need addressing.


next up previous contents
Next: Prior Elicitation and Non-informative Up: Bayesian Approach Previous: Kernel Density Estimation

Cathal Walsh
Sat Jan 22 17:09:53 GMT 2000