4.2

Первые примеры

[16/0%]
LaTeX
Задача 4.2.1

Обращаясь к Примеру 1.5, предположим, что XX имеет биномиальное распределение b(p,n)b(p, n), а семейство априорных распределений для pp является семейством бета-распределений B(a,b)B(a, b).

?
(a)

Покажите, что маргинальное распределение XX является бета-биномиальным распределением с функцией вероятности

(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)

Покажите, что среднее и дисперсия бета-биномиального распределения задаются формулами

E[X]=naa+b и  var X=n(aa+b)(ba+b)(a+b+na+b+1). \mathbb {E}\left[X\right]=\frac{n a}{a+b} \quad \text{ и } \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).

[Указание: Для пункта (b) полезны тождества E[X]=E[E[X∣p]]\mathbb {E}\left[X\right]=E[\mathbb {E}\left[X \mid p\right]] и 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)].]

Задача 4.2.2

Для ситуации Примера 2.1, Lindley и Phillips (1976) приводят подробное описание влияния правил остановки, которое мы можем проиллюстрировать следующим образом. Пусть XX — число успехов в nn испытаниях Бернулли с вероятностью успеха pp.

?
(a)

Предположим, что число проведённых испытаний Бернулли является заранее заданным числом nn, так что мы имеем биномиальную модель выборки, 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. Вычислите байесовский риск байесовской оценки (1.1.12) и UMVU-оценки pp.

(b)

Предположим, что число проведённых испытаний Бернулли является случайной величиной NN. Значение N=nN=n было получено, когда было замечено заранее заданное число xx успехов, так что мы имеем модель отрицательной биномиальной выборки, 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. Вычислите байесовский риск байесовской оценки и UMVU-оценки pp.

(c)

Вычислите среднеквадратичные ошибки всех трёх оценок для каждой модели. Если неизвестно, какой механизм выборки породил данные, какую оценку вы предпочтёте в целом?

Задача 4.2.3

Покажите, что оценка (2.2.4) сходится по вероятности

?
(a)

к θ\theta при n→∞n \rightarrow \infty,

(b)

к μ\mu при b→0b \rightarrow 0, и

(c)

к θ\theta при b→∞b \rightarrow \infty.

Задача 4.2.4

Bickel и Mallows (1988) более подробно исследуют связь между несмещённостью и байесовским подходом, указывая условия, при которых эти свойства не могут выполняться одновременно. Кроме того, они показывают, что если априорное распределение несобственно, то апостериорное среднее может быть несмещённым. Пусть X∼1θf(x/θ),x>0X \sim \frac{1}{\theta } f(x / \theta ), x>0, где ∫0∞tf(t)dt=1\int_{0}^{\infty } t f(t) d t=1, и пусть π(θ)=1θ2dθ\pi (\theta )=\frac{1}{\theta^{2}} d \theta, θ>0\theta >0.

?
(a)

Покажите, что E[X∣θ]=θ\mathbb {E}\left[X \mid \theta \right]=\theta, так что XX несмещена.

(b)

Покажите, что π(θ∣x)=x2θ3f(x/θ)\pi (\theta \mid x)=\frac{x^{2}}{\theta^{3}} f(x / \theta ) является собственной плотностью.

(c)

Покажите, что E[θ∣x]=x\mathbb {E}\left[\theta \mid x\right]=x, и, следовательно, апостериорное среднее несмещено.

Задача 4.2.5

DasGupta (1994) приводит тождество, связывающее байесовский риск со смещением, которое иллюстрирует, что небольшое смещение может помочь достичь малого байесовского риска. Пусть X∼f(x∣θ)X \sim f(x \mid \theta ) и θ∼π(θ)\theta \sim \pi (\theta ). Байесовская оценка при квадратичной функции потерь равна δπ=E[θ∣x]\delta^{\pi }=\mathbb {E}\left[\theta \mid x\right]. Покажите, что байесовский риск δπ\delta^{\pi } может быть записан как

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

где b(θ)=E[δπ(X)∣θ]−θb(\theta )=\mathbb {E}\left[\delta^{\pi }(X) \mid \theta \right]-\theta — смещение δπ\delta^{\pi }.

?
Задача 4.2.6

Проверьте оценку (2.2.10).

?
Задача 4.2.7

В Примере 2.6 проверьте, что апостериорное распределение τ\tau есть Γ(r+g−1/2,1/(α+z))\Gamma (r+g-1 / 2,1 /(\alpha +z)).

?
Задача 4.2.8

В Примере 2.6 при α=g=0\alpha =g=0 покажите, что апостериорное распределение n(θ−Xˉ)/Z/(n−1)\sqrt{n}(\theta -\bar{X}) / \sqrt{Z /(n-1)} при заданных XX является tt-распределением Стьюдента с n−1n-1 степенями свободы.

?
Задача 4.2.9

В Примере 2.6 покажите, что апостериорное распределение θ\theta симметрично относительно xˉ\bar{x}, когда совместное априорное распределение θ\theta и σ\sigma имеет вид h(σ)dσdθh(\sigma ) d \sigma d \theta, где hh — произвольная плотность вероятности на (0,∞)(0, \infty ).

?
Задача 4.2.10

Rukhin (1978) исследует ситуацию, когда байесовская оценка одна и та же для каждой функции потерь из некоторого множества функций потерь, называя такие оценки универсальными байесовскими оценками. Для случая Примера 2.6, используя априорное распределение вида, указанного в Задаче 2.9, покажите, что Xˉ\bar{X} является байесовской оценкой при любой чётной функции потерь.

?
Задача 4.2.11

Пусть XX и YY независимо распределены согласно распределениям PξP_{\xi } и QηQ_{\eta } соответственно. Предположим, что ξ\xi и η\eta являются вещественнозначными и независимы согласно некоторым априорным распределениям Λ\Lambda и Λ′\Lambda^{\prime }. Если при квадратичной функции потерь δΛ\delta_{\Lambda } — байесовская оценка ξ\xi на основе XX, а δΛ′′\delta_{\Lambda^{\prime }}^{\prime } — байесовская оценка η\eta на основе YY,

?
(a)

покажите, что δΛ′′−δΛ\delta_{\Lambda^{\prime }}^{\prime }-\delta_{\Lambda } является байесовской оценкой η−ξ\eta -\xi на основе (X,Y)(X, Y);

(b)

если η>0\eta >0 и δΛ′∗\delta_{\Lambda^{\prime }}^{*} — байесовская оценка 1/η1 / \eta на основе YY, покажите, что δΛ⋅δΛ′∗\delta_{\Lambda } \cdot \delta_{\Lambda^{\prime }}^{*} является байесовской оценкой ξ/η\xi / \eta на основе (X,Y)(X, Y).

Задача 4.2.12

Для плотности (2.2.13) и несобственного априорного распределения (dσ/σ)⋅(dσA/σA)(d \sigma / \sigma ) \cdot \left(d \sigma_{A} / \sigma_{A}\right) покажите, что апостериорное распределение (σ,σA)\left(\sigma , \sigma_{A}\right) остаётся несобственным.

?
Задача 4.2.13
?
(a)

В Примере 2.7 получите априорное распределение Джеффриса для (σ,τ)(\sigma , \tau ).

(b)

Покажите, что для априорного распределения из пункта (a) апостериорное распределение (σ,τ)(\sigma , \tau ) является собственным.

Задача 4.2.14

Проверьте байесовскую оценку (2.2.14).

?
Задача 4.2.15

Пусть X∼N(θ,1)X \sim N(\theta , 1) и L(θ,δ)=(θ−δ)2L(\theta , \delta )=(\theta -\delta )^{2}.

?
(a)

Покажите, что XX является пределом байесовских оценок δπn\delta^{\pi_{n}}, где πn\pi_{n} есть N(0,1)N(0,1). Следовательно, XX является одновременно обобщённой байесовской оценкой и пределом байесовских оценок.

(b)

Для априорной меры π(θ)=eaθ,a>0\pi (\theta )=e^{a \theta }, a>0, покажите, что обобщённая байесовская оценка равна X+aX+a.

(c)

Для a>0a>0 покажите, что не существует последовательности собственных априорных распределений, для которой δπn→X+a\delta^{\pi_{n}} \rightarrow X+a.

Этот пример принадлежит Farrell; см. Kiefer 1966. Heath и Sudderth (1989), опираясь на работу Stone (1976), показали, что выводы из этой модели некогерентны, и установили, когда обобщённые байесовские оценки приводят к когерентным (то есть непротиворечивым) выводам. Их работа связана с теорией «аппроксимируемости собственными априорными распределениями», развитой Stein (1965) и Stone (1965, 1970, 1976), которая показывает, когда обобщённые байесовские оценки можно рассматривать как байесовские оценки.

Задача 4.2.16
?
(a)

Для ситуации Примера 2.8 проверьте, что δ(x)=x/n\delta (x)=x / n является обобщённой байесовской оценкой.

(b)

Если X∼N(0,1)X \sim N(0,1) и L(θ,δ)=(θ−δ)2L(\theta , \delta )=(\theta -\delta )^{2}, покажите, что XX является обобщённой байесовской оценкой при несобственном априорном распределении π(θ)=1\pi (\theta )=1.