
- •Базисные определения матричной алгебры
- •Матричные операции
- •Частные методики построения линейных моделей
- •Простая линейная регрессия
- •Множественная линейная регрессия
- •Однопутевая классификация
- •Двухпутевые классификации
- •Двухпутевая перекрестная модель
- •Число наблюдений
- •Итоговые приросты живой массы по субклассам
- •Двухпутевая гнездовая классификация
- •Модели взаимодействия
- •Модели максимального заполнения
- •Число наблюдений
- •Суммы субклассов (прироста живой массы)
- •Модели с «потерянными» данными
- •Число наблюдений
- •Итоги субклассовых наблюдений
- •Ковариантный анализ
- •Однопутевая классификация
- •Внутриклассовая регрессия
- •Смешанные модели
- •Двухпутевая перекрестная классификация
- •Способ решения смешанных моделей с заданным значением соотношения варианс
- •Представление случайных эффектов как фиксированных
- •Решение смешанных моделей при игнорировании случайных эффектов
- •Решение смешанных моделей на основе итерационного вычисления случайных эффектов
- •Число наблюдений
- •Суммарные наблюдения в субклассах (удой, кг)
- •Двухпутевая гнездовая классификация
- •Решение, основанное на априорном знании соотношения варианс
- •Решение с заменой случайных эффектов на фиксированные
- •Решение, основанное на игнорировании случайных эффектов
- •Решение, основанное на определении истинного коэффициента γ в процессе вычислений
Решение с заменой случайных эффектов на фиксированные
Решения для фиксированных эффектов (когда соотношение заранее неизвестно) представляются в виде:
.
Формализуя каждый элемент матрицы, получаем:
.
Ранг
[ХZ]
матрицы равен количеству градаций в
s-параметре
за минусом единицы, т.е.
,
поэтому для получения инверсии обычно
представляют
.
Очевидно, что при решении системы корни
уравнений не являются уникальными и
зависят от вида полученной инверсии.
Различия же между s-параметрами
не зависят от принятой структуры
инверсии и могут быть описаны при помощи
оценочной функции
с прогнозом вариансы ошибки:
Прямое
сравнение двух уровней параметра «b»
под фиксированную модель не может быть
выполнено, т.е.
не является оценочной функцией. Вместе
с тем оценочная функция представляется
как
и, очевидно, зависит от выбора М′.
Наиболее часто для оценки различий
между уровнями «b»
применяется матрица М′ следующего
представления:
.
.
Общая формула для оценки вариансы выглядит как:
,
где .
Тогда:
.
Оценка факториальной и ошибочной варианс проводится согласно:
;
.
Методика
определения
предполагает:
Вычисление
, где М′ предусматривает попарное сравнение всех s-параметров внутри каждого уровня i.
Вычисление
с
.
Тогда:
– промежуточная форма представления.
Приравняв
, получаем:
.
Математическое ожидание от МS для (S|b) определяется как
, где k1 – число степеней свободы, равное коэффициенту при .
Следовательно,
.
Решение для факториальной вариансы:
.
Подставляя
полученные численные значения
и
,
получаем оценку
.
Решение, основанное на игнорировании случайных эффектов
Фиксированные эффекты могут быть оценены с помощью операционной модели:
,
где рэндомизированные эффекты игнорируются.
В матричной форме решение представляется в виде:
.
Переходя к скалярным значениям, имеем:
.
Рант
инверсии равен (i-1)
и, положив
,
получаем решение для каждого
.
Оценить разницу между любыми
можно при помощи оценочной функции
.
Варианса разности определяется:
,
где .
В
окончательном виде
,
где K1
и K2
– числовые коэффициенты.
В данном случае величины и неизвестны, однако оценить их можно согласно следующей методики:
Определить
, где
– ранг инверсии для модели с фиксированными эффектами.
Используя модель с представлением случайных эффектов как фиксированных и наложив методику расчета, описанную в п. 2.7.2.2, получаем:
,
t1 и t2 – определенные коэффициенты.
В результате, такой способ представления эффектов дает очень приближенное значение прогноза отдельных параметров и их влияния на результирующую наблюдаемую переменную. Поэтому его применение в практике очень ограничено.