Section 4.2

First Examples

[16/0%]
LaTeX
Problem 4.2.1

Referring to Example 1.5, suppose that XX has the binomial distribution b(p,n)b(p, n) and the family of prior distributions for pp is the family of beta distributions B(a,b)B(a, b).

?
(a)

Show that the marginal distribution of XX is the beta-binomial distribution with mass function

(nx)Γ(a+b)Γ(a)Γ(b)Γ(x+a)Γ(n−x+b)Γ(n+a+b). \binom {n}{x} \frac{\Gamma (a+b)}{\Gamma (a) \Gamma (b)} \frac{\Gamma (x+a) \Gamma (n-x+b)}{\Gamma (n+a+b)}.
(b)

Show that the mean and variance of the beta-binomial is given by

E[X]=naa+b and  var X=n(aa+b)(ba+b)(a+b+na+b+1). \mathbb {E}\left[X\right]=\frac{n a}{a+b} \quad \text{ and } \quad \text{ var } X=n\left(\frac{a}{a+b}\right)\left(\frac{b}{a+b}\right)\left(\frac{a+b+n}{a+b+1}\right).

[Hint: For part (b), the identities E[X]=E[E[X∣p]]\mathbb {E}\left[X\right]=E[\mathbb {E}\left[X \mid p\right]] and var⁡X=var⁡[E[X∣p]]+E[var⁡(X∣p)]\operatorname {var} X=\operatorname {var}[\mathbb {E}\left[X \mid p\right]]+ E[\operatorname {var}(X \mid p)] are helpful.]

Problem 4.2.2

For the situation of Example 2.1, Lindley and Phillips (1976) give a detailed account of the effect of stopping rules, which we can illustrate as follows. Let XX be the number of successes in nn Bernoulli trials with success probability pp.

?
(a)

Suppose that the number of Bernoulli trials performed is a prespecifed number nn, so that we have the binomial sampling model, P(X=x)=(nx)px(1−p)n−x\mathbb {P}\left(X=x\right)=\binom {n}{x} p^{x}(1-p)^{n-x}, x=0,1,…,nx=0,1, \ldots , n. Calculate the Bayes risk of the Bayes estimator (1.1.12) and the UMVU estimator of pp.

(b)

Suppose that the number of Bernoulli trials performed is a random variable NN. The value N=nN=n was obtained when a prespecified number, xx, of successes was observed so that we have the negative binomial sample model, P(N=n)=(n−1x−1)px(1−p)n−x\mathbb {P}\left(N=n\right)=\binom {n-1}{x-1} p^{x}(1- p)^{n-x}, n=xn=x. Calculate the Bayes risk of the Bayes estimator and the UMVU estimator of pp.

(c)

Calculate the mean squared errors of all three estimators under each model. If it is unknown which sampling mechanism generated the data, which estimator do you prefer overall?

Problem 4.2.3

Show that the estimator (2.2.4) tends in probability

?
(a)

to θ\theta as n→∞n \rightarrow \infty,

(b)

to μ\mu as b→0b \rightarrow 0, and

(c)

to θ\theta as b→∞b \rightarrow \infty.

Problem 4.2.4

Bickel and Mallows (1988) further investigate the relationship between unbiasedness and Bayes, specifying conditions under which these properties cannot hold simultaneously. In addition, they show that if a prior distribution is improper, then a posterior mean can be unbiased. Let X∼1θf(x/θ),x>0X \sim \frac{1}{\theta } f(x / \theta ), x>0, where ∫0∞tf(t)dt=1\int_{0}^{\infty } t f(t) d t=1, and let π(θ)=1θ2dθ\pi (\theta )=\frac{1}{\theta^{2}} d \theta, θ>0\theta >0.

?
(a)

Show that E[X∣θ]=θ\mathbb {E}\left[X \mid \theta \right]=\theta, so XX is unbiased.

(b)

Show that π(θ∣x)=x2θ3f(x/θ)\pi (\theta \mid x)=\frac{x^{2}}{\theta^{3}} f(x / \theta ) is a proper density.

(c)

Show that E[θ∣x]=x\mathbb {E}\left[\theta \mid x\right]=x, and hence the posterior mean, is unbiased.

Problem 4.2.5

DasGupta (1994) presents an identity relating the Bayes risk to bias, which illustrates that a small bias can help achieve a small Bayes risk. Let X∼f(x∣θ)X \sim f(x \mid \theta ) and θ∼π(θ)\theta \sim \pi (\theta ). The Bayes estimator under squared error loss is δπ=E[θ∣x]\delta^{\pi }=\mathbb {E}\left[\theta \mid x\right]. Show that the Bayes risk of δπ\delta^{\pi } can be written

r(π,δπ)=∫Θ∫X[θ−δπ(x)]2f(x∣θ)π(θ)dxdθ=∫Θθb(θ)π(θ)dθ r\left(\pi , \delta ^{\pi }\right)=\int _{\Theta } \int _{\mathcal{X}}\left[\theta -\delta ^{\pi }(x)\right]^{2} f(x \mid \theta ) \pi (\theta ) d x d \theta =\int _{\Theta } \theta b(\theta ) \pi (\theta ) d \theta

where b(θ)=E[δπ(X)∣θ]−θb(\theta )=\mathbb {E}\left[\delta^{\pi }(X) \mid \theta \right]-\theta is the bias of δπ\delta^{\pi }.

?
Problem 4.2.6

Verify the estimator (2.2.10).

?
Problem 4.2.7

In Example 2.6, verify that the posterior distribution of τ\tau is Γ(r+g−1/2,1/(α+z))\Gamma (r+g-1 / 2,1 /(\alpha +z)).

?
Problem 4.2.8

In Example 2.6 with α=g=0\alpha =g=0, show that the posterior distribution given the XX 's of n(θ−Xˉ)/Z/(n−1)\sqrt{n}(\theta -\bar{X}) / \sqrt{Z /(n-1)} is Student's tt-distribution with n−1n-1 degrees of freedom.

?
Problem 4.2.9

In Example 2.6, show that the posterior distribution of θ\theta is symmetric about xˉ\bar{x} when the joint prior of θ\theta and σ\sigma is of the form h(σ)dσdθh(\sigma ) d \sigma d \theta, where hh is an arbitrary probability density on (0,∞)(0, \infty ).

?
Problem 4.2.10

Rukhin (1978) investigates the situation when the Bayes estimator is the same for every loss function in a certain set of loss functions, calling such estimators universal Bayes estimators. For the case of Example 2.6, using the prior of the form of Problem 2.9, show that Xˉ\bar{X} is the Bayes estimator under every even loss function.

?
Problem 4.2.11

Let XX and YY be independently distributed according to distributions PξP_{\xi } and QηQ_{\eta }, respectively. Suppose that ξ\xi and η\eta are real-valued and independent according to some prior distributions Λ\Lambda and Λ′\Lambda^{\prime }. If, with squared error loss, δΛ\delta_{\Lambda } is the Bayes estimator of ξ\xi on the basis of XX, and δΛ′′\delta_{\Lambda^{\prime }}^{\prime } is that of η\eta on the basis of YY,

?
(a)

show that δΛ′′−δΛ\delta_{\Lambda^{\prime }}^{\prime }-\delta_{\Lambda } is the Bayes estimator of η−ξ\eta -\xi on the basis of (X,Y)(X, Y);

(b)

if η>0\eta >0 and δΛ′∗\delta_{\Lambda^{\prime }}^{*} is the Bayes estimator of 1/η1 / \eta on the basis of YY, show that δΛ⋅δΛ′∗\delta_{\Lambda } \cdot \delta_{\Lambda^{\prime }}^{*} is the Bayes estimator of ξ/η\xi / \eta on the basis of (X,Y)(X, Y).

Problem 4.2.12

For the density (2.2.13) and improper prior (dσ/σ)⋅(dσA/σA)(d \sigma / \sigma ) \cdot \left(d \sigma_{A} / \sigma_{A}\right), show that the posterior distribution of (σ,σA)\left(\sigma , \sigma_{A}\right) continues to be improper.

?
Problem 4.2.13
?
(a)

In Example 2.7, obtain the Jeffreys prior distribution of (σ,τ)(\sigma , \tau ).

(b)

Show that for the prior of part (a), the posterior distribution of (σ,τ)(\sigma , \tau ) is proper.

Problem 4.2.14

Verify the Bayes estimator (2.2.14).

?
Problem 4.2.15

Let X∼N(θ,1)X \sim N(\theta , 1) and L(θ,δ)=(θ−δ)2L(\theta , \delta )=(\theta -\delta )^{2}.

?
(a)

Show that XX is the limit of the Bayes estimators δπn\delta^{\pi_{n}}, where πn\pi_{n} is N(0,1)N(0,1). Hence, XX is both generalized Bayes and a limit of Bayes estimators.

(b)

For the prior measure π(θ)=eaθ,a>0\pi (\theta )=e^{a \theta }, a>0, show that the generalized Bayes estimator is X+aX+a.

(c)

For a>0a>0, show that there is no sequence of proper priors for which δπn→X+a\delta^{\pi_{n}} \rightarrow X+a.

This example is due to Farrell; see Kiefer 1966. Heath and Sudderth (1989), building on the work of Stone (1976), showed that inferences from this model are incoherent, and established when generalized Bayes estimators will lead to coherent (that is, noncontradictory) inferences. Their work is connected to the theory of "approximable by proper priors," developed by Stein (1965) and Stone (1965, 1970, 1976), which shows when generalized Bayes estimators can be looked upon as Bayes estimators.

Problem 4.2.16
?
(a)

For the situation of Example 2.8, verify that δ(x)=x/n\delta (x)=x / n is a generalized Bayes estimator.

(b)

If X∼N(0,1)X \sim N(0,1) and L(θ,δ)=(θ−δ)2L(\theta , \delta )=(\theta -\delta )^{2}, show that XX is generalized Bayes under the improper prior π(θ)=1\pi (\theta )=1.