3.5

Модели со случайными и смешанными эффектами

[17/0%]
Показать
LaTeX
Задача 3.5.1

В Примере 5.1:

?
(a)

Покажите, что совместная плотность ZijZ_{i j} задаётся формулой (5.2).

(b)

Получите совместную многомерную нормальную плотность XijX_{i j} напрямую, вычислив их ковариационную матрицу, а затем обратив её.

Задача 3.5.2

Пусть A=(aij)A=\left(a_{i j}\right) — невырожденная матрица размера n×nn \times n с aii=aa_{i i}=a и aij=ba_{i j}=b для всех i≠ji \neq j. Определите элементы A−1A^{-1}.

?
Задача 3.5.3

Проверьте UMVU-оценку σA2/σ2\sigma_{A}^{2} / \sigma^{2}, приведённую в Примере 5.1.

?
Задача 3.5.4

Получите совместную плотность XijX_{i j} в Примере 5.1 в несбалансированном случае, когда j=1,…,nij=1, \ldots , n_{i}, причём nin_{i} не все равны, и определите минимальный набор достаточных статистик (который зависит от числа различных значений nin_{i}).

?
Задача 3.5.5

В сбалансированной однофакторной схеме из Примера 5.1 определите lim⁡P(σ^A2<0)\lim \mathbb {P}\left(\hat{\sigma }_{A}^{2}<0\right) при n→∞n \rightarrow \infty для σA2/σ2=0,0.2,0.5,1\sigma_{A}^{2} / \sigma^{2}=0,0.2,0.5,1 и s=3,4,5,6s=3,4,5,6.

?
Задача 3.5.6

В предыдущей задаче вычислите значения P(σ^A2<0)\mathbb {P}\left(\hat{\sigma }_{A}^{2}<0\right) для конечного nn. Когда следует ожидать, что отрицательные оценки станут проблемой?

?
Примечание.
?

Вероятность P(σ^A2<0)\mathbb {P}\left(\hat{\sigma }_{A}^{2}<0\right), включающая случайную величину FF, может быть также выражена через неполную бета-функцию, значения которой легко доступны как в подробных таблицах, так и в компьютерных пакетах. Searle и др. (1992, Раздел 3.5d) рассматривают эту задачу подробно.

Задача 3.5.7

Следующая задача показывает, что в Примерах 5.1–5.3 каждая несмещённая оценка компонент дисперсии (кроме σ2\sigma^{2}) принимает отрицательные значения. (Некоторые связанные результаты см. в Pukelsheim 1981.) Пусть XX имеет распределение P∈PP \in \mathcal{P}, и предположим, что TT — полная достаточная статистика для P\mathcal{P}. Если g(P)g(P) — произвольная UU-оцениваемая функция, определённая на P\mathcal{P}, и её UMVU-оценка η(T)\eta (T) принимает отрицательные значения с вероятностью >0>0, то покажите, что это верно для каждой несмещённой оценки g(P)g(P).

?
Задача 3.5.8

Модифицируйте иллюстрацию с автомобилями из Примера 5.1 так, чтобы она иллюстрировала (5.5).

?
Задача 3.5.9

В Примере 5.2 определите линейное преобразование XijkX_{i j k}, приводящее к совместному распределению ZijkZ_{i j k}, указанному в связи с (5.6), и проверьте полные достаточные статистики (5.7).

?
Задача 3.5.10

В Примере 5.2 получите UMVU-оценки компонент дисперсии σA2,σB2\sigma_{A}^{2}, \sigma_{B}^{2} и σ2\sigma^{2}, когда σC2=0\sigma_{C}^{2}=0, и сравните их с оценками, полученными без этого предположения.

?
Задача 3.5.11

Для XijkX_{i j k}, заданных в (5.8), определите преобразование, переводящее их в переменные ZijkZ_{i j k} с распределением, указанным в Примере 5.3.

?
Задача 3.5.12

В Примере 5.3 получите UMVU-оценки компонент дисперсии σA2,σB2\sigma_{A}^{2}, \sigma_{B}^{2} и σ2\sigma^{2}.

?
Задача 3.5.13

В Примере 5.3 получите UMVU-оценки σA2\sigma_{A}^{2} и σ2\sigma^{2}, когда σB2=0\sigma_{B}^{2}=0, так что члены BB в (5.8) исчезают, и сравните их с оценками из Задачи 5.12.

?
Задача 3.5.14

В Примере 5.4:

?
(a)

Приведите преобразование, переводящее переменные XijkX_{i j k} в WijkW_{i j k} с плотностью (5.11).

(b)

Получите UMVU-оценки μ,αi,σB2\mu , \alpha_{i}, \sigma_{B}^{2} и σ2\sigma^{2}.

Задача 3.5.15

Общий класс моделей, содержащий линейные модели Типа I и II, а также смешанные модели как частные случаи, предполагает, что вектор наблюдений X\mathbf{X} размера 1×n1 \times n нормально распределён со средним θA\boldsymbol {\theta } A, как в (4.13), и с ковариационной матрицей Σi=1mγiVi\Sigma_{i=1}^{m} \gamma_{i} V_{i}, где γ\gamma — компоненты дисперсии, а ViV_{i} — известные симметричные положительно полуопределённые матрицы размера n×nn \times n. Покажите, что следующие модели относятся к этому типу, и в каждом случае укажите γ\gamma и VV:

?
(a)

(5.1 );

(b)

(5.5);

(c)

(5.5) без членов CijC_{i j};

(d)

(5.8);

(e)

(5.10).

Задача 3.5.16

Рассмотрим вложенную трёхфакторную схему с

Xijkl=μ+αi+bij+cijk+Uijkl X_{i j k l}=\mu +\alpha _{i}+b_{i j}+c_{i j k}+U_{i j k l}

(i=1,…,I;j=1,…,J;k=1,…,K;l=1,…,n)(i=1, \ldots , I ; j=1, \ldots , J ; k=1, \ldots , K ; l=1, \ldots , n) в вариантах

?
(a)

ai=αi,bij=βij,cijk=γijka_{i}=\alpha_{i}, b_{i j}=\beta_{i j}, c_{i j k}=\gamma_{i j k};

(b)

ai=αi,bij=βij,cijk=Cijka_{i}=\alpha_{i}, b_{i j}=\beta_{i j}, c_{i j k}=C_{i j k};

(c)

ai=αi,bij=Bij,cijk=Cijka_{i}=\alpha_{i}, b_{i j}=B_{i j}, c_{i j k}=C_{i j k};

(d)

ai=Ai,bij=Bij,cijk=Cijka_{i}=A_{i}, b_{i j}=B_{i j}, c_{i j k}=C_{i j k};

где α\alpha, β\beta и γ\gamma — неизвестные константы, однозначно определяемые обычными соглашениями, а AA, BB, CC и UU — ненаблюдаемые случайные величины, независимо и нормально распределённые со средними, равными нулю, и с дисперсиями σA2,σB2,σC2\sigma_{A}^{2}, \sigma_{B}^{2}, \sigma_{C}^{2} и σ2\sigma^{2}. В каждом случае преобразуйте XijklX_{i j k l} в независимые переменные ZijklZ_{i j k l} и получите UMVU-оценки неизвестных параметров.

Задача 3.5.17

Для ситуации Примера 5.5 ослабьте предположение нормальности, предположив лишь, что AiA_{i} и UijU_{i j} имеют нулевые средние и конечные вторые моменты. Покажите, что среди всех линейных оценок (вида ∑cijxij,cij\sum c_{i j} x_{i j}, c_{i j} известны) UMVU-оценка μ+αi\mu +\alpha_{i} (наилучший линейный предиктор) задаётся формулой (5.14).

?
Примечание.
?

Это теорема Гаусса—Маркова для предсказания в смешанных моделях. См. Harville (1976) для обобщений.