Допустимость и минимаксность в экспоненциальных семействах
[25/0%]Лемма 2.1 была расширена Berger (1990a) для случая, когда оцениваемая функция не обязана быть ограничена конечным интервалом, а вместо этого достигает максимума или минимума в конечной точке параметра.
Лемма 8.1. Пусть оцениваемая функция непостоянна и достигает глобального максимума или минимума в точке , для которой п.в. (относительно доминирующей меры ), и пусть функция потерь удовлетворяет предположениям Леммы 2.1. Тогда любая оценка , принимающая значения выше максимума или ниже минимума, недопустима.
Покажите, что если минимизирует и является несмещённой оценкой , то существует такое, что множество удовлетворяет . Аналогичный вывод верен, если является максимумом.
Предположим, что является минимумом. (Случай максимума рассматривается аналогично.) Покажите, что оценка
удовлетворяет .
Для ситуации из Примера 2.3 примените Лемму 8.1, чтобы установить недопустимость UMVU-оценки . Также объясните, почему предположения Леммы 8.1 не выполняются для оценивания .
Определите байесовский риск оценки (2.4), когда имеет априорное распределение .
Докажите утверждение (d) во втором доказательстве Примера 2.8: что существует последовательность значений , для которой .
Покажите, что оценка величины недопустима (при квадратичной функции потерь) при каждом из следующих условий:
если для всех и ;
если для всех и .
Покажите, что оценка величины недопустима (при квадратичной функции потерь) при каждом из следующих условий:
если и ,
если и .
Покажите, что если , то оценка недопустима (при квадратичной функции потерь) при каждом из следующих условий:
если для всех и ;
если для всех и (Gupta 1966).
Brown (1986a) указывает на связь между информационным неравенством и несмещённой оценкой риска оценок стейновского типа.
Покажите, что из (2.7) следует
и, следовательно, если , то .
Покажите, что нетривиальное решение привело бы к улучшенной оценке при в Следствии 4.7.2.
Функция плотности называется понижающей вариацию порядка , если для любой функции с переменами знака (без учёта нулей) математическое ожидание имеет не более перемен знака. Если имеет ровно перемен знака, они происходят в том же порядке. Покажите, что обладает свойством тогда и только тогда, когда она имеет монотонное отношение правдоподобия. (См. TSH2, Лемма 2, Раздел 3.3, для доказательства импликации «если».) Brown и др. (1981) дают подробное введение в эту тему, включая характеризации для многих семейств распределений (экспоненциальное семейство обладает свойством , как и с параметром , и нецентральное по параметру ). Существует эквивалентность между и — полной положительностью порядка по Karlin (1968), а именно .
Для ситуации из Примера 2.9 покажите, что:
без потери общности ограничение можно свести к , .
Если — априорное распределение, помещающее массу 1/2 в каждую из точек , то байесовская оценка при квадратичной функции потерь есть
Для
и, следовательно, по Следствию 1.6, является минимаксной.
Для оценка из пункта (b) больше не является минимаксной. Объясните, почему это так, и предложите альтернативную оценку в этом случае.
Для ситуации из Примера 2.10 покажите, что:
.
Оценка с является линейной минимаксной оценкой для всех с минимаксным риском .
является линейной минимаксной оценкой для .
Предположим, что имеет распределение , а имеет распределение , где и изменяются независимо. Если известно, что , то любая оценка может быть улучшена с помощью
В Примере 2.13 докажите, что оценка недопустима при .
Пусть независимые одинаково распределённые случайные величины с плотностью , и пусть . Нас интересует оценивание при квадратичной функции потерь с помощью линейных оценок , где и — константы. Покажите, что:
допустимость оценки в Примере 2.13 эквивалентна допустимости при подходящем выборе и .
риск задаётся формулой
при , когда , и, следовательно, оценка в Примере 2.13 недопустима при .
Это упражнение иллюстрирует тот факт, что константы не обязательно являются допустимыми оценками.
Для ситуации из Примера 2.15 пусть .
Покажите, что риск оценки при квадратичной функции потерь минимизируется при выборе
Stein (1964) показал, что для всех . Предполагая это, выведите, что доминирует над при квадратичной функции потерь, где
Покажите, что , и, следовательно, задаётся формулой (2.31).
Наилучшая эквивариантная оценка вида была получена Brewster и Zidek (1974) и задаётся формулой
где математическое ожидание вычисляется в предположении и . Покажите, что является обобщённой байесовской оценкой относительно априорного распределения
Brewster и Zidek первоначально получили свою оценку не как байесовскую, а, напротив, сначала нашли оценку, а затем нашли априорное распределение. Brown (1968) рассмотрел семейство оценок, аналогичных оценкам Stein (1964), которые принимали различные значения в зависимости от порогового значения для . Brewster и Zidek (1974) показали, что число пороговых значений может быть сколь угодно большим. Они построили последовательность оценок с убывающими рисками и всё более плотными пороговыми значениями, предел которой являлся наилучшей эквивалентной оценкой.
Покажите эквивалентность следующих соотношений:
(2.26) и (2.27),
(2.34) и (2.35) при , и
(2.38) и (2.39).
В Примере 2.17 покажите, что оценка недопустима для всех вне треугольника (2.39).
Докажите допустимость оценок, соответствующих внутренности треугольника (2.39), применив Теорему 2.4 и используя результаты Примера 4.1.5.
Используя Теорему 2.14, приведите альтернативное доказательство допустимости оценки , удовлетворяющей (2.6), в Примере 2.5.
Определите, какие оценки допустимы для оценивания при квадратичной функции потерь в следующих ситуациях:
имеет распределение Пуассона.
имеет отрицательное биномиальное распределение (Gupta 1966).
Пусть имеет распределение , и рассмотрим оценивание при функции потерь с ограничением , где известно.
Используя рассуждение, аналогичное рассуждению из Примера 2.9, покажите, что не является минимаксной оценкой и что наименее благоприятное априорное распределение должно иметь множество [из (1.5)], состоящее из конечного числа точек.
Пусть — априорное распределение, помещающее массу , в точки параметра . Покажите, что байесовская оценка, соответствующая этому априорному распределению, есть
Пусть — решение уравнения . Покажите, что при , , одноточечное априорное распределение даёт минимаксную оценку. Вычислите минимаксный риск и сравните его с риском .
Пусть — первый положительный нуль уравнения , где — двухточечное априорное распределение . Покажите, что при , , двухточечное априорное распределение даёт минимаксную оценку (используйте Следствие 1.6). Вычислите минимаксный риск и сравните его с риском .
С ростом ситуация усложняется, и точные минимаксные решения становятся труднодостижимыми. В этих случаях линейные приближения могут быть вполне удовлетворительными. См. Johnstone и MacGibbon 1992, 1993.
Покажите, что условия (2.41) и (2.42) из Примера 2.22 являются не только достаточными, но и необходимыми для допустимости (2.40).
Пусть и независимо распределены по закону Пуассона с и соответственно. Покажите, что допустима для оценивания при квадратичной функции потерь тогда и только тогда, когда либо , либо (Makani 1972).
Пусть имеет плотность .
Покажите, что .
Покажите, что допустима для оценивания при квадратичной функции потерь тогда и только тогда, когда .
Пусть распределена как , и пусть имеет несобственную априорную плотность . При квадратичной функции потерь формальная байесовская оценка равна , и она не является ни минимаксной, ни допустимой. (См. также Задачу 2.15.) Условия, при которых формальная байесовская оценка, соответствующая несобственному априорному распределению для в Примере 3.4, является допустимой, даны в Zidek (1970).
Покажите, что естественное параметрическое пространство семейства (2.16) есть для нормального (с известной дисперсией), биномиального распределений и распределения Пуассона, но не в случае гамма-распределения или отрицательного биномиального распределения.