Асимптотическая оптимальность
[169/0%]Пусть — независимые одинаково распределённые случайные величины с .
Если имеют конечный четвёртый момент, установите (1.3)
Для положительного целого покажите, что и , если они существуют, оба являются .
Для фиксированного опишите относительную ошибку в Примере 1.3 как функцию от .
Теорема 1.1 утверждает: пусть — независимые одинаково распределённые случайные величины с , и конечным четвёртым моментом, и пусть — функция вещественной переменной, первые четыре производные которой существуют для всех , где — интервал с . Предположим далее, что для всех при некотором . Тогда,
и если, кроме того, четвёртая производная от также ограничена,
где остаточный член в обоих случаях есть .
Теорема 1.5 утверждает: предположим, что предположения Теоремы 1.1 выполнены, и что — последовательность констант, удовлетворяющая . Тогда дисперсия удовлетворяет
Докажите Теорему 1.5.
Пусть — независимые одинаково распределённые случайные величины с распределением , известно, и пусть . Определите, с точностью до членов порядка ,
дисперсию UMVU-оценки ;
смещение MLE для .
Пусть — независимые одинаково распределённые случайные величины с распределением , известно. Для чётного определите дисперсию UMVU-оценки (2.2.4) величины с точностью до членов порядка .
Решите предыдущую задачу для случая, когда неизвестно.
Для оценивания в Примере 3.3.1 определите, с точностью до порядка ,
дисперсию UMVU-оценки (2.3.2);
смещение MLE.
Решите предыдущую задачу, если заменить оцениваемым параметром из Задачи 2.3.3.
Пусть — независимые одинаково распределённые случайные величины с распределением Пуассона .
Определите UMVU-оценку .
Вычислите дисперсию оценки из пункта (a) с точностью до членов порядка .
Решите пункт (b) предыдущей задачи для оценки (2.3.22).
В условиях Задачи 1.1 покажите, что . [Указание: Используйте тот факт, что , вместе с результатом Задачи 1.1.]
Получите вариант Теоремы 1.1, который требует существования и ограниченности только вместо , но где есть лишь .
Чтобы убедиться, что Теорема 1.1 не обязательно верна без ограниченности четвёртой (или некоторой более высокой) производной, предположим, что распределены как , и пусть . Тогда все моменты и все производные существуют.
Покажите, что математическое ожидание не существует ни при каком , и, следовательно, что при всех значениях .
С другой стороны, покажите, что имеет предельное нормальное распределение с конечной дисперсией, и определите эту дисперсию.
Пусть — независимые одинаково распределённые случайные величины из показательного распределения с плотностью , и .
Используя Теорему 1.1, найдите приближения для и .
Проверьте точное вычисление
и покажите, что .
Согласуйте результаты пунктов (a) и (b). Объясните, почему, хотя Теорема 1.1 неприменима, она дала верный ответ.
Покажите, что аналогичный вывод верен для .
Пусть — независимые одинаково распределённые случайные величины с распределением . Определите дисперсию UMVU-оценки , где — целое число, .
В условиях Задачи 1.15 найдите MLE для и сравните её ожидаемую квадратичную ошибку с дисперсией UMVU-оценки.
Пусть — независимые одинаково распределённые случайные величины с распределением , пусть , и пусть — функция, удовлетворяющая условиям Теоремы 1.1. Покажите, что
и
Примените результаты Задачи 1.17, чтобы получить приближённые ответы к Задачам 1.15 и 1.16, и сравните ответы с точными решениями.
Если таковы, как в Теореме 1.1, и если первые пять производных существуют, а пятая производная ограничена, покажите, что
и, если пятая производная также ограничена,
где .
В условиях предыдущей задачи проведите вычисление дисперсии (1.16) до членов порядка и сравните результат с результатом предыдущей задачи.
Проведите вычисление из Задачи 1.4 до членов порядка .
Для оцениваемых параметров из Задачи 1.4 вычислите ожидаемую квадратичную ошибку MLE до членов порядка и сравните её с дисперсией, вычисленной в Задаче 1.21.
Вычислите дисперсию (1.18) до членов порядка и сравните её с ожидаемой квадратичной ошибкой MLE, вычисленной до того же порядка.
Найдите дисперсию оценки (2.3.17) с точностью до членов порядка .
Для ситуации Примера 1.12 покажите, что UMVU-оценка является MLE с поправкой на смещение, где MLE есть .
Для оценок из Примера 1.13:
Вычислите их точные дисперсии.
Используя результат пункта (a), проверьте (1.27).
В условиях Теоремы 1.5, если все четвёртые моменты конечны, покажите, что и что все третьи и четвёртые моменты и так далее имеют порядок .
Если, кроме того, все производные суммарного порядка существуют, а производные порядка 4 равномерно ограничены, то
и если производные порядка 4 также ограничены,
где остаточный член в обоих случаях есть .
На основе выборки из пусть — вероятность того, что UMVU-оценка величины ( известно) отрицательна.
Покажите, что — убывающая функция от .
Покажите, что при для любых фиксированных и .
Определите значение .
Используя -распределение, найдите значение из предыдущей задачи для UMVU-оценки , когда неизвестно, для характерных значений .
Теорема 1.9 утверждает: если последовательность случайных величин стремится к по вероятности, и если непрерывна в точке , то . Чтобы показать, что , достаточно заметить, что в силу непрерывности разность будет сколь угодно мала, если близко к , и что для любого вероятность того, что , стремится к при .
Заполните детали доказательства Теоремы 1.9. (См. также Задачу 1.8.8.)
В Примере 8.13 при покажите, что не имеет в точности распределение .
В Примере 8.13 пусть , что является улучшением по сравнению с .
Покажите, что имеет то же предельное распределение, что и , при .
Опишите предельное распределение при .
Пусть имеет биномиальное распределение , и пусть . UMVU-оценкой является . Определите предельное распределение и , когда и соответственно.
Пусть — независимые одинаково распределённые случайные величины с распределением . Определите предельное поведение распределения UMVU-оценки .
Определите предельное поведение оценки (2.3.22) при .
Пусть — независимые одинаково распределённые случайные величины с распределением , и предположим, что является UMVU-оценкой на основе . Если существует и несмещённая оценка , имеющая конечную дисперсию при всех , то состоятельна для .
Пусть имеет распределение с вероятностью и распределение с вероятностью . Если и , определите, при каких значениях последовательность имеет предельное распределение, а при каких — нет.
В Задаче 1.37 определите, к каким значениям может стремиться при , если и , но в остальном обе последовательности произвольны.
Используя (a), покажите, что предел дисперсии не обязан совпадать с дисперсией предельного распределения.
Пусть — двойная последовательность вещественных чисел, которая при каждом фиксированном не убывает по . Покажите, что и , при условии что указанные пределы существуют (они могут быть бесконечны), и где и обозначают, соответственно, наименьшую и наибольшую предельные точки, достижимые последовательностью , с и .
Пусть — независимые одинаково распределённые случайные величины с распределением . Рассмотрим две оценки
где , и с — наибольшим целым числом .
Покажите, что асимптотическая эффективность относительно равна нулю.
Покажите, что для любого фиксированного , но .
Для больших значений что можно сказать о двух вероятностях из пункта (b), когда заменяется на ? (Basu 1956).
Если для некоторой последовательности , покажите, что тот же результат верен, если заменить на , где .
Предположим, что распределение сходится к распределению со средним 0 и дисперсией . Используя лемму Фату (Лемма 1.2.6), установите, что для всех .
Если — выборка из однопараметрического экспоненциального семейства (1.5.2), то является минимальной достаточной статистикой, и . Покажите, что для любой функции , для которой верна Теорема 1.8.12, асимптотически несмещена для .
Если — независимые одинаково распределённые случайные величины с распределением , покажите, что является асимптотически несмещённой оценкой .
Пусть — независимые одинаково распределённые случайные величины с распределением . Согласно Примеру 2.1.14, является UMVU-оценкой , тогда как MLE есть . Определите предельное распределение
и
.
Прокомментируйте асимптотическое смещение этих оценок.
Для ситуации Задачи 2.6:
Вычислите средние квадратичные ошибки и как оценок .
Покажите, что
Проверьте асимптотическое распределение, заявленное для в Примере 2.5.
Пусть — произвольная оценка, удовлетворяющая (2.2) с . Постройте последовательность , такую что с при и .
В предыдущей задаче постройте , такую что для всех и , и при и .
Постройте последовательность , удовлетворяющую (2.2), но для которой смещение не стремится к нулю.
В Примере 2.7 с , заданной формулой (2.11), покажите, что при и что .
Пусть — смещение оценки из Примера 2.5.
Покажите, что
Покажите, что при любом и .
Используя (b), объясните, как оценка Ходжеса может нарушать (2.7), не нарушая при этом информационное неравенство.
В Примере 2.7 покажите, что если , то .
Пусть имеет биномиальное распределение . Определите MLE для
обычным методом математического анализа, определяя максимум функции;
показав, что .
В предыдущей задаче покажите, что MLE не существует, когда ограничено условием , а или .
Пусть — независимые одинаково распределённые случайные величины с распределением . Определите MLE для
, когда известно,
, когда известно, и
, когда оба неизвестны.
Предположим, что — независимые одинаково распределённые случайные величины с распределением , где . Покажите, что MLE есть , когда , и не существует, когда .
Пусть принимает значения 0 и 1 с вероятностями и соответственно. Если известно, что ,
найдите MLE и
покажите, что ожидаемая квадратичная ошибка MLE равномерно больше, чем у .
[Аналогичная задача оценивания возникает в опросах с рандомизированным ответом. См. Пример 5.2.2.]
Когда конечно, покажите, что MLE состоятельна тогда и только тогда, когда она удовлетворяет (3.2).
Покажите, что Теорема 3.2 остаётся верной, если предположение A1 ослабить до A1': Существует непустое множество , такое что , и содержится в носителе каждого .
Докажите существование и единственность , удовлетворяющих (3.4).
Докажите (3.9).
В Примере 3.6 с определите состоятельную оценку .
Проверьте характер корней в Примере 3.9.
Пусть имеет распределение . Покажите, что при условии переменная стремится по вероятности к при .
Рассмотрим выборку из распределения Пуассона, обусловленного положительностью, так что при . Покажите, что уравнение правдоподобия имеет единственный корень при всех значениях .
Пусть имеет отрицательное биномиальное распределение (2.3.3). Найдите ELE для .
Плотность называется строго унимодальной, или, что то же самое, логарифмически вогнутой, если является вогнутой функцией. Покажите, что такая плотность имеет единственную моду.
Пусть — независимые одинаково распределённые случайные величины с плотностью . Покажите, что функция правдоподобия имеет единственный корень, если монотонна, и что этот корень является максимумом, если убывает. Следовательно, логарифмически вогнутые плотности дают единственные MLE.
Пусть — положительные случайные величины (или симметрично распределённые относительно нуля) с совместной плотностью . Покажите, что уравнение правдоподобия имеет единственный максимум, если строго убывает при .
Если — независимые одинаково распределённые случайные величины с плотностью , где унимодальна, и если уравнение правдоподобия имеет единственный корень, покажите, что уравнение правдоподобия также имеет единственный корень, когда плотность каждой равна , при известном .
Для каждой из следующих плотностей определите,
(a) является ли она строго унимодальной, и
(b) является ли строго убывающей при .
Соответственно прокомментируйте, имеют ли соответствующие параметры сдвига и масштаба единственные MLE:
Если — независимые одинаково распределённые случайные величины с плотностью или , и — логистическая плотность , уравнение правдоподобия имеет единственные решения и как в случае сдвига, так и в случае масштаба. Определите предельное распределение и .
В Задаче 3.15(b), при — плотности Коши , уравнение правдоподобия имеет единственный корень , и .
Если — независимые одинаково распределённые случайные величины с распределением , то при любом фиксированном существует положительная вероятность того,
что уравнение правдоподобия имеет корней, и
что уравнение правдоподобия имеет единственный корень.
Если — независимые одинаково распределённые случайные величины с гамма-распределением , уравнение правдоподобия имеет единственный корень.
Пусть — независимые одинаково распределённые случайные величины с распределением Вейбулла с плотностью
которое не принадлежит ни экспоненциальному семейству, ни семейству сдвига, ни семейству масштаба. Тем не менее покажите, что существует единственный внутренний максимум функции правдоподобия.
В условиях Теоремы 3.2 покажите, что
стремится по распределению к .
Пусть — независимые одинаково распределённые случайные величины с распределением , где — известная положительная константа.
Найдите явное выражение для ELE параметра .
Определите, существует ли MRE-оценка при подходящей группе преобразований.
Этот случай рассматривался в Berk (1972).
Проверьте, что предположения Теоремы 3.10 выполнены в Примере 3.12.
Для — независимых одинаково распределённых случайных величин с распределением покажите, что
выборочная медиана является MLE для , и
выборочная медиана асимптотически нормальна с дисперсией — границей информационного неравенства.
В Примере 3.12 покажите непосредственно, что является асимптотически эффективной оценкой , рассмотрев её предельное распределение.
Пусть — независимые одинаково распределённые случайные величины с распределением , где — пара заданных плотностей вероятности относительно , и .
Приведите один пример , для которого предположения Теоремы 3.10 выполнены, и один, для которого они не выполнены.
Обсудите существование и характер корней уравнения правдоподобия при , 2, 3.
В условиях Теоремы 3.7 предположим, что и — две состоятельные последовательности корней уравнения правдоподобия. Докажите, что при . [Указание:
Пусть . Для всех существует между и , такое что . Для всех пусть — общее значение и . Тогда — состоятельная последовательность корней уравнения правдоподобия.
по вероятности, и, следовательно, по вероятности.
Пусть , и пусть
Тогда . С другой стороны, на , так что содержится в дополнении к (Huzurbazar 1948).]
Чтобы установить измеримость последовательности корней из Теоремы 3.7, можно следовать доказательству Serfling (1980, Раздел 4.2.2), где доказана измеримость аналогичной последовательности.
Для определённости определим как значение, минимизирующее при условии
Покажите, что измеримо.
Покажите, что — корень, ближайший к , — измерим.
Пусть
Покажите, что при подходящем функция непрерывна и бесконечно дифференцируема при .
Покажите, что плотность (4.1) с удовлетворяет всем условиям Теоремы 3.10, за исключением условия (d) Теоремы 2.6.
Покажите, что плотность (4.4) с удовлетворяет всем условиям Теоремы 3.10.
В Примере 4.5 вычислите оценки (4.8) и (4.14) для случая Коши, используя в качестве выборочную медиану.
В Примере 4.7 покажите, что вогнута.
В Примере 4.7, если , покажите, как получить -состоятельную оценку, приравнивая выборочный и теоретический вторые моменты.
В Теореме 4.8 покажите, что .
Не используя Теорему 4.8, покажите в Примере 4.13, что EM-последовательность сходится к MLE.
Рассмотрим следующие 12 наблюдений из двумерного нормального распределения с параметрами :
| 1 | 1 | -1 | -1 | 2 | 2 | -2 | -2 | * | * | * | * | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 1 | -1 | 1 | -1 | * | * | * | * | 2 | 2 | -2 | -2 |
где «» обозначает пропущенное значение.
Покажите, что функция правдоподобия имеет глобальные максимумы при , и седловую точку при .
Покажите, что если EM-последовательность начинается с , то она остаётся при на всех последующих итерациях.
Покажите, что если EM-последовательность начинается с , отделённого от нуля, то она сходится к максимуму.
Эта задача принадлежит Murray (1977) и обсуждается в Wu (1983).
Покажите, что если плотность полных данных EM из (4.21) принадлежит искривлённому экспоненциальному семейству, то предположения Теоремы 4.12 выполнены.
В EM-алгоритме вычисление E-шага, то есть вычисление математического ожидания, может быть сложным. В таких случаях E-шаг иногда можно заменить методом Монте-Карло, что приводит к MCEM-алгоритму (Wei and Tanner 1990). Рассмотрим следующее MCEM-вычисление : При заданном
-
Сгенерируйте — независимые одинаково распределённые случайные величины из ,
-
Положите
а затем вычислите как значение, максимизирующее .
Покажите, что при .
Какие условия обеспечат, что при достаточно больших ? Достаточны ли для этого предположения Теоремы 4.12?
Для смеси распределений из Примера 4.7, то есть
где и известны, для нахождения ML-оценки можно использовать EM-алгоритм. Пусть , где указывает, из какого распределения получена , так что
Покажите, что функцию правдоподобия полных данных можно записать в виде
Покажите, что , и, следовательно, EM-последовательность задаётся формулой
Покажите, что — ML-оценка .
Для ситуации Примера 4.10:
Покажите, что M-шаг EM-алгоритма задаётся формулами
Покажите, что E-шаг EM-алгоритма задаётся формулой
При ограничении покажите, что EM-последовательность сходится к , где .
При ограничении покажите, что EM-последовательность сходится к , где .
Для общей однофакторной схемы с обработками и наблюдениями на обработку покажите, как использовать EM-алгоритм для дополнения данных так, чтобы каждая обработка имела наблюдений. Запишите EM-последовательность и покажите, к чему она сходится при ограничениях пунктов (c) и (d).
[Ограничения пунктов (c) и (d) встречались в Примере 3.4.9, где они привели, соответственно, к анализу невзвешенных средних и анализу взвешенных средних.]
В двухфакторной схеме (см. Пример 3.4.11) EM-алгоритм может быть очень полезен при вычислении ML-оценок в несбалансированном случае. Предположим, что мы наблюдаем
где . Данные будут дополнены так, чтобы полные данные содержали наблюдений на ячейку.
Покажите, как вычислить и E-шаг, и M-шаг EM-алгоритма.
При ограничении покажите, что EM-последовательность сходится к ML-оценкам, соответствующим анализу невзвешенных средних.
При ограничении покажите, что EM-последовательность сходится к ML-оценкам, соответствующим анализу взвешенных средних.
Для однофакторной схемы со случайными эффектами (Пример 3.5.1) EM-алгоритм полезен для вычисления ML-оценок. (На самом деле, он весьма полезен во многих смешанных моделях; см. Searle и др. 1992, Глава 8.) Предположим, что имеет место модель
где и — независимые нормальные случайные величины с нулевым средним и известной дисперсией. Для вычисления ML-оценок и обычно применяют EM-алгоритм, используя ненаблюдаемые в качестве дополненных данных. Выпишите E-шаг и M-шаг и покажите, что EM-последовательность сходится к ML-оценкам.
Оценивание максимального правдоподобия в пробит-модели из Раздела 3.6 можно реализовать с помощью EM-алгоритма. Мы наблюдаем независимые бернуллиевские величины , которые зависят от ненаблюдаемых величин , независимо распределённых как , где
Предполагая, что известно, мы хотим получить ML-оценки и .
Покажите, что функция правдоподобия равна , где
Если считать полными данными, то правдоподобие полных данных равно
а ожидаемый логарифм правдоподобия полных данных равен
Покажите, что EM-последовательность задаётся формулами
где
Покажите, что
где
Покажите, что и — ML-оценки и .
Проверьте (4.30).
EM-алгоритм также можно применить в байесовской иерархической модели для нахождения моды апостериорного распределения. Вспомним модель (4.5.5.1),
где интерес представляет оценивание величин, связанных с . Поскольку
где , EM-алгоритм является одним из возможных методов нахождения моды , где используется в качестве дополненных данных.
Определите и покажите, что
Если последовательность удовлетворяет
покажите, что . При каких условиях последовательность будет сходиться к моде ?
Для иерархии
покажите, как использовать EM-алгоритм для вычисления моды апостериорного распределения .
Между EM-алгоритмом и сэмплированием Гиббса существует связь, поскольку оба метода основаны на теории марковских цепей. Один из способов убедиться в этом — показать, что правдоподобие неполных данных является решением интегрального уравнения выборки методом последовательных подстановок (см. Задачи 4.5.9–4.5.11), и что затем сэмплирование Гиббса можно использовать для вычисления функции правдоподобия. Если — правдоподобие неполных данных, а — правдоподобие полных данных, определим
Покажите, что является решением уравнения
где, как обычно, .
Покажите, как последовательность , полученная из итерации Гиббса,
будет сходиться к случайной величине с плотностью при . Как это можно использовать для вычисления функции правдоподобия ?
[Используя функции и , EM-алгоритм даёт нам ML-оценку из , тогда как сэмплер Гиббса даёт нам всю функцию целиком. Эта реализация сэмплера Гиббса для вычисления правдоподобия использовалась Casella и Berger (1994), а также описана у Smith и Roberts (1993). Версия EM-алгоритма, в которой связь с марковскими цепями вполне очевидна, была предложена Baum и Petrie (1966) и Baum и др. (1970).]
Если вектор в сходится по вероятности к постоянному вектору , а — непрерывная функция, определённая на , покажите, что по вероятности.
Используя (a), покажите, что элементы сходятся по вероятности к элементам , как утверждается в доказательстве Леммы 5.2.
Покажите, что (5.26) без учёта остаточного члена имеет тот же вид, что и (5.15), и определите .
Покажите, что получаемые из Леммы 5.2 совпадают с теми, что в (5.23).
Покажите, что остаточным членом в (5.26) можно пренебречь в доказательстве Теоремы 5.3.
Пусть — независимые одинаково распределённые случайные величины с распределением .
Покажите, что уравнения правдоподобия имеют единственный корень.
Покажите непосредственно (т.е. не прибегая к Теореме 5.1), что MLE и асимптотически эффективны.
Пусть имеет полиномиальное распределение .
Покажите, что уравнения правдоподобия имеют единственный корень.
Покажите непосредственно, что MLE асимптотически эффективны.
Теорема 5.3 утверждает: предположим, что выполнены предположения Теоремы 5.1 (стандартные условия многопараметрической регулярности правдоподобия и состоятельности для независимой одинаково распределённой выборки из регулярного параметрического семейства с параметром ), и что — -состоятельная оценка при . Пусть — решение уравнений
где — логарифм правдоподобия, , а . Тогда асимптотически эффективна.
Следствие 5.4 утверждает: предположим, что выполнены предположения Теоремы 5.3 и что элементы информационной матрицы величин непрерывны. Тогда решения уравнений
асимптотически эффективны.
Докажите Следствие 5.4.
Покажите, что существует функция двух переменных, для которой уравнения 0 и имеют единственное решение, причём это решение является локальным, но не глобальным максимумом .
В Примере 6.1 покажите, что уравнения правдоподобия задаются формулами (6.2) и (6.3).
В Примере 6.1 проверьте уравнение (6.4).
Проверьте (6.5).
Если и если
то асимптотическое распределение при предположениях Теоремы 5.1 не зависит от того, известны ли .
Пусть — независимые одинаково распределённые случайные величины с гамма-распределением и плотностью .
Вычислите информационную матрицу для обычной параметризации .
Запишите плотность через параметры . Вычислите информационную матрицу для параметризации и покажите, что она диагональна, а значит, параметры ортогональны.
Если MLE в пункте (a) равны , покажите, что . Таким образом, обе модели одинаково хорошо оценивают среднее.
(Теорию и другие примеры ортогональности параметров см. в Cox и Reid 1987.)
В Примере 6.4 проверьте MLE и , когда неизвестны.
В Примере 6.4 покажите, что , заданные формулой (6.15), не зависят от и имеют то же совместное распределение, что и статистики (6.13), в которых заменено на .
Проверьте матрицы
(6.17) и
(6.18).
Рассмотрим ситуацию, приводящую к (6.20), где , независимы и одинаково распределены по двумерному нормальному закону с и неизвестным коэффициентом корреляции .
Покажите, что уравнение правдоподобия является кубическим, для которого вероятность единственного корня стремится к 1 при . [Указание: Для кубического уравнения положим и . Тогда условие единственности вещественного корня есть .]
Покажите, что если — состоятельное решение уравнения правдоподобия, то оно удовлетворяет (6.20).
Покажите, что — состоятельная оценка и что , а значит, менее эффективна, чем MLE для .
Проверьте предельное распределение, утверждаемое в (6.21).
Пусть — независимые одинаково распределённые по распределению Пуассона . Найдите ARE Число по отношению к как оценок .
Покажите, что эффективность (6.27) стремится к 0 при .
Для ситуации Примера 6.9 рассмотрим в качестве другого семейства распределений семейство загрязнённых нормальных смесей, предложенное Tukey (1960) в качестве модели для наблюдений, которые обычно подчиняются нормальному распределению, но иногда в эксперименте или при регистрации его результатов что-то идёт не так, так что полученное наблюдение оказывается грубой ошибкой. В модели Tukey функция распределения имеет вид
То есть в случаях грубой ошибки предполагается, что наблюдения нормально распределены с тем же средним , но с другой (большей) дисперсией
Покажите, что если имеют распределение , предельное распределение не меняется.
Покажите, что предельное распределение нормально со средним нуль и дисперсией .
Сравните асимптотическую относительную эффективность и .
Пусть независимы и одинаково распределены как .
Покажите, что является состоятельной оценкой тогда и только тогда, когда .
Определите ARE оценки при по отношению к MLE .
Пусть независимы и одинаково распределены с и , и рассмотрим несмещённые оценки и параметра .
Определите ARE оценки по отношению к .
Покажите, что , если симметричны относительно .
Найдите распределение , для которого .
Свойство асимптотической относительной эффективности было определено (Определение 6.6) для оценок, сходящихся к нормальному распределению со скоростью . Это определение, а также Теорема 6.7, можно обобщить, включив другие распределения и скорости сходимости.
Как отметил Stigler (1973), подобные модели для распределений с тяжёлыми хвостами уже предлагались значительно раньше в забытой работе Newcomb (1882, 1886). Теорема 9.1 Пусть — две последовательности оценок , такие что
где распределение величины имеет носитель на интервале со строго возрастающей функцией распределения на . Тогда ARE по отношению к существует и равна
В Примере 6.10 покажите, что условия Теоремы 5.1 выполнены.
Рассмотрим независимых выборок с индексом объёмов , таких что при стремлении общего объёма выборки имеем , где и каждое .
Теорема 7.1 утверждает: предположим, что предположения Теоремы 5.1 выполнены для каждой из плотностей . Пусть обозначает информационную матрицу, соответствующую , и пусть . Логарифм правдоподобия задаётся формулой , а уравнения правдоподобия — . При этих обозначениях выводы Теоремы 5.1 остаются в силе.
Докажите Теорему 7.1.
Для ситуации Примера 7.3 при :
Покажите, что необходимым условием сходимости (7.5) к является 0, где и , для и из (7.4).
Используя тот факт, что имеет -распределение, покажите, что .
Покажите, что полная MLE задаётся решением
и выведите её асимптотическую эффективность из Теоремы 5.1.
В Примере 7.4 определите совместное распределение
и
.
Рассмотрим выборки и из двух двумерных нормальных распределений с нулевыми средними и дисперсиями-ковариациями и соответственно. Используя Теорему 7.1 и Примеры 6.5 и 6.8, найдите предельное распределение
и , когда известно, что
, когда известно, что и .
В предыдущей задаче найдите выигрыш в эффективности (если он есть)
в пункте (a), обусловленный знанием того, что
в пункте (b), обусловленный знанием того, что и .
Покажите, что уравнения правдоподобия (7.11) имеют не более одного решения.
В Примере 7.6 предположим, что и что и строго вогнуты. Тогда уравнения правдоподобия имеют не более одного решения.
Если функция распределения симметрична и строго вогнута, то и строго вогнута.
Покажите, что строго вогнута, когда сильно унимодальна, но не когда — распределение Коши.
В Примере 7.7 покажите, что менее информативна, чем .
Покажите, что оценка из Примера 7.7 удовлетворяет (7.14).
Найдите подходящие нормирующие константы для из Примера 7.7, когда
,
, и
.
Пусть независимы, нормальны, с дисперсией 1 и средним (где известны). Обсудите оценивание по аналогии с Примером 7.7.
Обобщите предыдущую задачу на ситуацию, в которой
и , и
и , где и — неизвестные параметры, подлежащие оцениванию.
Пусть независимо распределены с плотностями ( вещественнозначный), пусть — информация о , содержащаяся в , и пусть — полная информация о в выборке. Предположим, что — состоятельный корень уравнения правдоподобия и что, обобщая (3.18)–(3.20),
и
Покажите, что
Докажите, что последовательность из Примера 7.8 стационарна при условии, что она удовлетворяет (7.17).
В Примере 7.8 покажите, что уравнение правдоподобия имеет единственное решение, что оно является MLE и что оно имеет то же асимптотическое распределение, что и .
Покажите непосредственно, что — состоятельная оценка .
В Примере 7.8:
Покажите, что при математическое ожидание условной информации (при заданном ), которую содержит о , равно .
Определите информацию о , содержащуюся в .
Покажите, что при в (7.21) MLE существует и состоятельна.
Предположим, что в (7.21) сами являются случайными величинами, независимо одинаково распределёнными как .
Покажите, что совместная плотность является плотностью выборки из двумерного нормального распределения, и определите параметры этого распределения.
В модели пункта (a) найдите асимптотически эффективные оценки параметров и .
Проверьте корни (7.22).
Покажите, что правдоподобие (7.21) неограничено.
Покажите, что если определяется формулой (7.24), то и всюду непрерывны.
Пусть имеет дифференцируемую плотность , и пусть .
Используя интегрирование по частям, запишите знаменатель (7.27) в виде .
Покажите, что , применив неравенство Шварца к пункту (a).
Следующие три задачи посвящены исследованию технических условий, необходимых для состоятельности и асимптотической нормальности M-оценок, как отмечено в (7.26).
Для состоятельности -оценок достаточным условием является единственность и изолированность корня оценивающей функции. Установите следующую теорему.
Теорема 9.2 Предположим, что выполнены условия (A0)–(A3). Пусть — изолированный корень уравнения , где монотонна по и непрерывна в окрестности . Если — решение уравнения , то сходится к по вероятности.
Теорема 9.3 При условиях Теоремы 9.2, если, кроме того
-
конечно и отлично от нуля,
-
,
то
где . [Заметим, что это небольшое обобщение (7.27).]
Для каждой из следующих оценок выпишите определяющую её функцию и покажите, что оценка состоятельна и асимптотически нормальна при условиях Теорем 9.2 и 9.3.
Оценка наименьших квадратов, минимизирующая .
Оценка наименьших абсолютных отклонений, минимизирующая .
Усечённое среднее Хьюбера, минимизирующее (7.24).
В Примере 7.12 сравните
асимптотические распределения и ;
нормированную ожидаемую квадратичную ошибку и .
В Примере 7.12 покажите, что
и
.
В Примере 7.13 покажите, что
и независимы и имеют указанные распределения;
и — полные достаточные статистики на основе выборки из (7.33).
В Примере 7.13 определите UMVU-оценки и и асимптотические распределения этих оценок.
В предыдущей задаче сравните
асимптотическое распределение MLE и UMVU-оценки ;
нормированную ожидаемую квадратичную ошибку этих двух оценок.
В Примере 7.15,
проверьте уравнение (7.39),
покажите, что выбор даёт оценку с наилучшей эффективностью второго порядка,
покажите, что предельное отношение риска к равно 2, и
обсудите поведение этой оценки при малых объёмах выборки.
Пусть — независимые одинаково распределённые по трёхпараметрическому логнормальному распределению (7.37). Покажите, что
где
при . [Указание: (b) Для , достаточно близкого к ,
и, следовательно,
Правая часть стремится к бесконечности при (Hill 1963.]
Производные всех порядков плотности (7.37) стремятся к нулю при .
Определите предельное распределение байесовской оценки, соответствующей квадратичной функции потерь, и проверьте, что она асимптотически эффективна, в каждом из следующих случаев:
Наблюдения независимы и одинаково распределены как , где известно, а оцениваемый параметр — . Априорное распределение для — сопряжённое нормальное распределение, скажем, . (См. Пример 4.2.2.)
Наблюдения имеют гамма-распределение , оцениваемый параметр — , а имеет сопряжённую априорную плотность .
Наблюдения и априорное распределение такие же, как в Задаче 4.1.9, а оцениваемый параметр — .
Наблюдения имеют отрицательное биномиальное распределение (4.3), имеет априорную плотность , а оцениваемый параметр — (a) и (b) .
Обращаясь к Примеру 8.1, рассмотрим вместо этого минимаксную оценку параметра , заданную формулой (1.11), которая соответствует последовательности бета-априорных распределений с . Тогда,
Покажите, что предельное распределение есть , так что имеет ту же асимптотическую дисперсию, что и , но при является асимптотически смещённой.
Покажите, что ARE оценки относительно не существует, кроме случая , когда она равна 1.
Условия регулярности Теоремы 2.6 таковы: пусть независимы и одинаково распределены, каждая с плотностью относительно -конечной меры , где вещественнозначный, и предположим:
параметрическое пространство является открытым интервалом (не обязательно конечным);
распределения величин имеют общий носитель, так что не зависит от ;
для каждого плотность дважды дифференцируема по , и вторая производная непрерывна по ;
интеграл можно дважды дифференцировать под знаком интеграла;
информация Фишера удовлетворяет условию ;
для любого заданного существуют и функция (возможно, зависящие от ), такие что для всех , , и .
Если — логарифм правдоподобия, а — истинное значение, то (8.1) и (8.2) — это утверждения
где
Предположения Теоремы 2.6 влекут (8.1) и (8.2).
В Примере 8.5 апостериорная плотность после одного наблюдения равна ; это собственная плотность, и она удовлетворяет (B5) при условии .
Пусть — независимые положительные величины, каждая с плотностью , а имеет несобственную плотность . Апостериорная плотность после одного наблюдения является собственной плотностью и удовлетворяет (B5) при условии .
Приведите пример, в котором апостериорная плотность является собственной (с вероятностью 1) после двух наблюдений, но не после одного.
Докажите результат, сформулированный перед Примером 8.6.
Пусть независимы и одинаково распределены как , и рассмотрим несобственную плотность . Тогда апостериорное распределение будет несобственным при всех .
Пусть независимы и одинаково распределены с плотностью (относительно ), где вещественнозначный, а параметрическое пространство — открытый интервал; пусть обозначает истинное значение, — логарифм правдоподобия, а — score в точке . Предположим, что для любой последовательности ,
где при . Пусть
и определим
Лемма 8.7 утверждает: величина равна
Докажите Лемму 8.7.
Если и при , то , где супремум берётся по общему множеству .
Используя (a), покажите, что (8.22) и (8.23) влекут (8.21).
Условия регулярности Теоремы 2.6 таковы: пусть независимы и одинаково распределены, каждая с плотностью относительно -конечной меры , где вещественнозначный, и предположим:
параметрическое пространство является открытым интервалом (не обязательно конечным);
распределения величин имеют общий носитель, так что не зависит от ;
для каждого плотность дважды дифференцируема по , и вторая производная непрерывна по ;
интеграл можно дважды дифференцировать под знаком интеграла;
информация Фишера удовлетворяет условию ;
для любого заданного существуют и функция (возможно, зависящие от ), такие что для всех , , и .
Предположение (B1) таково: логарифм правдоподобия удовлетворяет этим условиям.
Пусть обозначает истинное значение, — score в точке , — информацию Фишера, и
Тогда (8.24) — это утверждение , а (8.26) — утверждение, что ограничена по вероятности.
Покажите, что (B1) влечёт
-
(8.24) и
-
(8.26).