Глава 1

Предварительные сведения

[173/0%]
Показать
LaTeX
§
Задача 1.1.1

Если (x1,y1),…,(xn,yn)\left(x_{1}, y_{1}\right), \ldots ,\left(x_{n}, y_{n}\right) — это nn точек на плоскости, найдите наилучшую в смысле метода наименьших квадратов прямую y=α+βxy= \alpha +\beta x, то есть найдите значения α\alpha и β\beta, минимизирующие Σ[yi−(α+βxi)]2\Sigma \left[y_{i}-\left(\alpha +\beta x_{i}\right)\right]^{2}.

?
Задача 1.1.2

Пусть X1,…,XnX_{1}, \ldots , X_{n} — некоррелированные случайные величины с общим математическим ожиданием θ\theta и дисперсией σ2\sigma^{2}. Тогда среди всех линейных оценок ΣαiXi\Sigma \alpha_{i} X_{i} параметра θ\theta, удовлетворяющих условию Σαi=1\Sigma \alpha_{i}=1, среднее Xˉ\bar{X} имеет наименьшую дисперсию.

?
Задача 1.1.3

В предыдущей задаче минимизируйте дисперсию ΣαiXi(Σαi=1)\Sigma \alpha_{i} X_{i}\left(\Sigma \alpha_{i}=1\right)

?
(a)

Когда дисперсия XiX_{i} равна σ2/αi\sigma^{2} / \alpha_{i} ( αi\alpha_{i} известно).

(b)

Когда XiX_{i} имеют общую дисперсию σ2\sigma^{2}, но коррелированы с общим коэффициентом корреляции ρ\rho.

(Обобщения этих результатов см., например, Watson 1967 и Kruskal 1968.)

Задача 1.1.4

Пусть XX и YY имеют общее математическое ожидание θ\theta, дисперсии σ2\sigma^{2} и τ2\tau^{2} и коэффициент корреляции ρ\rho. Определите условия на σ,τ\sigma , \tau и ρ\rho, при которых

?
(a)

var⁡(X)<var⁡[(X+Y)/2]\operatorname {var}(X)<\operatorname {var}[(X+Y) / 2].

(b)

Значение α\alpha, минимизирующее var⁡[αX+(1−α)Y]\operatorname {var}[\alpha X+(1-\alpha ) Y], отрицательно.

Дайте интуитивное объяснение полученным результатам.

Задача 1.1.5

Пусть Xi(i=1,2)X_{i}(i=1,2) независимо распределены согласно плотностям Коши C(ai,bi)C\left(a_{i}, b_{i}\right). Тогда X1+X2X_{1}+X_{2} распределена как C(a1+a2,b1+b2)C\left(a_{1}+a_{2}, b_{1}+b_{2}\right). [Указание: Перейдите к новым переменным Y1=X1+X2,Y2=X2Y_{1}=X_{1}+X_{2}, Y_{2}=X_{2}.]

?
Задача 1.1.6

Если X1,…,XnX_{1}, \ldots , X_{n} независимые одинаково распределённые как C(a,b)C(a, b), то распределение Xˉ\bar{X} снова является C(a,b)C(a, b). [Указание: Докажите по индукции, используя Задачу 5.]

?
Задача 1.1.7

Медианой XX называется любое значение mm, такое что P(X≤m)≥1/2P(X \leq m) \geq 1 / 2 и P(X≥m)≥1/2P(X \geq m) \geq 1 / 2.

?
(a)

Покажите, что это эквивалентно условиям P(X<m)≤1/2P(X<m) \leq 1 / 2 и P(X>m)≤1/2P(X>m) \leq 1 / 2.

(b)

Покажите, что множество медиан всегда является замкнутым интервалом m0≤m≤m1m_{0} \leq m \leq m_{1}.

Задача 1.1.8

Если ϕ(a)=E∣X−a∣<∞\phi (a)=E\left|X-a\right|<\infty для некоторого aa, покажите, что ϕ(a)\phi (a) минимизируется любой медианой XX. [Указание: Если m0≤m≤m1m_{0} \leq m \leq m_{1} (в обозначениях Задачи 1.7) и m1<cm_{1}<c, то

E∣X−c∣−E∣X−m∣=(c−m)[P(X≤m)−P(X>m)]+2∫m<x<c(c−x)dP(x) E\left|X-c\right|-E\left|X-m\right|=(c-m)[P(X \leq m)-P(X>m)]+2 \int _{m<x<c}(c-x) d P(x)
?
Задача 1.1.9
?
(a)

Медианой любого набора различных действительных чисел x1,…,xnx_{1}, \ldots , x_{n} называется средний элемент упорядоченных xx, если nn нечётно, и любое значение между двумя средними упорядоченными xx, если nn чётно. Покажите, что это также медиана случайной величины XX, принимающей каждое из значений x1,…,xnx_{1}, \ldots , x_{n} с вероятностью 1/n1 / n.

(b)

Для любого набора различных действительных чисел x1,…,xnx_{1}, \ldots , x_{n} сумма абсолютных отклонений Σ∣xi−a∣\Sigma \left|x_{i}-a\right| минимизируется любой медианой чисел xx.

(c)

Для nn заданных точек (xi,yi),i=1,…,n\left(x_{i}, y_{i}\right), i=1, \ldots , n, найдите значение bb, минимизирующее Σ∣yi−bxi∣\Sigma \left|y_{i}-b x_{i}\right|. [Указание: Сведите задачу к частному случаю Задачи 8.]

Задача 1.1.10

Для любого набора чисел x1,⋯ ,xnx_{1}, \cdots , x_{n} и монотонной функции h(⋅)h(\cdot ) покажите, что значение aa, минимизирующее ∑i=1n[h(xi)−h(a)]2\sum_{i=1}^{n}\left[h\left(x_{i}\right)-h(a)\right]^{2}, задаётся формулой a=h−1(∑i=1nh(xi)/n)a=h^{-1}\left(\sum_{i=1}^{n} h\left(x_{i}\right) / n\right). Найдите функции hh, при которых минимизирующее значение даёт среднее арифметическое, среднее геометрическое и среднее гармоническое.

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

Если две оценки δ1,δ2\delta_{1}, \delta_{2} имеют непрерывные симметричные плотности fi(x−θ),i=1,2f_{i}(x-\theta ), i=1,2, и f1(0)>f2(0)f_{1}(0)>f_{2}(0), то

P[∣δ1−θ∣<c]>P[∣δ2−θ∣<c] для некоторого c>0 P\left[\left|\delta _{1}-\theta \right|<c\right]>P\left[\left|\delta _{2}-\theta \right|<c\right] \text{ для некоторого } c>0

и, следовательно, δ1\delta_{1} будет ближе к θ\theta, чем δ2\delta_{2}, по отношению к мере (1.5).

(b)

Пусть X,YX, Y независимо распределены с общей непрерывной симметричной плотностью ff, и пусть δ1=X,δ2=(X+Y)/2\delta_{1}=X, \delta_{2}=(X+Y) / 2. Неравенство из пункта (a) будет выполняться, если 2∫f2(x)dx<f(0)2 \int f^{2}(x) d x<f(0) (Edgeworth 1883, Stigler 1980).

Задача 1.1.12
?
(a)

Пусть f(x)=(1/2)(k−1)/(1+∣x∣)k,k≥2f(x)=(1 / 2)(k-1) /(1+\left|x\right|)^{k}, k \geq 2. Покажите, что ff является плотностью вероятности и что все её моменты порядка <k−1<k-1 конечны.

(b)

Плотность из пункта (a) удовлетворяет неравенству Задачи 1.11(b).

Задача 1.1.13
?
(a)

Если XX биномиальна b(p,n)b(p, n), покажите, что

E∣xn−p∣=2(n−1k−1)pk(1−p)n−k+1 для k−1n≤p≤kn. E\left|\frac{x}{n}-p\right|=2\binom {n-1}{k-1} p^{k}(1-p)^{n-k+1} \quad \text{ для } \frac{k-1}{n} \leq p \leq \frac{k}{n}.
(b)

Постройте график функции риска из пункта (i) для n=4n=4 и n=5n=5.

§
Задача 1.2.1

Если A1,A2,…A_{1}, A_{2}, \ldots — элементы σ\sigma-алгебры A\mathcal{A} (множества AA не обязательно непересекающиеся), то их объединение и пересечение также принадлежат A\mathcal{A}.

?
Задача 1.2.2

Для любых a<ba<b следующие множества являются борелевскими

?
(a)

{x:a<x}\left\{ x: a<x\right\} и

(b)

{x:a≤x≤b}\left\{ x: a \leq x \leq b\right\}.

Задача 1.2.3

В условиях Задачи 2.1 пусть

A‾=lim inf⁡An={x:x∈An для всех, кроме конечного числа n },Aˉ=lim sup⁡An={x:x∈An для бесконечного числа n}. \begin{align} \underline{A} & =\liminf A_{n}=\left\{ x: x \in A_{n} \text{ для всех, кроме конечного числа } n\text{ }\right\} , \\ \bar{A} & =\limsup A_{n}=\left\{ x: x \in A_{n} \text{ для бесконечного числа } n\right\} . \end{align}

Тогда A‾\underline{A} и Aˉ\bar{A} принадлежат A\mathcal{A}.

?
Задача 1.2.4

Покажите, что

?
(a)

Если A1⊂A2⊂⋯A_{1} \subset A_{2} \subset \cdots, то A‾=Aˉ=∪An\underline{A}=\bar{A}=\cup A_{n}.

(b)

Если A1⊃A2⊃⋯A_{1} \supset A_{2} \supset \cdots, то A‾=Aˉ=∩An\underline{A}=\bar{A}=\cap A_{n}.

Задача 1.2.5

Для любой последовательности действительных чисел a1,a2,…a_{1}, a_{2}, \ldots покажите, что множество всех предельных точек подпоследовательностей замкнуто. Наименьшая и наибольшая из таких предельных точек (которые могут быть бесконечными) обозначаются через lim inf⁡ak\liminf a_{k} и lim sup⁡ak\limsup a_{k} соответственно.

?
Задача 1.2.6

В условиях Задач 2.1 и 2.3 покажите, что

IA‾(x)=lim inf⁡IAk(x) и IAˉ(x)=lim⁡sup⁡IAk(x) I_{\underline{A}}(x)=\liminf I_{A_{k}}(x) \quad \text{ и } \quad I_{\bar{A}}(x)=\lim \sup I_{A_{k}}(x)

где IA(x)I_{A}(x) обозначает индикатор множества AA.

?
Задача 1.2.7

Пусть (X,A,μ)(\mathcal{X}, \mathcal{A}, \mu ) — пространство с мерой, и пусть B\mathcal{B} — класс всех множеств A∪CA \cup C, где A∈AA \in \mathcal{A}, а CC — подмножество некоторого множества A′∈AA^{\prime } \in \mathcal{A} с μ(A′)=0\mu \left(A^{\prime }\right)=0. Покажите, что B\mathcal{B} является σ\sigma-алгеброй.

?
Задача 1.2.8

Если ff и gg — измеримые функции, то измеримыми являются и

  1. f+gf+g, и

  2. max⁡(f,g)\max (f, g).

?
Задача 1.2.9

Если ff интегрируема относительно μ\mu, то и ∣f∣\left|f\right| интегрируема, причём ∣∫fdμ∣≤∫∣f∣dμ\left|\int f d \mu \right| \leq \int \left|f\right| d \mu. [Указание: Выразите ∣f∣\left|f\right| через f+f^{+} и f−f^{-}.]

?
Задача 1.2.10

Пусть X={x1,x2,…},μ=\mathcal{X}=\left\{ x_{1}, x_{2}, \ldots \right\} , \mu = считающая мера на X\mathcal{X}, и ff интегрируема. Тогда ∫fdμ=Σf(xi)\int f d \mu = \Sigma f\left(x_{i}\right). [Указание: Предположим сначала, что f≥0f \geq 0, и пусть sn(x)s_{n}(x) — простая функция, равная f(x)f(x) при x=x1,…,xnx=x_{1}, \ldots , x_{n}, и 0 в остальных случаях.]

?
Задача 1.2.11

Пусть f(x)=1f(x)=1 или 0 в зависимости от того, рационально xx или иррационально. Покажите, что интеграл Римана от ff не существует.

?
§
Задача 1.3.1

Пусть XX имеет стандартное нормальное распределение, и пусть Y=2XY=2 X. Определите, является ли

?
(a)

функция распределения F(x,y)F(x, y) пары (X,Y)(X, Y) непрерывной.

(b)

распределение (X,Y)(X, Y) абсолютно непрерывным относительно меры Лебега на плоскости (x,y)(x, y).

Задача 1.3.2

Покажите, что любая функция ff, удовлетворяющая (3.7), непрерывна.

?
Задача 1.3.3

Пусть XX — измеримое преобразование из (E,B)(\mathcal{E}, \mathcal{B}) в (X,A)(\mathcal{X}, \mathcal{A}) (то есть такое, что для любого A∈AA \in \mathcal{A} множество {e:X(e)∈A}\left\{ e: X(e) \in A\right\} принадлежит B\mathcal{B}), и пусть YY — измеримое преобразование из (X,A)(\mathcal{X}, \mathcal{A}) в (Y,C)(\mathcal{Y}, \mathcal{C}). Тогда Y[X(e)]Y[X(e)] является измеримым преобразованием из (E,B)(\mathcal{E}, \mathcal{B}) в (Y,C)(\mathcal{Y}, \mathcal{C}).

?
Задача 1.3.4

В Примере 3.1 покажите, что носитель PP равен [a,b][a, b] тогда и только тогда, когда FF строго возрастает на [a,b][a, b].

?
Задача 1.3.5

Пусть SS — носитель распределения на евклидовом пространстве (X,A)(\mathcal{X}, \mathcal{A}). Тогда,

  1. SS замкнуто;

  2. P(S)=1P(S)=1;

  3. SS является пересечением всех замкнутых множеств CC с P(C)=1P(C)=1.

?
Задача 1.3.6

Если PP и QQ — две эквивалентные вероятностные меры на одном и том же евклидовом пространстве, то они имеют один и тот же носитель.

?
Задача 1.3.7

Пусть PP и QQ задают вероятности

P:P(X=1n)=pn>0,n=1,2,…(Σpn=1),Q:P(X=0)=12;P(X=1n)=qn>0;n=1,2,…(Σqn=12). \begin{align} & P: P\left(X=\frac{1}{n}\right)=p_{n}>0, \quad n=1,2, \ldots \quad \left(\Sigma p_{n}=1\right), \\ & Q: P(X=0)=\frac{1}{2} ; \quad P\left(X=\frac{1}{n}\right)=q_{n}>0 ; \quad n=1,2, \ldots \quad \left(\Sigma q_{n}=\frac{1}{2}\right). \end{align}

Тогда покажите, что PP и QQ имеют один и тот же носитель, но не являются эквивалентными.

?
Задача 1.3.8

Пусть XX и YY — независимые случайные величины с X∼E(λ,1)X \sim E(\lambda , 1) и Y∼E(μ,1)Y \sim E(\mu , 1). Непосредственные наблюдения XX и YY получить невозможно. Вместо этого мы наблюдаем случайные величины ZZ и WW, где

Z=min⁡{X,Y} и W={1 если Z=X0 если Z=Y. Z=\min \left\{ X, Y\right\} \quad \text{ и } \quad W= \begin{cases} 1 & \text{ если } Z=X \\ 0 & \text{ если } Z=Y.\end{cases}

Найдите совместное распределение ZZ и WW и покажите, что они независимы. (Переменные XX и YY цензурированы — ситуация, часто встречающаяся в медицинских экспериментах. Предположим, что XX измеряет время выживания при некотором лечении, а пациент выбывает из исследования по причине, не связанной с этим. Мы не получаем измерения XX, а только нижнюю границу.)

?
§
Задача 1.4.1

Если распределения положительной случайной величины XX образуют семейство масштаба, покажите, что распределения log⁡X\log X образуют семейство сдвига.

?
Задача 1.4.2

Если XX распределена согласно равномерному распределению U(0,θ)U(0, \theta ), покажите, что распределение −log⁡X-\log X является показательным.

?
Задача 1.4.3

Пусть UU равномерно распределена на (0,1)(0,1), и рассмотрим переменные X=Uα,0<αX=U^{\alpha }, 0<\alpha. Покажите, что это определяет групповое семейство, и найдите плотность XX.

?
Задача 1.4.4

Покажите, что группа преобразований является группой.

?
Задача 1.4.5

Если g0g_{0} — произвольный элемент группы GG, покажите, что когда gg пробегает GG, то и gg0g g_{0} пробегает GG.

?
Задача 1.4.6

Покажите, что при p=2p=2 плотность (4.15) переходит в (4.16).

?
Задача 1.4.7

Покажите, что семейство преобразований (4.12) с невырожденной нижнетреугольной матрицей BB образует группу GG.

?
Задача 1.4.8

Покажите, что совокупность всех невырожденных многомерных нормальных распределений может быть получена с помощью подгруппы GG из (4.12), описанной в Задаче 4.7.

?
Задача 1.4.9

В предыдущей задаче покажите, что GG можно заменить подгруппой G0G_{0} нижнетреугольных матриц B=(bij)B=\left(b_{i j}\right), в которых диагональные элементы b11,…,bppb_{11}, \ldots , b_{p p} все положительны, но никакая собственная подгруппа G0G_{0} для этого не подойдёт.

?
Задача 1.4.10

Покажите, что семейство всех непрерывных распределений, носитель которых представляет собой интервал с положительной нижней границей, является групповым семейством. [Указание: Пусть UU равномерно распределена на интервале (2,3)(2,3), и пусть X=b[g(U)]aX=b[g(U)]^{a}, где α,b>0\alpha , b>0, а gg непрерывна и взаимно однозначна из (2,3)(2,3) в (2,3)(2,3).]

?
Задача 1.4.11

Найдите модификацию группы преобразований (4.22), порождающую случайную выборку из генеральной совокупности {y1,…,yN}\left\{ y_{1}, \ldots , y_{N}\right\}, в которой значения yy, вместо того чтобы быть произвольными, ограничены условием

?
(a)

быть положительными и

(b)

удовлетворять 0<yi<10<y_{i}<1.

Задача 1.4.12

Обобщите группу преобразований из Примера 4.10 на случай ss генеральных совокупностей {yij,j=1,…,Ni},i=1,…,s\left\{ y_{i j}, j=1, \ldots , N_{i}\right\} , i=1, \ldots , s, когда из ii-й совокупности берётся случайная выборка объёма nin_{i}.

?
Задача 1.4.13

Пусть UU — положительная случайная величина, и пусть

X=bU1/c,b>0,c>0. X=b U^{1 / c}, \quad b>0, \quad c>0.
?
(a)

Покажите, что это определяет групповое семейство.

(b)

Если UU распределена как E(0,1)E(0,1), то XX распределена согласно распределению Вейбулла с плотностью

cb(xb)c−1e−(x/b)c,x>0 \frac{c}{b}\left(\frac{x}{b}\right)^{c-1} e^{-(x / b)^{c}}, \quad x>0
Задача 1.4.14

Если FF и F0F_{0} — две непрерывные строго возрастающие функции распределения на вещественной прямой, и если функция распределения UU равна F0F_{0}, а gg строго возрастает, покажите, что функция распределения g(U)g(U) равна FF тогда и только тогда, когда g=F−1(F0)g=F^{-1}\left(F_{0}\right).

?
Задача 1.4.15

Следующие два семейства распределений не являются групповыми семействами:

?
(a)

Класс биномиальных распределений b(p,n)b(p, n) с фиксированным nn и 0<p<10<p<1.

(b)

Класс распределений Пуассона P(λ),0<λP(\lambda ), 0<\lambda.

Задача 1.4.16

Пусть X1,…,XrX_{1}, \ldots , X_{r} имеют многомерное нормальное распределение с E(Xi)=ξiE\left(X_{i}\right)=\xi_{i} и ковариационной матрицей Σ\Sigma. Если XX — вектор-столбец с элементами XiX_{i}, а BB — матрица констант размера r×rr \times r, то BXB X имеет многомерное нормальное распределение со средним BξB \xi и ковариационной матрицей BΣB′B \Sigma B^{\prime }.

?
§
Задача 1.5.1

Найдите естественное параметрическое пространство для (5.2), когда s=1,T1(x)=x,μs=1, T_{1}(x)=x, \mu — мера Лебега, а h(x)h(x) равно

  1. e−∣x∣e^{-\left|x\right|} и

  2. e−∣x∣/(1+x2)e^{-\left|x\right|} /\left(1+x^{2}\right).

?
Задача 1.5.2

Предположим, что в (5.2) s=2s=2 и T2(x)=T1(x)T_{2}(x)=T_{1}(x). Объясните, почему невозможно оценить η1\eta_{1}. [Указание: Сравните эту модель с моделью, полученной при подстановке η1′=η1+c,η2′=η2−c\eta_{1}^{\prime }=\eta_{1}+c, \eta_{2}^{\prime }=\eta_{2}-c.]

?
Задача 1.5.3

Покажите, что распределение выборки из pp-мерной нормальной плотности (4.15) образует ss-мерное экспоненциальное семейство. Найдите ss и определите функции ηi\eta_{i}, TiT_{i} и BB из (5.1).

?
Задача 1.5.4

Efron (1975) даёт весьма общие определения кривизны, обобщающие (10.1) и (10.2). Для ss-мерного семейства (5.1) с ковариационной матрицей Σθ\Sigma_{\theta }, если θ\theta — скаляр, определим статистическую кривизну как γθ=(∣Mθ∣/m113)1/2\gamma_{\theta }=\left(\left|M_{\theta }\right| / m_{11}^{3}\right)^{1 / 2}, где

Mθ=(m11m12m21m22)=(η˙θ′Σ˙θηθη˙θ′Σθη¨θη¨θ′Σθη˙θη¨θ′Σθη¨θ), M_{\theta }=\left(\begin{array}{ll} m_{11} & m_{12} \\ m_{21} & m_{22} \\ \end{array}\right)=\left(\begin{array}{ll} \dot{\eta }_{\theta }^{\prime } \dot{\Sigma }_{\theta } \eta _{\theta } & \dot{\eta }_{\theta }^{\prime } \Sigma _{\theta } \ddot{\eta }_{\theta } \\ \ddot{\eta }_{\theta }^{\prime } \Sigma _{\theta } \dot{\eta }_{\theta } & \ddot{\eta }_{\theta }^{\prime } \Sigma _{\theta } \ddot{\eta }_{\theta } \end{array}\right),

причём η(θ)={ηi(θ)},η˙(θ)={ηi′(θ)}\eta (\theta )=\left\{ \eta_{i}(\theta )\right\} , \dot{\eta }(\theta )=\left\{ \eta_{i}^{\prime }(\theta )\right\} и η¨(θ)={ηi′′(θ)}\ddot{\eta }(\theta )=\left\{ \eta_{i}^{\prime \prime }(\theta )\right\}. Вычислите кривизну семейства (см. Пример 6.19) Cexp⁡[−∑i=1n(xi−θ)m]C \exp \left[-\sum_{i=1}^{n}\left(x_{i}-\theta \right)^{m}\right] для m=2,3,4m=2,3,4. Упорядочены ли значения γθ\gamma_{\theta } так, как вы ожидали?

?
Задача 1.5.5

Пусть (X1,X2)\left(X_{1}, X_{2}\right) имеет двумерное нормальное распределение с вектором средних ξ=(ξ1,ξ2)\xi =\left(\xi_{1}, \xi_{2}\right) и единичной ковариационной матрицей. В каждой из следующих ситуаций проверьте кривизну γθ\gamma_{\theta } семейства.

?
(a)

ξ=(θ,θ),γθ=0\xi =(\theta , \theta ), \gamma_{\theta }=0.

(b)

ξ=(θ1,θ2),θ12+θ22=r2,γθ=1/r\xi =\left(\theta_{1}, \theta_{2}\right), \theta_{1}^{2}+\theta_{2}^{2}=r^{2}, \gamma_{\theta }=1 / r.

Задача 1.5.6

В плотности (5.1)

?
(a)

Для s=1s=1 покажите, что Eθ[T(X)]=B′(θ)/η′(θ)E_{\theta }[T(X)]=B^{\prime }(\theta ) / \eta^{\prime }(\theta ) и var⁡θ[T(X)]=B′′(θ)[η′(θ)]2−η′′(θ)B′(θ)[η′(θ)]3\operatorname {var}_{\theta }[T(X)]=\frac{B^{\prime \prime }(\theta )}{\left[\eta^{\prime }(\theta )\right]^{2}}-\frac{\eta^{\prime \prime }(\theta ) B^{\prime }(\theta )}{\left[\eta^{\prime }(\theta )\right]^{3}}.

(b)

Для s>1s>1 покажите, что Eθ[T(X)]=J−1∇BE_{\theta }[T(X)]=J^{-1} \nabla B, где JJ — матрица Якоби, определённая как J={∂ηj∂θi}J=\left\{ \frac{\partial \eta_{j}}{\partial \theta_{i}}\right\}, а ∇B\nabla B — вектор градиента ∇B={∂∂θiB(θ)}\nabla B=\left\{ \frac{\partial }{\partial \theta_{i}} B(\theta )\right\}.

(Общее рассмотрение этих тождеств см. в Johnson, Ladalla, and Liu (1979).)

Задача 1.5.7

Проверьте соотношения

?
(a)

(5.22) и

(b)

(5.26).

Задача 1.5.8

Для биномиального распределения (5.28) проверьте

?
(a)

производящую функцию моментов (5.30) и

(b)

моменты (5.31).

Задача 1.5.9

Для распределения Пуассона (5.32) проверьте моменты (5.35).

?
Задача 1.5.10

В последовательности испытаний Бернулли с вероятностью успеха pp пусть X+mX+m — число испытаний, необходимых для достижения mm успехов.

?
(a)

Покажите, что распределение XX — отрицательное биномиальное распределение — задаётся Таблицей 5.1.

(b)

Проверьте, что вероятности отрицательного биномиального распределения в сумме дают 1, разложив (1p−qp)−m=pm(1−q)−m\left(\frac{1}{p}-\frac{q}{p}\right)^{-m} =p^{m}(1-q)^{-m}.

(c)

Покажите, что распределения из пункта (a) образуют однопараметрическое экспоненциальное семейство.

(d)

Покажите, что производящая функция моментов XX равна MX(u)=pm/(1−qeu)mM_{X}(u)=p^{m} /\left(1-q e^{u}\right)^{m}.

(e)

Покажите, что E(X)=mq/pE(X)=m q / p и var⁡(X)=mq/p2\operatorname {var}(X)=m q / p^{2}.

(f)

Разложив KX(u)K_{X}(u), покажите, что первые четыре кумулянта XX равны k1=mq/pk_{1}=m q / p, k2=mq/p2,k3=mq(1+q)/p3k_{2}=m q / p^{2}, k_{3}=m q(1+q) / p^{3} и k4=mq(1+4q+q2)/p4k_{4}=m q\left(1+4 q+q^{2}\right) / p^{4}.

Задача 1.5.11

В предыдущей задаче пусть Xi+1X_{i}+1 — число испытаний, требуемых после (i−1)(i-1)-го успеха до наступления следующего успеха. Используя тот факт, что X=Σi=1mXiX=\Sigma_{i=1}^{m} X_{i}, найдите альтернативный вывод среднего и дисперсии из пункта (e).

?
Задача 1.5.12

Дискретная случайная величина с вероятностями

P(X=x)=a(x)θx/C(θ),x=0,1,…;a(x)≥0;θ>0, P(X=x)=a(x) \theta ^{x} / C(\theta ), \quad x=0,1, \ldots ; a(x) \geq 0 ; \theta >0,

называется распределением степенного ряда. Это экспоненциальное семейство (5.1) с s=1,η=log⁡θs=1, \eta =\log \theta и T=XT=X. Производящая функция моментов равна MX(u)=C(θeu)/C(θ)M_{X}(u)=C\left(\theta e^{u}\right) / C(\theta ).

?
Задача 1.5.13

Покажите, что биномиальное, отрицательное биномиальное распределения и распределение Пуассона являются частными случаями распределения степенного ряда из Задачи 5.12, и найдите θ\theta и C(θ)C(\theta ).

?
Задача 1.5.14

Распределение из Задачи 5.12 с a(x)=1/xa(x)=1 / x и C(θ)=−log⁡(1−θ),x=1,2,…;0<θ<1C(\theta )=-\log (1-\theta ), x= 1,2, \ldots ; 0<\theta <1, является логарифмическим распределением. Покажите, что производящая функция моментов равна log⁡(1−θeu)/log⁡(1−θ)\log \left(1-\theta e^{u}\right) / \log (1-\theta ), и найдите E(X)E(X) и var⁡(X)\operatorname {var}(X).

?
Задача 1.5.15

Для полиномиального распределения (5.4) проверьте формулы моментов (5.16).

?
Задача 1.5.16

В качестве альтернативы использованию (5.14) и (5.15), получите моменты (5.16), представив каждую XiX_{i} в виде суммы nn индикаторов, как это было сделано в (5.5):

?
Задача 1.5.17

Для гамма-распределения (5.41).

?
(a)

проверьте формулы (5.42), (5.43) и (5.44);

(b)

покажите, что (5.43) без среднего члена выполняется не только для всех положительных целых rr, но и для всех действительных r>−αr>-\alpha.

Задача 1.5.18

Лемма 5.15 (тождество Стейна) утверждает: если XX распределена с плотностью (5.2) и gg — произвольная дифференцируемая функция, такая что E∣g′(X)∣<∞E\left|g'(X)\right| < \infty, то

E[(h′(X)h(X)+∑i=1sηiTi′(X))g(X)]=−Eg′(X), E\left[\left(\frac{h'(X)}{h(X)}+\sum _{i=1}^{s} \eta _{i} T_{i}'(X)\right) g(X)\right]=-E g'(X),

при условии, что носитель XX — это (−∞,∞)(-\infty , \infty ). Если носитель XX — ограниченный интервал (a,b)(a, b), то тождество выполняется, если exp⁡{∑iηiTi(x)}h(x)→0\exp \left\{ \sum_{i} \eta_{i} T_{i}(x)\right\} h(x) \rightarrow 0 при x→ax \rightarrow a или bb.

?
(a)

Докажите Лемму 5.15. (Используйте интегрирование по частям.)

(b)

Выбирая g(x)g(x) равным x2x^{2} и x3x^{3}, используйте тождество Стейна для вычисления третьего и четвёртого моментов распределения N(μ,σ2)N\left(\mu , \sigma^{2}\right).

Задача 1.5.19

Используя Лемму 5.15:

?
(a)

Выведите вид тождества для X∼Gamma⁡(α,b)X \sim \operatorname {Gamma}(\alpha , b) и используйте его для проверки моментов, приведённых в (5.44).

(b)

Выведите вид тождества для X∼Beta⁡(a,b)X \sim \operatorname {Beta}(a, b) и используйте его для проверки того, что E(X)=a/(a+b)E(X)= a /(a+b) и var⁡(X)=ab/(a+b)2(a+b+1)\operatorname {var}(X)=a b /(a+b)^{2}(a+b+1).

Задача 1.5.20

В качестве альтернативы подходу Задачи 5.19(b) для вычисления моментов X∼B(a,b)X \sim B(a, b) можно вывести общую формулу для EXkE X^{k} (аналогичную уравнению (5.43)). Сделайте это и используйте полученный результат для проверки среднего и дисперсии XX, найденных в Задаче 5.19. [Указание: Запишите EXkE X^{k} в виде интеграла от xc−1(1−x)d−1x^{c-1}(1-x)^{d-1} и используйте константу B(c,d)B(c, d) из Таблицы 5.1. Заметьте, что аналогичный подход применим и ко многим другим распределениям, включая распределения χ2\chi^{2}, Стьюдента tt и FF.]

?
Задача 1.5.21

Тождество Стейна также может быть применено к дискретным экспоненциальным семействам, как показано у Hudson (1978) и обобщено у Hwang (1982a). Если XX принимает значения в N={0,1,…}N=\left\{ 0,1, \ldots \right\} с функцией вероятности

pθ(x)=exp⁡[θx−B(θ)]h(x), p_{\theta }(x)=\exp [\theta x-B(\theta )] h(x),

то для любой g:N→ℜg: N \rightarrow \Re с Eθ∣g(X)∣<∞E_{\theta }\left|g(X)\right|<\infty имеет место тождество

Eg(X)=e−θE{t(X)g(X−1)} E g(X)=e^{-\theta } E\left\{ t(X) g(X-1)\right\}

где t(0)=0t(0)=0 и t(x)=h(x−1)/h(x)t(x)=h(x-1) / h(x) при x>0x>0.

?
(a)

Докажите это тождество.

(b)

Используйте тождество для вычисления первых четырёх моментов биномиального распределения (5.31).

(c)

Используйте тождество для вычисления первых четырёх моментов распределения Пуассона (5.35).

Задача 1.5.22

Обратное гауссовское распределение IG(λ,μ)I G(\lambda , \mu ) имеет функцию плотности

λ2πe(λμ)1/2x−3/2e−12(λx+μx),x>0,λ,μ>0. \sqrt{\frac{\lambda }{2 \pi }} e^{(\lambda \mu )^{1 / 2}} x^{-3 / 2} e^{-\frac{1}{2}\left(\frac{\lambda }{x}+\mu x\right)}, \quad x>0, \lambda , \mu >0.
?
(a)

Покажите, что эта плотность образует экспоненциальное семейство.

(b)

Покажите, что эта плотность является семейством масштаба (в смысле определения из Примера 4.1).

(c)

Покажите, что статистики Xˉ=(1/n)Σxi\bar{X}=(1 / n) \Sigma x_{i} и S∗=Σ(1/xi−1/xˉ)S^{*}=\Sigma \left(1 / x_{i}-1 / \bar{x}\right) являются полными достаточными статистиками.

(d)

Покажите, что Xˉ∼IG(nλ,nμ)\bar{X} \sim I G(n \lambda , n \mu ) и S∗∼(1/λ)χn−12S^{*} \sim (1 / \lambda ) \chi_{n-1}^{2}.

Замечание: Наряду с нормальным и гамма-распределением, обратное гауссовское распределение завершает тройку семейств, являющихся одновременно экспоненциальным и групповым семейством распределений. Этот факт играет важную роль в теории распределений, основанной на седловых приближениях (Daniels 1983) или теории правдоподобия (Barndorff-Nielsen 1983).

Задача 1.5.23

В Примере 5.14 покажите, что

?
(a)

χ12\chi_{1}^{2} — это распределение Y2Y^{2}, где YY распределена как N(0,1)N(0,1);

(b)

χn2\chi_{n}^{2} — это распределение Y12+…+Yn2Y_{1}^{2}+\ldots +Y_{n}^{2}, где YiY_{i} независимы и распределены как N(0,1)N(0,1).

Задача 1.5.24

Найдите значения α\alpha, при которых плотность (5.41) является

?
(a)

убывающей функцией xx на (0,∞)(0, \infty ) и

(b)

возрастает при x<x0x<x_{0} и убывает при x>x0(0<x0)x>x_{0}\left(0<x_{0}\right).

В случае (b) найдите моду плотности.

Задача 1.5.25

Случайная величина XX имеет распределение Парето P(c,k)P(c, k), если её функция распределения равна 1−(k/x)c1-(k / x)^{c}, x>k>0,c>0x>k>0, c>0.

?
(a)

Распределения P(c,1)P(c, 1) образуют однопараметрическое экспоненциальное семейство (5.2) с η=−c\eta =-c и T=log⁡XT=\log X.

(b)

Статистика TT распределена как E(log⁡k,1/c)E(\log k, 1 / c).

(c)

Семейство P(c,k)(0<k,0<c)P(c, k)(0<k, 0<c) является групповым семейством.

Задача 1.5.26

Если (X,Y)(X, Y) распределена согласно двумерному нормальному распределению (4.16) с ξ=η=0\xi =\eta =0:

?
(a)

Покажите, что производящая функция моментов (X,Y)(X, Y) равна

MX,Y(u1,u2)=e−[u12σ2+2ρστu1u2+u22τ2]/2. M_{X, Y}\left(u_{1}, u_{2}\right)=e^{-\left[u_{1}^{2} \sigma ^{2}+2 \rho \sigma \tau u_{1} u_{2}+u_{2}^{2} \tau ^{2}\right] / 2}.
(b)

Используя (a), покажите, что

μ12=μ21=0,μ11=ρστ,μ13=3ρστ3,μ31=3ρσ3τ,μ22=(1+2ρ2)σ2τ2. \begin{align} & \mu _{12}=\mu _{21}=0, \quad \mu _{11}=\rho \sigma \tau , \\ & \mu _{13}=3 \rho \sigma \tau ^{3}, \quad \mu _{31}=3 \rho \sigma ^{3} \tau , \quad \mu _{22}=\left(1+2 \rho ^{2}\right) \sigma ^{2} \tau ^{2}. \end{align}
Задача 1.5.27
?
(a)

Если X — случайный вектор-столбец с математическим ожиданием ξ\xi, то ковариационная матрица X\mathbf{X} равна cov⁡(X)=E[(X′−ξ)(X′−ξ′)]\operatorname {cov}(\mathbf{X})=E\left[\left(\mathbf{X}^{\prime }-\xi \right)\left(\mathbf{X}^{\prime }-\xi^{\prime }\right)\right].

(b)

Если плотность X равна (4.15), то ξ=a\xi =\mathbf{a} и cov⁡(X)=Σ\operatorname {cov}(\mathbf{X})=\Sigma.

Задача 1.5.28
?
(a)

Пусть XX распределена с плотностью pθ(x)p_{\theta }(x), заданной (5.1), и пусть AA — произвольное фиксированное подмножество выборочного пространства. Тогда распределения XX, усечённые на AA, то есть распределения с плотностью pθ(x)IA(x)/Pθ(A)p_{\theta }(x) I_{A}(x) / P_{\theta }(A), снова образуют экспоненциальное семейство.

(b)

Приведите пример, в котором естественное параметрическое пространство исходного экспоненциального семейства является собственным подмножеством естественного параметрического пространства усечённого семейства.

Задача 1.5.29

Если XiX_{i} независимо распределены согласно Γ(αi,b)\Gamma \left(\alpha_{i}, b\right), покажите, что ΣXi\Sigma X_{i} распределена как Γ(Σαi,b)\Gamma \left(\Sigma \alpha_{i}, b\right). [Указание: Метод 1. Докажите это сначала для суммы двух гамма-величин с помощью перехода к новым переменным Y1=X1+X2,Y2=X1/X2Y_{1}=X_{1}+X_{2}, Y_{2}=X_{1} / X_{2}, а затем используйте индукцию. Метод 2. Найдите производящую функцию моментов ΣXi\Sigma X_{i} и используйте тот факт, что распределение однозначно определяется своей производящей функцией моментов, если последняя существует хотя бы для некоторого u≠0u \neq 0.]

?
Задача 1.5.30

Если XiX_{i} независимо распределены согласно распределениям Пуассона P(λi)P\left(\lambda_{i}\right), найдите распределение ΣXi\Sigma X_{i}.

?
Задача 1.5.31

Пусть X1,…,XnX_{1}, \ldots , X_{n} независимо распределены как Γ(α,b)\Gamma (\alpha , b). Покажите, что совместное распределение является двухпараметрическим экспоненциальным семейством, и определите функции ηi,Ti\eta_{i}, T_{i} и BB из (5.1).

?
Задача 1.5.32

Если YY распределена как Γ(α,b)\Gamma (\alpha , b), найдите распределение clog⁡Yc \log Y и покажите, что при фиксированном α\alpha и меняющемся bb оно определяет экспоненциальное семейство.

?
Задача 1.5.33

Morris (1982, 1983b) исследовал свойства естественных экспоненциальных семейств с квадратичными функциями дисперсии. Существует лишь шесть таких семейств: нормальное, биномиальное, гамма-распределение, распределение Пуассона, отрицательное биномиальное и менее известное обобщённое распределение гиперболического секанса, являющееся плотностью X=1πlog⁡(Y1−Y)X=\frac{1}{\pi } \log \left(\frac{Y}{1-Y}\right), где Y∼Beta⁡(12+θπ,12−θπ),∣θ∣<π2Y \sim \operatorname {Beta}\left(\frac{1}{2}+\frac{\theta }{\pi }, \frac{1}{2}-\frac{\theta }{\pi }\right),\left|\theta \right|<\frac{\pi }{2}.

?
(a)

Найдите плотность XX и покажите, что она образует экспоненциальное семейство.

(b)

Найдите среднее и дисперсию XX и покажите, что дисперсия равна 1+μ21+\mu^{2}, где μ\mu — среднее.

Дальнейшие работы по квадратичным и другим семействам со степенной дисперсией принадлежат Bar-Lev and Enis (1986, 1988), Bar-Lev and Bshouty (1989) и Letac and Mora (1990).

§
Задача 1.6.1

Обобщите Пример 6.2 на случай, когда X1,…,XrX_{1}, \ldots , X_{r} независимо распределены с распределениями Пуассона P(λi)P\left(\lambda_{i}\right), где λi=aiλ(ai>0\lambda_{i}=a_{i} \lambda \left(a_{i}>0\right., известно )).

?
Задача 1.6.2

Пусть X1,…,XnX_{1}, \ldots , X_{n} независимые одинаково распределённые согласно распределению FF и плотности вероятности ff. Покажите, что условное распределение, при условии X(i)=aX_{(i)}=a, значений в количестве i−1i-1, лежащих слева от aa, и значений в количестве n−in-i, лежащих справа от aa, совпадает с распределением i−1i-1 переменных, независимо распределённых с плотностью вероятности f(x)/F(a)f(x) / F(a), и n−in-i переменных, независимо распределённых с плотностью f(x)/[1−F(a)]f(x) /[1-F(a)], соответственно, причём эти два набора (условно) независимы друг от друга.

?
Задача 1.6.3

Пусть ff — положительная интегрируемая функция на (0,∞)(0, \infty ), и пусть pθ(x)p_{\theta }(x) — плотность на (0,θ)(0, \theta ), определённая как pθ(x)=c(θ)f(x)p_{\theta }(x)=c(\theta ) f(x), если 0<x<θ0<x<\theta, и 0 в противном случае. Если X1,…,XnX_{1}, \ldots , X_{n} независимые одинаково распределённые с плотностью pθp_{\theta }, покажите, что X(n)X_{(n)} достаточна для θ\theta.

?
Задача 1.6.4

Пусть ff — положительная интегрируемая функция, определённая на (−∞,∞)(-\infty , \infty ), и пусть pξ,η(x)p_{\xi , \eta }(x) — плотность вероятности, определённая как pξ,η(x)=c(ξ,η)f(x)p_{\xi , \eta }(x)=c(\xi , \eta ) f(x), если ξ<x<η\xi <x<\eta, и 0 в противном случае. Если X1,…,XnX_{1}, \ldots , X_{n} независимые одинаково распределённые с плотностью pξ,ηp_{\xi , \eta }, покажите, что (X(1),X(n))\left(X_{(1)}, X_{(n)}\right) достаточна для (ξ,n)(\xi , n).

?
Задача 1.6.5

Покажите, что каждая из статистик T1−T4T_{1}-T_{4} Примера 6.11 является достаточной.

?
Задача 1.6.6

Теорема 6.5 (критерий факторизации) утверждает: необходимым и достаточным условием того, что статистика TT достаточна для семейства P={Pθ,θ∈Ω}\mathcal{P}=\left\{ P_{\theta }, \theta \in \Omega \right\} распределений XX, доминируемого σ\sigma-конечной мерой μ\mu, является существование неотрицательных функций gθg_{\theta } и hh, таких что плотности pθp_{\theta } распределений PθP_{\theta } удовлетворяют pθ(x)=gθ[T(x)]h(x)p_{\theta }(x)=g_{\theta }[T(x)] h(x) (п.в. μ\mu).

Следствие 6.13 утверждает: в условиях Теоремы 6.5 необходимым и достаточным условием достаточности статистики UU является то, что для любых фиксированных θ\theta и θ0\theta_{0} отношение pθ(x)/pθ0(x)p_{\theta }(x) / p_{\theta_{0}}(x) является функцией только от U(x)U(x).

Докажите Следствие 6.13.

?
Задача 1.6.7

Пусть X1,…,XmX_{1}, \ldots , X_{m} и Y1,…,YnY_{1}, \ldots , Y_{n} независимо распределены согласно N(ξ,σ2)N\left(\xi , \sigma^{2}\right) и N(η,τ2)N\left(\eta , \tau^{2}\right) соответственно. Найдите минимальные достаточные статистики для следующих случаев:

?
(a)

ξ,η,σ,τ\xi , \eta , \sigma , \tau произвольны: −∞<ξ,η<∞,0<σ,τ-\infty <\xi , \eta <\infty , 0<\sigma , \tau.

(b)

σ=τ\sigma =\tau, а ξ,η,σ\xi , \eta , \sigma произвольны.

(c)

ξ=η\xi =\eta, а ξ,σ,τ\xi , \sigma , \tau произвольны.

Задача 1.6.8

Пусть X1,…,XnX_{1}, \ldots , X_{n} независимые одинаково распределённые согласно N(σ,σ2),0<σN\left(\sigma , \sigma^{2}\right), 0<\sigma. Найдите минимальный набор достаточных статистик.

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

Если (x1,…,xn)\left(x_{1}, \ldots , x_{n}\right) и (y1,…,yn)\left(y_{1}, \ldots , y_{n}\right) имеют одинаковые элементарные симметрические функции Σxi=Σyi,Σi≠jxiyj=Σi≠jyiyj,…,x1⋯xn=y1⋯yn\Sigma x_{i}=\Sigma y_{i}, \Sigma_{i \neq j} x_{i} y_{j}=\Sigma_{i \neq j} y_{i} y_{j}, \ldots , x_{1} \cdots x_{n}=y_{1} \cdots y_{n}, то yy являются перестановкой xx.

(b)

В обозначениях Примера 6.10 покажите, что UU эквивалентна VV. [Указание: Сравните коэффициенты и корни многочленов P(x)=Π(x−ui)P(x)=\Pi \left(x-u_{i}\right) и Q(x)=Π(x−vi)Q(x)= \Pi \left(x-v_{i}\right).]

Задача 1.6.10

Покажите, что порядковые статистики являются минимальными достаточными для семейства сдвига (6.7), когда ff — плотность

?
(a)

двойного экспоненциального распределения D(0,1)D(0,1).

(b)

распределения Коши C(0,1)C(0,1).

Задача 1.6.11

Докажите следующее обобщение Теоремы 6.12 на семейства без общего носителя.

Теорема 9.1. Пусть P\mathcal{P} — конечное семейство с плотностями pi,i=0,…,kp_{i}, i=0, \ldots , k, и для любого xx пусть S(x)S(x) — множество пар индексов (i,j)(i, j), для которых pi(x)+pj(x)>0p_{i}(x)+p_{j}(x)>0. Тогда статистика

T(X)={pj(X)pi(X),i<j и (i,j)∈S(X)} T(X)=\left\{ \frac{p_{j}(X)}{p_{i}(X)}, \quad i<j \text{ и }(i, j) \in S(X)\right\}

является минимальной достаточной. Здесь pj(x)/pi(x)=∞p_{j}(x) / p_{i}(x)=\infty, если pi(x)=0p_{i}(x)=0 и pj(x)>0p_{j}(x)>0.

?
Задача 1.6.12

В Задаче 6.11 недостаточно заменить pi(X)p_{i}(X) на p0(X)p_{0}(X). Чтобы убедиться в этом, положим k=2k=2 и p0=U(−1,0),p1=U(0,1)p_{0}=U(-1,0), p_{1}=U(0,1), и p2(x)=2x,0<x<1p_{2}(x)=2 x, 0<x<1.

?
Задача 1.6.13

Пусть k=1k=1 и Pi=U(i,i+1),i=0,1P_{i}=U(i, i+1), i=0,1.

?
(a)

Покажите, что минимальная достаточная статистика для P={P0,P1}\mathcal{P}=\left\{ P_{0}, P_{1}\right\} есть T(X)=iT(X)=i, если i<X<i+1i<X<i+1, i=0,1i=0,1.

(b)

Пусть X1X_{1} и X2X_{2} независимые одинаково распределённые согласно распределению из P\mathcal{P}. Покажите, что каждая из двух статистик T1=T(X1)T_{1}=T\left(X_{1}\right) и T2=T(X2)T_{2}=T\left(X_{2}\right) достаточна для (X1,X2)\left(X_{1}, X_{2}\right).

(c)

Покажите, что T(X1)T\left(X_{1}\right) и T(X2)T\left(X_{2}\right) эквивалентны.

Задача 1.6.14

В Лемме 6.14 покажите, что предположение об общем носителе можно заменить более слабым предположением о том, что каждое P0\mathcal{P}_{0}-нулевое множество также является P\mathcal{P}-нулевым множеством, так что (п.в. P0\mathcal{P}_{0}) эквивалентно (п.в. P\mathcal{P}).

?
Задача 1.6.15

Пусть X1,…,XnX_{1}, \ldots , X_{n} независимые одинаково распределённые согласно распределению из P={U(0,θ),θ>0}\mathcal{P}=\left\{ U(0, \theta ), \theta >0\right\}, и пусть P0\mathcal{P}_{0} — подсемейство P\mathcal{P}, для которого θ\theta рационально. Покажите, что каждое P0\mathcal{P}_{0}-нулевое множество в выборочном пространстве также является P\mathcal{P}-нулевым множеством.

?
Задача 1.6.16

Пусть X1,…,XnX_{1}, \ldots , X_{n} независимые одинаково распределённые согласно распределению из семейства P\mathcal{P}. Покажите, что TT является минимальной достаточной в следующих случаях:

?
(a)

P={U(0,θ),θ>0};T=X(n)\mathcal{P}=\left\{ U(0, \theta ), \theta >0\right\} ; T=X_{(n)}.

(b)

P={U(θ1,θ2),−∞<θ1<θ2<∞};T=(X(1),X(n))\mathcal{P}=\left\{ U\left(\theta_{1}, \theta_{2}\right),-\infty <\theta_{1}<\theta_{2}<\infty \right\} ; T=\left(X_{(1)}, X_{(n)}\right).

(c)

P={U(θ−1/2,θ+1/2),−∞<θ<∞};T=(X(1),X(n))\mathcal{P}=\left\{ U(\theta -1 / 2, \theta +1 / 2),-\infty <\theta <\infty \right\} ; T=\left(X_{(1)}, X_{(n)}\right).

Задача 1.6.17

Решите предыдущую задачу для следующих случаев:

?
(a)

P={E(θ,1),−∞<θ<∞};T=X(1)\mathcal{P}=\left\{ E(\theta , 1),-\infty <\theta <\infty \right\} ; T=X_{(1)}.

(b)

P={E(0,b),0<b};T=ΣXi\mathcal{P}=\left\{ E(0, b), 0<b\right\} ; T=\Sigma X_{i}.

(c)

P={E(a,b),−∞<a<∞,0<b};T=(X(1),Σ[Xi−X(1)])\mathcal{P}=\left\{ E(a, b),-\infty <a<\infty , 0<b\right\} ; T=\left(X_{(1)}, \Sigma \left[X_{i}-X_{(1)}\right]\right).

Задача 1.6.18

Покажите, что статистики X(1)X_{(1)} и Σ[Xi−X(1)]\Sigma \left[X_{i}-X_{(1)}\right] из Задачи 6.17(c) независимо распределены как E(a,b/n)E(a, b / n) и bb Gamma (n−2,1)(n-2,1) соответственно.

?
Задача 1.6.19

Покажите, что достаточные статистики

  1. Задачи 6.3 и

  2. Задачи 6.4

являются минимальными достаточными.

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

Покажите, что в искривлённом экспоненциальном семействе N(θ,θ)N(\theta , \theta ) достаточная статистика T=(∑xi,∑xi2)T= (\sum x_{i}, \sum x_{i}^{2}) не является минимальной.

(b)

Для плотности из Примера 6.19 покажите, что T=(∑xi,∑xi2,∑xi3)T=\left(\sum x_{i}, \sum x_{i}^{2}, \sum x_{i}^{3}\right) является минимальной достаточной статистикой.

Задача 1.6.21

Для ситуации из Примера 6.25(ii) найдите несмещённую оценку ξ\xi, основанную на ∑Xi\sum X_{i}, и ещё одну, основанную на ∑Xi2\sum X_{i}^{2}; отсюда сделайте вывод, что T=(∑Xi,∑Xi2)T=\left(\sum X_{i}, \sum X_{i}^{2}\right) не является полной.

?
Задача 1.6.22

Для ситуации из Примера 6.26 покажите, что X\mathbf{X} является минимальной достаточной и полной.

?
Задача 1.6.23

Для ситуации из Примера 6.27:

?
(a)

Покажите, что X=(X1,X2)\mathbf{X}=\left(X_{1}, X_{2}\right) является минимальной достаточной для семейства (6.16) с ограничением (6.17).

(b)

Установите (6.18) и, следовательно, покажите, что минимальная достаточная статистика из пункта (a) не является полной.

Задача 1.6.24

Покажите, что для общей модели «доза-эффект»

pθ(x)=∏i=1m(nixi)[ηθ(di)]xi[1−ηθ(di)]ni−xi, p_{\theta }(\mathbf{x})=\prod _{i=1}^{m}\binom {n_{i}}{x_{i}}\left[\eta _{\theta }\left(d_{i}\right)\right]^{x_{i}}\left[1-\eta _{\theta }\left(d_{i}\right)\right]^{n_{i}-x_{i}},

статистика X=(X1,X2,…,Xm)\mathbf{X}=\left(X_{1}, X_{2}, \ldots , X_{m}\right) является минимальной достаточной, если существуют векторы θ1,θ2,⋯ ,θm)\left.\theta_{1}, \theta_{2}, \cdots , \theta_{m}\right), такие что матрица размера m×mm \times m

P={log⁡(ηθj(di)[1−ηθ0(di)]ηθ0(di)[1−ηθj(di)])} P=\left\{ \log \left(\frac{\eta _{\theta _{j}}\left(d_{i}\right)\left[1-\eta _{\theta _{0}}\left(d_{i}\right)\right]}{\eta _{\theta _{0}}\left(d_{i}\right)\left[1-\eta _{\theta _{j}}\left(d_{i}\right)\right]}\right)\right\}

обратима. (Указание: Теорема 6.12.)

?
Задача 1.6.25

Пусть (Xi,Yi),i=1,…,n\left(X_{i}, Y_{i}\right), i=1, \ldots , n, независимые одинаково распределённые согласно равномерному распределению на множестве RR на плоскости (x,y)(x, y), и пусть P\mathcal{P} — семейство распределений, получаемое при изменении RR в некотором классе R\mathcal{R} множеств RR. Найдите минимальную достаточную статистику для следующих случаев:

?
(a)

R\mathcal{R} — множество всех прямоугольников a1<x<a2,b1<y<b2,−∞<a1<a2<∞a_{1}<x<a_{2}, b_{1}<y<b_{2},-\infty <a_{1}<a_{2}<\infty, −∞<b1<b2<∞-\infty <b_{1}<b_{2}<\infty.

(b)

R′\mathcal{R}^{\prime } — подмножество R\mathcal{R}, для которого a2−a1=b2−b1a_{2}-a_{1}=b_{2}-b_{1}.

(c)

R′′\mathcal{R}^{\prime \prime } — подмножество R′\mathcal{R}^{\prime }, для которого a2−a1=b2−b1=1a_{2}-a_{1}=b_{2}-b_{1}=1.

Задача 1.6.26

Решите предыдущую задачу, если

?
(a)

R\mathcal{R} — множество всех треугольников со сторонами, параллельными оси xx, оси yy и прямой y=xy=x соответственно.

(b)

R′\mathcal{R}^{\prime } — подмножество R\mathcal{R}, в котором стороны, параллельные осям xx и yy, равны.

Задача 1.6.27

Сформулируйте общий результат, частными случаями которого являются Задачи 6.25(a) и 6.26(a).

?
Задача 1.6.28

Если YY распределена как E(η,1)E(\eta , 1), то распределение X=e−YX=e^{-Y} есть U(0,e−η)U\left(0, e^{-\eta }\right). (Этот результат полезен при компьютерной генерации случайных величин; см. Задачу 4.4.14.)

?
Задача 1.6.29

Если минимальная достаточная статистика существует, то необходимым условием полноты достаточной статистики является её минимальность. [Указание: Предположим, что T=h(U)T=h(U) является минимальной достаточной, а UU полна. Чтобы показать, что UU эквивалентна TT, заметьте, что в противном случае существует ψ\psi, такая что ψ(U)≠η[h(U)]\psi (U) \neq \eta [h(U)] с положительной вероятностью, где η(t)=E[ψ(U)∣t]\eta (t)=E[\psi (U) \mid t].]

?
Задача 1.6.30

Покажите, что минимальные достаточные статистики T=(X(1),X(n))T=\left(X_{(1)}, X_{(n)}\right) из Задачи 6.16(b) полны. [Указание: Используйте подход Примера 6.24.]

?
Задача 1.6.31

Для каждой из следующих задач определите, является ли минимальная достаточная статистика полной:

?
(a)

Задача 6.7(a)-(c);

(b)

Задача 6.25(a)-(c);

(c)

Задача 6.26(a) и (b).

Задача 1.6.32
?
(a)

Покажите, что если P0,P1\mathcal{P}_{0}, \mathcal{P}_{1} — два семейства распределений, такие что P0∈P1\mathcal{P}_{0} \in \mathcal{P}_{1} и каждое нулевое множество P0\mathcal{P}_{0} также является нулевым множеством P1\mathcal{P}_{1}, то достаточная статистика TT, полная для P0\mathcal{P}_{0}, также полна для P1\mathcal{P}_{1}.

(b)

Пусть P0\mathcal{P}_{0} — класс биномиальных распределений b(p,n),0<p<1,n=b(p, n), 0<p<1, n= фиксировано, и пусть P1=P0∪{Q}\mathcal{P}_{1}=\mathcal{P}_{0} \cup \left\{ Q\right\}, где QQ — распределение Пуассона с математическим ожиданием 1. Тогда P0\mathcal{P}_{0} полно, а P1\mathcal{P}_{1} — нет.

Задача 1.6.33

Пусть X1,…,XnX_{1}, \ldots , X_{n} независимые одинаково распределённые, каждая с плотностью f(x)f(x) (относительно меры Лебега), которая неизвестна. Покажите, что порядковые статистики полны.

?
Задача 1.6.34

Пусть X1,…,XnX_{1}, \ldots , X_{n} — независимая одинаково распределённая выборка из семейства сдвига-масштаба с функцией распределения F((x−a)/b)F((x-a) / b).

?
(a)

Если bb известно, покажите, что разности (X1−Xi)/b,i=2,…,n\left(X_{1}-X_{i}\right) / b, i=2, \ldots , n, являются вспомогательными.

(b)

Если aa известно, покажите, что отношения (X1−a)/(Xi−a),i=2,…,n\left(X_{1}-a\right) /\left(X_{i}-a\right), i=2, \ldots , n, являются вспомогательными.

(c)

Если ни aa, ни bb неизвестны, покажите, что величины (X1−Xi)/(X2−Xi),i=3,…,n\left(X_{1}-X_{i}\right) /\left(X_{2}-X_{i}\right), i= 3, \ldots , n, являются вспомогательными.

Задача 1.6.35

Используя теорему Басу, докажите независимость следующих пар статистик:

?
(a)

Xˉ\bar{X} и Σ(Xi−Xˉ)2\Sigma \left(X_{i}-\bar{X}\right)^{2}, где XX независимые одинаково распределённые как N(ξ,σ2)N\left(\xi , \sigma^{2}\right).

(b)

X(1)X_{(1)} и Σ[Xi−X(1)]\Sigma \left[X_{i}-X_{(1)}\right] в Задаче 6.18.

Задача 1.6.36
?
(a)

В условиях Задачи 6.18 отношения Zi=[X(n)−X(i)]/[X(n)−X(n−1)],i=1,…,n−2Z_{i}=\left[X_{(n)}-X_{(i)}\right] / \left[X_{(n)}-X_{(n-1)}\right], i=1, \ldots , n-2, независимы от {X(1),Σ[Xi−X(1)]}\left\{ X_{(1)}, \Sigma \left[X_{i}-X_{(1)}\right]\right\}.

(b)

В условиях Задач 6.16(b) и 6.30 отношения Zi=[X(i)−X(1)]/[X(n)−X(1)],i=2,…,n−1Z_{i}=\left[X_{(i)}-X_{(1)}\right] / \left[X_{(n)}-X_{(1)}\right], i=2, \ldots , n-1, независимы от {X(1),X(n)}\left\{ X_{(1)}, X_{(n)}\right\}.

Задача 1.6.37

В условиях Теоремы 6.5 пусть AA — произвольное фиксированное множество в выборочном пространстве, Pθ∗P_{\theta }^{*} — распределение PθP_{\theta }, усечённое на AA, и P∗={Pθ∗,θ∈Ω}\mathcal{P}^{*}=\left\{ P_{\theta }^{*}, \theta \in \Omega \right\}. Тогда докажите

?
(a)

если TT достаточна для P\mathcal{P}, то она достаточна для P∗\mathcal{P}^{*}.

(b)

если, кроме того, TT полна для P\mathcal{P}, то она также полна для P∗\mathcal{P}^{*}.

Обобщения этого результата были получены Tukey в 1940-х годах, а также Smith (1957). Аналогичная задача для наблюдений, которые цензурированы, а не усечены, рассматривается у Bhattacharyya, Johnson, and Mehrotra (1977).

Задача 1.6.38

Если X1,…,XnX_{1}, \ldots , X_{n} независимые одинаково распределённые как B(a,b)B(a, b),

?
(a)

Покажите, что [ΠXi,Π(1−Xi)]\left[\Pi X_{i}, \Pi \left(1-X_{i}\right)\right] является минимальной достаточной для (a,b)(a, b).

(b)

Найдите минимальную достаточную статистику при a=ba=b.

§
Задача 1.7.1

Проверьте выпуклость функций (i)-(vi) из Примера 7.3.

?
Задача 1.7.2

Покажите, что xpx^{p} вогнута на (0,∞)(0, \infty ), если 0<p<10<p<1.

?
Задача 1.7.3

Приведите пример, показывающий, что выпуклая функция не обязательно непрерывна на замкнутом интервале.

?
Задача 1.7.4

Если ϕ\phi выпукла на (a,b)(a, b), а ψ\psi выпукла и не убывает на области значений ϕ\phi, покажите, что функция ψ[ϕ(x)]\psi [\phi (x)] выпукла на (a,b)(a, b).

?
Задача 1.7.5

Докажите или опровергните контрпримером каждое из следующих утверждений. Если ϕ\phi выпукла на (a,b)(a, b), то выпукла и

  1. eϕ(x)e^{\phi (x)} и

  2. log⁡ϕ(x)\log \phi (x), если ϕ>0\phi >0.

?
Задача 1.7.6

Покажите, что если в (7.1) достигается равенство при некотором 0<γ<10<\gamma <1, то ϕ\phi линейна на [x,y][x, y].

?
Задача 1.7.7

Установите следующую лемму, полезную при исследовании функций риска некоторых оценок. (Дальнейшее обсуждение см. в Casella 1990).

Лемма 9.2. Пусть r:[0,∞)→[0,∞)r:[0, \infty ) \rightarrow [0, \infty ) вогнута. Тогда,

  1. r(t)r(t) не убывает и

  2. r(t)/tr(t) / t не возрастает.

?
Задача 1.7.8

Докажите неравенство Йенсена для случая, когда XX принимает значения x1,…,xnx_{1}, \ldots , x_{n} с вероятностями γ1,…,γn(Σγi=1)\gamma_{1}, \ldots , \gamma_{n}\left(\Sigma \gamma_{i}=1\right), непосредственно из (7.1) по индукции по nn.

?
Задача 1.7.9

Несколько иная форма теоремы Рао—Блэквелла, применимая только к дисперсии оценки, а не к произвольной выпуклой функции потерь, может быть установлена без неравенства Йенсена.

?
(a)

Для произвольной оценки δ(x)\delta (x) с var⁡[δ(X)]<∞\operatorname {var}[\delta (X)]<\infty и произвольной статистики TT покажите, что

var⁡[δ(X)]=var⁡[E(δ(X)∣T)]+E[var⁡(δ(X)∣T)]. \operatorname {var}[\delta (X)]=\operatorname {var}[E(\delta (X) \mid T)]+E[\operatorname {var}(\delta (X) \mid T)].
(b)

Опираясь на тождество из пункта (a), сформулируйте и докажите теорему типа Рао—Блэквелла для дисперсий.

(c)

Тождество из пункта (a) играет важную роль как в теоретической, так и в прикладной статистике. Например, объясните, как уравнение (1.2) можно интерпретировать как частный случай этого тождества.

Задача 1.7.10

Пусть UU равномерно распределена на (0,1)(0,1), и пусть FF — функция распределения на вещественной прямой.

?
(a)

Если FF непрерывна и строго возрастает, покажите, что F−1(U)F^{-1}(U) имеет функцию распределения FF.

(b)

Для произвольной FF покажите, что F−1(U)F^{-1}(U) по-прежнему имеет функцию распределения FF.

Задача 1.7.11

Покажите, что kk-мерный шар Σi=1kx12≤c\Sigma_{i=1}^{k} x_{1}^{2} \leq c выпуклый.

?
Задача 1.7.12

Покажите, что f(a)=∣x−a∣+∣y−a∣f(a)=\sqrt{\left|x-a\right|}+\sqrt{\left|y-a\right|} минимизируется при a=xa=x и a=ya=y.

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

Покажите, что ϕ(x)=eΣxi\phi (\mathbf{x})=e^{\Sigma x_{i}} выпукла, показав, что её матрица Гессе положительно полуопределена.

(b)

Покажите, что результат Задачи 7.4 остаётся справедливым, если ϕ\phi — выпуклая функция, определённая на открытом выпуклом множестве в EkE_{k}.

(c)

Используя (b), получите альтернативное доказательство результата из пункта (a).

Задача 1.7.14

Определите, являются ли следующие функции супер- или субгармоническими:

?
(a)

∑i=1kxip,p<1,xi>0\sum_{i=1}^{k} x_{i}^{p}, \quad p<1, x_{i}>0.

(b)

e−∑i=1kxi2e^{-\sum_{i=1}^{k} x_{i}^{2}}.

(c)

log⁡(∏i=1kxi)\log \left(\prod_{i=1}^{k} x_{i}\right).

Задача 1.7.15

Функция называется полунепрерывной снизу в точке yy, если f(y)≤lim inf⁡x→yf(x)f(y) \leq \liminf_{x \rightarrow y} f(x). Определение супергармоничности можно распространить с непрерывных функций на полунепрерывные снизу.

?
(a)

Покажите, что непрерывная функция полунепрерывна снизу.

(b)

Функция f(x)=I(a<x<b)f(x)=I(a<x<b) супергармонична на (−∞,∞)(-\infty , \infty ).

(c)

Для оценки dd параметра θ\theta покажите, что функция потерь

L(θ,d)={0 если ∣d−θ∣≤k2 если ∣d−θ∣>k L(\theta , d)= \begin{cases} 0 & \text{ если }\left|d-\theta \right| \leq k \\ 2 & \text{ если }\left|d-\theta \right|>k\end{cases}

субгармонична.

Задача 1.7.16
?
(a)

Если f:ℜp→ℜf: \Re^{p} \rightarrow \Re супергармонична, то φ(f(⋅))\varphi (f(\cdot )) также супергармонична, где φ:ℜ→ℜ\varphi : \Re \rightarrow \Re — дважды дифференцируемая возрастающая вогнутая функция.

(b)

Если hh супергармонична, то h∗(x)=∫g(x−y)h(y)dyh^{*}(x)=\int g(x-y) h(y) d y также супергармонична, где g(⋅)g(\cdot ) — плотность.

(c)

Если hγh_{\gamma } супергармонична, то и h∗(x)=∫hγ(x)dG(γ)h^{*}(x)=\int h_{\gamma }(x) d G(\gamma ) супергармонична, где G(γ)G(\gamma ) — функция распределения.

(Предполагается, что все необходимые интегралы существуют и что дифференцирование можно вносить под знак интеграла.)

Задача 1.7.17

Используя выпуклость функции ϕ\phi из Задачи 7.13, покажите, что естественное параметрическое пространство экспоненциального семейства (5.2) выпукло.

?
Задача 1.7.18

Покажите, что если ff определена и ограничена на (−∞,∞)(-\infty , \infty ) или (0,∞)(0, \infty ), то ff не может быть выпуклой (если только она не постоянна).

?
Задача 1.7.19

Покажите, что ϕ(x,y)=−xy\phi (x, y)=-\sqrt{x y} выпукла при x>0,y>0x>0, y>0.

?
Задача 1.7.20

Если ff и gg — вещественнозначные функции, такие что f2,g2f^{2}, g^{2} измеримы относительно σ\sigma-конечной меры μ\mu, докажите неравенство Шварца

(∫fgdμ)2≤∫f2dμ∫g2dμ \left(\int f g d \mu \right)^{2} \leq \int f^{2} d \mu \int g^{2} d \mu
?
Задача 1.7.21

Покажите, что функции потерь (7.24) непрерывно дифференцируемы.

?
Задача 1.7.22

Докажите утверждения, сделанные в Примере 7.20(i) и (ii).

?
Задача 1.7.23

Пусть ff — унимодальная плотность, симметричная относительно 0, и пусть L(θ,d)=ρ(d−θ)L(\theta , d)=\rho (d-\theta ) — функция потерь с ρ\rho, не убывающей на (0,∞)(0, \infty ) и симметричной относительно 0.

?
(a)

Функция ϕ(a)=E[ρ(X−a)]\phi (a)=E[\rho (X-a)], определённая в Теореме 7.15, принимает минимум при 0.

(b)

Если

Sa={x:[ρ(x+a)−ρ(x−a)][f(x+a)−f(x−a)]≠0}, S_{a}=\left\{ x:[\rho (x+a)-\rho (x-a)][f(x+a)-f(x-a)] \neq 0\right\} ,

то ϕ(a)\phi (a) принимает своё единственное минимальное значение при a=0a=0 тогда и только тогда, когда существует a0a_{0}, такое что ϕ(a0)<∞\phi \left(a_{0}\right)<\infty, и μ(Sa)>0\mu \left(S_{a}\right)>0 для всех aa. [Указание: Заметьте, что ϕ(0)≤1/2[ϕ(2a)+ϕ(−2a)]\phi (0) \leq 1 / 2[\phi (2 a)+\phi (-2 a)], причём строгое неравенство имеет место тогда и только тогда, когда μ(Sa)>0\mu \left(S_{a}\right)>0 для всех aa.]

Задача 1.7.24
?
(a)

Предположим, что ff и ρ\rho удовлетворяют условиям Задачи 7.23 и что ff строго убывает на [0,∞)[0, \infty ). Тогда, если ϕ(a0)<∞\phi \left(a_{0}\right)<\infty для некоторого a0a_{0}, ϕ(a)\phi (a) имеет единственный минимум в нуле, если только не существует c≤dc \leq d, такие что

ρ(0)=c и ρ(x)=d для всех x≠0. \rho (0)=c \quad \text{ и } \quad \rho (x)=d \text{ для всех } \quad x \neq 0.
(b)

Если ρ\rho симметрична относительно 0, строго возрастает на [0,∞)[0, \infty ), и ϕ(a0)<∞\phi \left(a_{0}\right)<\infty для некоторого a0a_{0}, то ϕ(a)\phi (a) имеет единственный минимум в (0) для всех симметричных унимодальных ff.

[Задачи 7.23 и 7.24 сообщены Dr. W.Y. Loh.]

Задача 1.7.25

Пусть ρ\rho — вещественнозначная функция, удовлетворяющая

0≤ρ(t)≤M<∞ и ρ(t)→M при t→±∞, 0 \leq \rho (t) \leq M<\infty \text{ и } \rho (t) \rightarrow M \text{ при } t \rightarrow \pm \infty ,

и пусть XX — случайная величина с непрерывной плотностью вероятности ff. Тогда ϕ(a)=E[ρ(X−1)]\phi (a)= E[\rho (X-1)] достигает своего минимума.

?
Задача 1.7.26

Пусть ϕ\phi — строго выпуклая функция, определённая на интервале II (конечном или бесконечном). Если существует значение a0a_{0} в II, минимизирующее ϕ(a)\phi (a), то a0a_{0} единственно.

?
Задача 1.7.27

Следствие 7.19 утверждает: в условиях Теоремы 7.15 (пусть ρ\rho — выпуклая функция, определённая на (−∞,∞)(-\infty , \infty ), а XX — случайная величина, такая что ϕ(a)=E[ρ(X−a)]\phi (a)=E[\rho (X-a)] конечна при некотором aa), предположим, что ρ\rho чётна, а XX симметрична относительно μ\mu. Тогда ϕ(a)\phi (a) достигает минимума при a=μa=\mu.

Обобщите Следствие 7.19 на случай, когда X\mathbf{X} и μ\boldsymbol {\mu } являются векторами.

?
§
Задача 1.8.1

Лемма 9.3 утверждает: достаточным условием сходимости YnY_{n} по вероятности к cc является E(Yn−c)2→0E\left(Y_{n}-c\right)^{2} \rightarrow 0.

?
(a)

Докажите неравенство Чебышёва: для любой случайной величины XX и неотрицательной функции g(⋅)g(\cdot ),

P(g(X)≥ε)≤1εEg(X) P(g(X) \geq \varepsilon ) \leq \frac{1}{\varepsilon } E g(X)

для каждого ε>0\varepsilon >0. (Во многих статистических приложениях полезно брать g(x)=(x−a)2/b2g(x)= (x-a)^{2} / b^{2} для некоторых констант aa и bb.)

(b)

Докажите Лемму 9.3. [Указание: Примените неравенство Чебышёва.]

Задача 1.8.2

Чтобы убедиться, что обратное утверждение к Теореме 8.2 не верно, пусть X1,…,XnX_{1}, \ldots , X_{n} независимые одинаково распределённые с E(Xi)=θ,var⁡(Xi)=σ2<∞E\left(X_{i}\right)=\theta , \operatorname {var}\left(X_{i}\right)=\sigma^{2}<\infty, и пусть δn=Xˉ\delta_{n}=\bar{X} с вероятностью 1−εn1-\varepsilon_{n} и δn=An\delta_{n}=A_{n} с вероятностью εn\varepsilon_{n}. Если εn\varepsilon_{n} и AnA_{n} — константы, удовлетворяющие

εn→0 и εnAn→∞, \varepsilon _{n} \rightarrow 0 \quad \text{ и } \quad \varepsilon _{n} A_{n} \rightarrow \infty ,

то δn\delta_{n} состоятельна для оценивания θ\theta, но E(δn−θ)2E\left(\delta_{n}-\theta \right)^{2} не стремится к нулю.

?
Задача 1.8.3

Предположим, что ρ(x)\rho (x) — чётная функция, не убывающая и неотрицательная при x≥0x \geq 0 и положительная при x>0x>0. Тогда E{ρ[δn−g(θ)]}→0E\left\{ \rho \left[\delta_{n}-g(\theta )\right]\right\} \rightarrow 0 для всех θ\theta влечёт, что δn\delta_{n} состоятельна для оценивания g(θ)g(\theta ).

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

Если An,BnA_{n}, B_{n} и YnY_{n} сходятся по вероятности к a,ba, b и yy соответственно, то An+BnYnA_{n}+B_{n} Y_{n} сходится по вероятности к a+bya+b y.

(b)

Если AnA_{n} принимает постоянное значение ana_{n} с вероятностью 1, и an→aa_{n} \rightarrow a, то An→aA_{n} \rightarrow a по вероятности.

Задача 1.8.5

Обращаясь к Примеру 8.4, покажите, что cnSn2→Pσ2c_{n} S_{n}^{2} \xrightarrow {P} \sigma^{2} для любой последовательности констант cn→1c_{n} \rightarrow 1. В частности, MLE σ^2=n−1nSn2\hat{\sigma }^{2}=\frac{n-1}{n} S_{n}^{2} является состоятельной оценкой σ2\sigma^{2}.

?
Задача 1.8.6

Проверьте уравнение (8.9).

?
Задача 1.8.7

Если {an}\left\{ a_{n}\right\} — последовательность действительных чисел, стремящаяся к aa, и если bn=(a1+…+an)/nb_{n}=\left(a_{1}+\ldots +a_{n}\right) / n, то bn→ab_{n} \rightarrow a.

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

Если δn\delta_{n} состоятельна для θ\theta, а gg непрерывна, то g(δn)g\left(\delta_{n}\right) состоятельна для g(θ)g(\theta ).

(b)

Пусть X1,…,XnX_{1}, \ldots , X_{n} независимые одинаково распределённые как N(θ,1)N(\theta , 1), и пусть g(θ)=0g(\theta )=0, если θ≠0\theta \neq 0, и g(0)=1g(0)=1. Найдите состоятельную оценку g(θ)g(\theta ).

Задача 1.8.9
?
(a)

В Примере 8.5 найдите cov⁡(Xi,Xj)\operatorname {cov}\left(X_{i}, X_{j}\right) для любых i≠ji \neq j.

(b)

Проверьте (8.10).

Задача 1.8.10
?
(a)

В Примере 8.5 найдите значение p1p_{1}, при котором pkp_{k} становится независимым от kk.

(b)

Если p1p_{1} имеет значение, найденное в (a), то для любых целых i1<⋯<iri_{1}<\cdots <i_{r} и kk совместное распределение Xi1,…,XirX_{i_{1}}, \ldots , X_{i_{r}} совпадает с распределением Xi1+k,…,Xir+kX_{i_{1}+k}, \ldots , X_{i_{r}+k}.

Задача 1.8.11

Предположим, что X1,…,XnX_{1}, \ldots , X_{n} имеют общее среднее ξ\xi и дисперсию σ2\sigma^{2}, и что cov⁡(Xi,Xj)=ρj−i\operatorname {cov}\left(X_{i}, X_{j}\right)= \rho_{j-i}. Для оценивания ξ\xi покажите, что:

?
(a)

Xˉ\bar{X} не состоятельна, если ρj−i=ρ≠0\rho_{j-i}=\rho \neq 0 для всех i≠ji \neq j;

(b)

Xˉ\bar{X} состоятельна, если ∣ρj−i∣≤Mγj−i\left|\rho_{j-i}\right| \leq M \gamma^{j-i} при ∣γ∣<1\left|\gamma \right|<1.

Задача 1.8.12

Предположим, что kn[δn−g(θ)]k_{n}\left[\delta_{n}-g(\theta )\right] сходится по распределению к непрерывному предельному распределению HH. Докажите, что:

?
(a)

Если kn′/kn→d≠0k_{n}^{\prime } / k_{n} \rightarrow d \neq 0 или ∞\infty, то kn′[δn−g(θ)]k_{n}^{\prime }\left[\delta_{n}-g(\theta )\right] также сходится к непрерывному предельному распределению.

(b)

Если kn′/kn→0k_{n}^{\prime } / k_{n} \rightarrow 0 или ∞\infty, то kn′[δn−g(θ)]k_{n}^{\prime }\left[\delta_{n}-g(\theta )\right] сходится по вероятности к нулю или бесконечности соответственно.

(c)

Если kn→∞k_{n} \rightarrow \infty, то δn→g(θ)\delta_{n} \rightarrow g(\theta ) по вероятности.

Задача 1.8.13

Покажите, что если Yn→cY_{n} \rightarrow c по вероятности, то она сходится по распределению к случайной величине YY, равной cc с вероятностью 1.

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

В Примере 8.7(i) и (ii) Yn→0Y_{n} \rightarrow 0 по вероятности. Покажите, что:

(b)

Если HnH_{n} обозначает функцию распределения YnY_{n} в Примере 8.7(i) и (ii), то Hn(a)→0H_{n}(a) \rightarrow 0 для всех a<0a<0 и Hn(a)→1H_{n}(a) \rightarrow 1 для всех a>0a>0.

(c)

Найдите lim⁡Hn(0)\lim H_{n}(0) для Примера 8.7(i) и (ii).

Задача 1.8.15

Если Tn>0T_{n}>0 удовлетворяет n[Tn−θ]→LN(0,τ2)\sqrt{n}\left[T_{n}-\theta \right] \xrightarrow {\mathcal{L}} N\left(0, \tau^{2}\right), найдите предельное распределение (надлежащим образом нормированных)

?
(a)

Tn\sqrt{T_{n}} и

(b)

log⁡Tn\log T_{n}.

Задача 1.8.16

Если TnT_{n} удовлетворяет n[Tn−θ]→LN(0,τ2)\sqrt{n}\left[T_{n}-\theta \right] \xrightarrow {\mathcal{L}} N\left(0, \tau^{2}\right), найдите предельное распределение (надлежащим образом нормированных)

?
(a)

Tn2T_{n}^{2},

(b)

log⁡∣Tn∣\log \left|T_{n}\right|,

(c)

1/Tn1 / T_{n}, и

(d)

eTne^{T_{n}}.

Задача 1.8.17

Дисперсионно-стабилизирующие преобразования — это преобразования, при которых полученная статистика имеет асимптотическую дисперсию, не зависящую от интересующих параметров. Для каждого из следующих случаев найдите асимптотическое распределение преобразованной статистики и покажите, что оно стабилизирует дисперсию.

?
(a)

Tn=1n∑i=1nXi,Xi∼Poisson⁡(λ),h(Tn)=TnT_{n}=\frac{1}{n} \sum_{i=1}^{n} X_{i}, X_{i} \sim \operatorname {Poisson}(\lambda ), h\left(T_{n}\right)=\sqrt{T_{n}}.

(b)

Tn=1n∑i=1nXi,Xi∼Bernoulli⁡(p),h(Tn)=arcsin⁡TnT_{n}=\frac{1}{n} \sum_{i=1}^{n} X_{i}, X_{i} \sim \operatorname {Bernoulli}(p), h\left(T_{n}\right)=\arcsin \sqrt{T_{n}}.

Задача 1.8.18
?
(a)

Функция v(⋅)v(\cdot ) является дисперсионно-стабилизирующим преобразованием, если оценка v(Tn)v\left(T_{n}\right) имеет асимптотическую дисперсию τ2(θ)[v′(θ)]2=c\tau^{2}(\theta )\left[v^{\prime }(\theta )\right]^{2}=c, где cc — константа, не зависящая от θ\theta.

(b)

Для произвольного положительного целого nn найдите дисперсионно-стабилизирующее преобразование, если τ2(θ)=θn\tau^{2}(\theta )=\theta^{n}. В частности, будьте внимательны в важном случае n=2n=2.

[Дисперсионно-стабилизирующее преобразование (если оно существует) является решением дифференциального уравнения, возникающего из приближения дельта-метода для дисперсии оценки (Теорема 8.12), и не зависит от распределения статистики (за исключением того, что распределение определяет вид дисперсии). Преобразования из пункта (b) известны как семейство степенных преобразований Бокса—Кокса и играют важную роль в прикладной статистике. Подробнее и для интересных обсуждений см. Bickel and Doksum 1981, Box and Cox 1982 и Hinkley and Runger 1984.]

Задача 1.8.19

Serfling (1980, Раздел 3.1) отмечает, что можно установить следующие вариации Теоремы 8.12. Покажите, что:

?
(a)

Если hh дифференцируема в окрестности θ\theta, а h′h^{\prime } непрерывна в θ\theta, то h′(θ)h^{\prime }(\theta ) можно заменить на h′(Tn)h^{\prime }\left(T_{n}\right), получив

n[h(Tn)−h(θ)]τh′(Tn)→LN(0,1). \sqrt{n} \frac{\left[h\left(T_{n}\right)-h(\theta )\right]}{\tau h^{\prime }\left(T_{n}\right)} \xrightarrow {\mathcal{L}} N(0,1).
(b)

Более того, если τ2\tau^{2} — непрерывная функция θ\theta, скажем, τ2(θ)\tau^{2}(\theta ), её можно заменить на τ2(Tn)\tau^{2}\left(T_{n}\right), получив

n[h(Tn)−h(θ)]τ(Tn)h′(Tn)→LN(0,1). \sqrt{n} \frac{\left[h\left(T_{n}\right)-h(\theta )\right]}{\tau \left(T_{n}\right) h^{\prime }\left(T_{n}\right)} \xrightarrow {\mathcal{L}} N(0,1).
Задача 1.8.20

Теорема 8.16 утверждает (она начинается в многомерной постановке, но её заключение одномерно): пусть (X1ν,…,Xsν)(X_{1 \nu }, \ldots , X_{s \nu }), ν=1,…,n\nu =1, \ldots , n, — nn независимых ss-наборов случайных величин с E(Xiν)=ξiE(X_{i \nu })=\xi_{i} и cov⁡(Xiν,Xjν)=σij\operatorname {cov}(X_{i \nu }, X_{j \nu })=\sigma_{i j}. Пусть Xˉi=∑νXiν/n\bar{X}_{i}=\sum_{\nu } X_{i \nu } / n, и предположим, что hh — вещественнозначная функция от ss аргументов с непрерывными первыми частными производными. Тогда,

n[h(Xˉ1,…,Xˉs)−h(ξ1,…,ξs)]→LN(0,v2), \sqrt{n}\left[h\left(\bar{X}_{1}, \ldots , \bar{X}_{s}\right)-h\left(\xi _{1}, \ldots , \xi _{s}\right)\right] \xrightarrow {\mathcal{L}} N\left(0, v^{2}\right),

где v2=∑i,jσij∂h∂ξi∂h∂ξjv^{2}=\sum_{i,j} \sigma_{i j} \dfrac {\partial h}{\partial \xi_{i}} \dfrac {\partial h}{\partial \xi_{j}}, при условии v2>0v^{2}>0.

Докажите Теорему 8.16.

?
Задача 1.8.21

Говорят, что последовательность чисел RnR_{n} есть o(1/kn)o\left(1 / k_{n}\right) при n→∞n \rightarrow \infty, если knRn→0k_{n} R_{n} \rightarrow 0, и есть O(1/kn)O\left(1 / k_{n}\right), если существуют MM и n0n_{0}, такие что ∣knRn∣<M\left|k_{n} R_{n}\right|<M для всех n>n0n>n_{0}, или, что эквивалентно, если knRnk_{n} R_{n} ограничена.

?
(a)

Если Rn=o(1/kn)R_{n}=o\left(1 / k_{n}\right), то Rn=0(1/kn)R_{n}=0\left(1 / k_{n}\right).

(b)

Rn=0(1)R_{n}=0(1) тогда и только тогда, когда RnR_{n} ограничена.

(c)

Rn=o(1)R_{n}=o(1) тогда и только тогда, когда Rn→0R_{n} \rightarrow 0.

(d)

Если RnR_{n} есть O(1/kn)O\left(1 / k_{n}\right), а kn′/knk_{n}^{\prime } / k_{n} стремится к конечному пределу, то RnR_{n} есть O(1/kn′)O\left(1 / k_{n}^{\prime }\right).

Задача 1.8.22
?
(a)

Если RnR_{n} и Rn′R_{n}^{\prime } оба есть O(1/kn)O\left(1 / k_{n}\right), то и Rn+Rn′R_{n}+R_{n}^{\prime } есть O(1/kn)O\left(1 / k_{n}\right).

(b)

Если RnR_{n} и Rn′R_{n}^{\prime } оба есть o(1/kn)o\left(1 / k_{n}\right), то и Rn+Rn′R_{n}+R_{n}^{\prime } есть o(1/kn)o\left(1 / k_{n}\right).

Задача 1.8.23

Предположим, что kn′/kn→∞k_{n}^{\prime } / k_{n} \rightarrow \infty.

?
(a)

Если Rn=0(1/kn)R_{n}=0\left(1 / k_{n}\right) и Rn′=0(1/kn′)R_{n}^{\prime }=0\left(1 / k_{n}^{\prime }\right), то Rn+Rn′=0(1/kn)R_{n}+R_{n}^{\prime }=0\left(1 / k_{n}\right).

(b)

Если Rn=o(1/kn)R_{n}=o\left(1 / k_{n}\right) и Rn′=o(1/kn′)R_{n}^{\prime }=o\left(1 / k_{n}^{\prime }\right), то Rn+Rn′=o(1/kn)R_{n}+R_{n}^{\prime }=o\left(1 / k_{n}\right).

Задача 1.8.24

Последовательность случайных величин YnY_{n} называется ограниченной по вероятности, если для любого ε>0\varepsilon >0 существуют MM и n0n_{0}, такие что P(∣Yn∣>M)<εP\left(\left|Y_{n}\right|>M\right)<\varepsilon для всех n>n0n>n_{0}. Покажите, что если YnY_{n} сходится по распределению, то YnY_{n} ограничена по вероятности.

?
Задача 1.8.25

В обобщение обозначений oo и OO, будем говорить, что Yn=op(1/kn)Y_{n}=o_{p}\left(1 / k_{n}\right), если knYn→0k_{n} Y_{n} \rightarrow 0 по вероятности, и что Yn=Op(1/kn)Y_{n}=O_{p}\left(1 / k_{n}\right), если knYnk_{n} Y_{n} ограничена по вероятности. Покажите, что результаты Задач 8.21-8.23 остаются справедливыми, если oo и OO заменить на opo_{p} и OpO_{p}.

?
Задача 1.8.26

Пусть (Xn,Yn)(X_{n}, Y_{n}) имеет двумерное нормальное распределение со средними E(Xn)=E(Yn)=0E\left(X_{n}\right)=E\left(Y_{n}\right)=0, дисперсиями E(Xn2)=E(Yn2)=1E\left(X_{n}^{2}\right)=E\left(Y_{n}^{2}\right)=1 и коэффициентом корреляции ρn\rho_{n}, стремящимся к 1 при n→∞n \rightarrow \infty.

?
(a)

Покажите, что (Xn,Yn)→L(X,Y)\left(X_{n}, Y_{n}\right) \xrightarrow {\mathcal{L}}(X, Y), где XX есть N(0,1)N(0,1) и P(X=Y)=1P(X=Y)=1.

(b)

Если S={(x,y):x=y}S=\left\{ (x, y): x=y\right\}, покажите, что (8.25) не выполняется.

Задача 1.8.27

Теорема 8.22 утверждает: предположим, что [n(Y1n−θ1),…,n(Yrn−θr)][\sqrt{n}(Y_{1n}-\theta_{1}), \ldots , \sqrt{n}(Y_{rn}-\theta_{r})] сходится по распределению к многомерному нормальному распределению с вектором средних 00 и ковариационной матрицей Σ\Sigma, и предположим, что h1,…,hrh_{1}, \ldots , h_{r} — rr вещественнозначных функций от θ=(θ1,…,θr)\theta =(\theta_{1}, \ldots , \theta_{r}), определённых и непрерывно дифференцируемых в окрестности ω\omega точки параметра θ\theta и таких, что матрица B=∥∂hi/∂θj∥B=\left\| \partial h_{i} / \partial \theta_{j}\right\| частных производных невырождена в ω\omega. Тогда [n[h1(Yn)−h1(θ)],…,n[hr(Yn)−hr(θ)]][\sqrt{n}[h_{1}(Y_{n})-h_{1}(\theta )], \ldots , \sqrt{n}[h_{r}(Y_{n})-h_{r}(\theta )]] сходится по распределению к многомерному нормальному распределению с вектором средних 00 и ковариационной матрицей BΣB′B\Sigma B'.

Докажите Теорему 8.22. [Указание: Проведите разложение Тейлора, как в доказательстве Теоремы 8.12, и используйте Задачу 4.16.]

?