Describe the Weak Sufficiency Principle (WSP) and the Strong Sufficiency Principle (SSP). Show that Bayesian inference with a fixed prior distribution respects WSP.
A parameter Φ has a prior distribution which is normal with mean 0 and precision (inverse variance) hΦ. Given Φ=ϕ, further parameters Θ:=(Θi:i=1,…,I) have independent normal distributions with mean ϕ and precision hΘ. Finally, given both Φ=ϕ and Θ=θ:=(θ1,…,θI), observables X:=(Xij:i=1,…,I;j=1,…,J) are independent, Xij being normal with mean θi, and precision hX. The precision parameters (hΦ,hΘ,hX) are all fixed and known. Let X:=(Xˉ1,…,XˉI), where Xˉi:=∑j=1JXij/J. Show, directly from the definition of sufficiency, that X is sufficient for (Φ,Θ). [You may assume without proof that, if Y1,…,Yn have independent normal distributions with the same variance, and Yˉ:=n−1∑i=1nYi, then the vector (Y1−Yˉ,…,Yn−Yˉ) is independent of Yˉ.]
For data-values x:=(xij:i=1,…,I;j=1,…,J), determine the joint distribution, Πϕ say, of Θ, given X=x and Φ=ϕ. What is the distribution of Φ, given Θ=θ and X=x?
Using these results, describe clearly how Gibbs sampling combined with RaoBlackwellisation could be applied to estimate the posterior joint distribution of Θ, given X=x.