Модели со случайными и смешанными эффектами
[17/0%]В Примере 5.1:
Покажите, что совместная плотность задаётся формулой (5.2).
Получите совместную многомерную нормальную плотность напрямую, вычислив их ковариационную матрицу, а затем обратив её.
Пусть — невырожденная матрица размера с и для всех . Определите элементы .
Проверьте UMVU-оценку , приведённую в Примере 5.1.
Получите совместную плотность в Примере 5.1 в несбалансированном случае, когда , причём не все равны, и определите минимальный набор достаточных статистик (который зависит от числа различных значений ).
В сбалансированной однофакторной схеме из Примера 5.1 определите при для и .
В предыдущей задаче вычислите значения для конечного . Когда следует ожидать, что отрицательные оценки станут проблемой?
Вероятность , включающая случайную величину , может быть также выражена через неполную бета-функцию, значения которой легко доступны как в подробных таблицах, так и в компьютерных пакетах. Searle и др. (1992, Раздел 3.5d) рассматривают эту задачу подробно.
Следующая задача показывает, что в Примерах 5.1–5.3 каждая несмещённая оценка компонент дисперсии (кроме ) принимает отрицательные значения. (Некоторые связанные результаты см. в Pukelsheim 1981.) Пусть имеет распределение , и предположим, что — полная достаточная статистика для . Если — произвольная -оцениваемая функция, определённая на , и её UMVU-оценка принимает отрицательные значения с вероятностью , то покажите, что это верно для каждой несмещённой оценки .
Модифицируйте иллюстрацию с автомобилями из Примера 5.1 так, чтобы она иллюстрировала (5.5).
В Примере 5.2 определите линейное преобразование , приводящее к совместному распределению , указанному в связи с (5.6), и проверьте полные достаточные статистики (5.7).
В Примере 5.2 получите UMVU-оценки компонент дисперсии и , когда , и сравните их с оценками, полученными без этого предположения.
Для , заданных в (5.8), определите преобразование, переводящее их в переменные с распределением, указанным в Примере 5.3.
В Примере 5.3 получите UMVU-оценки компонент дисперсии и .
В Примере 5.3 получите UMVU-оценки и , когда , так что члены в (5.8) исчезают, и сравните их с оценками из Задачи 5.12.
В Примере 5.4:
Приведите преобразование, переводящее переменные в с плотностью (5.11).
Получите UMVU-оценки и .
Общий класс моделей, содержащий линейные модели Типа I и II, а также смешанные модели как частные случаи, предполагает, что вектор наблюдений размера нормально распределён со средним , как в (4.13), и с ковариационной матрицей , где — компоненты дисперсии, а — известные симметричные положительно полуопределённые матрицы размера . Покажите, что следующие модели относятся к этому типу, и в каждом случае укажите и :
(5.1 );
(5.5);
(5.5) без членов ;
(5.8);
(5.10).
Рассмотрим вложенную трёхфакторную схему с
в вариантах
;
;
;
;
где , и — неизвестные константы, однозначно определяемые обычными соглашениями, а , , и — ненаблюдаемые случайные величины, независимо и нормально распределённые со средними, равными нулю, и с дисперсиями и . В каждом случае преобразуйте в независимые переменные и получите UMVU-оценки неизвестных параметров.
Для ситуации Примера 5.5 ослабьте предположение нормальности, предположив лишь, что и имеют нулевые средние и конечные вторые моменты. Покажите, что среди всех линейных оценок (вида известны) UMVU-оценка (наилучший линейный предиктор) задаётся формулой (5.14).
Это теорема Гаусса—Маркова для предсказания в смешанных моделях. См. Harville (1976) для обобщений.