16.2. ВЫБОРКИ ИЗ МНОГОМЕРНЫХ НОРМАЛЬНЫХ (MVN) РАСПРЕДЕЛЕНИЙ
Оценки некоторых параметров совокупности в многомерном случае можно получить по выборочным данным таким же образом, как и в одномерном, что было показано в разделе 16.1. Однако не все методы многомерного анализа представляют собой простые аналоги
методов, применяемых в одномерном случае. Большинство методов, рассматриваемых в данной главе, не имеет эквивалента в одномерном анализе.
16.2.1. ОЦЕНКИ МАКСИМАЛЬНОГО ПРАВДОПОДОБИЯ
Для заданной выборки объема к из многомерного нормального распределения [см. II, раздел 13.4] со средним вектором и матрицей ковариаций V (для обозначения этого распределения возьмем аббревиатуру функция правдоподобия как функция параметров и V имеет вид
так как векторные наблюдения независимы. (Заметим, что компоненты векторной случайной переменной не будут взаимонезависимыми, если V не является единичной матрицей I.) Лог-правдоподобие есть
Последний член, который зависит только от может быть записан в форме
откуда видно, что лог-правдоподобие максимизируется при значении следовательно, х — оценка максимального правдоподобия для
Проблема оценивания матрицы V не так проста. Поскольку матрица V — симметрична [см. I, раздел 6.7], имеется не элементов, которые нужно оценить, среди них дисперсий и ковариаций.
Уравнения правдоподобия удобнее переписать, взяв матрицу вместо V. Тогда определитель заменяется на Оценивая элементы в и используя многомерный аналог свойства инвариантности оценок максимального правдоподобия [см. раздел 6.2.6], можно показать, что оценкой максимального правдоподобия для V будет
где матрица А определена выражением (16.1.9). Эту оценку можно сделать несмещенной в результате умножения на Детали доказательства приведены в работах [Andersen (1958), с. 47—48] и [Morrison (1976), с. 99—100]. Рассмотренные результаты аналогичны результатам, полученным для оценок максимального правдоподобия в случае одномерного нормального распределения [см. пример 6.4.1].
Пример 16.2.1. Оценка дисперсионной матрицы двумерного нормального распределения. Рассмотрим случайную выборку
из двумерного нормального распределения с вектором средних
и матрицей ковариаций
Обратная матрица для V будет
Членом лог-правдоподсбяя, соответствующим последнему члену в (16.2.2), будет
Дифференцируем его сначала по а затем по Получаем, что максимум функции лог-правдоподобия достигается, когда
и
т. е. когда
и
Следовательно, оценка максимального правдоподобия (ОМП) для есть
а для
так что
Члены лог-правдоподобия, включающие можно представить в виде
Дифференцируя это выражение получим следующие уравнения правдоподобия ОМП для
Из свойства инвариантности ОМП вытекает, что ОМП для коэффициента корреляции
есть
и знак совпадает со знаком
Поскольку
из записанных выше уравнений следует, что
и
где — выборочные дисперсии, а — выборочная ковариация [см. (16.1.9)].