Примеры того, когда доверительный интервал и вероятный интервал совпадают


11

В статье в Википедии о Credible Interval говорится:

Для случая одного параметра и данных, которые могут быть сведены в одну достаточную статистику, можно показать, что вероятный интервал и доверительный интервал совпадут, если неизвестный параметр является параметром местоположения (т.е. функция прямой вероятности имеет вид Pr (x | μ) = f (x - μ)) с априором, который является равномерным плоским распределением; [5], а также, если неизвестный параметр является масштабным параметром (то есть прямая функция вероятности имеет вид Pr (x) | s) = f (x / s)), с предшествующим Джеффрисом [5] - последующим последующим, потому что при взятии логарифма такого масштабного параметра он превращается в параметр местоположения с равномерным распределением. Но это явно особые (хотя и важные) случаи; в общем, такая эквивалентность невозможна ».

Могут ли люди привести конкретные примеры этого? Когда 95% CI фактически соответствует «вероятности 95%», таким образом, «нарушая» общее определение CI?

Ответы:


13

нормальное распределение:

Возьмите нормальное распределение с известной дисперсией. Мы можем принять эту дисперсию равной 1, не теряя общности (просто разделив каждое наблюдение на квадратный корень из дисперсии). Это имеет распределение выборки:

п(Икс1,,,ИксN|μ)знак равно(2π)-N2ехр(-12Σязнак равно1N(Икся-μ)2)знак равноAехр(-N2(Икс¯-μ)2)

Где - это константа, которая зависит только от данных. Это показывает, что выборочное среднее является достаточной статистикой для среднего населения. Если мы используем равномерный априор, то апостериорное распределение для µ будет:Aμ

(μ|Икс1,,,ИксN)~NормaL(Икс¯,1N)(N(μ-Икс¯)|Икс1,,,ИксN)~NормaL(0,1)

Таким образом, вероятный интервал будет иметь вид:1-α

(Икс¯+1NLα,Икс¯+1NUα)

Где и U α выбираются так, чтобы стандартная нормальная случайная величина Z удовлетворяла:LαUαZ

пр(Lα<Z<Uα)знак равно1-α

Теперь мы можем начать с этой «ключевой величины» для построения доверительного интервала. Распределение выборки при фиксированномцявляется стандартным нормальным распределением, таким образоммы можем подставить это в приведенном выше вероятность:N(μ-Икс¯)μ

пр(Lα<N(μ-Икс¯)<Uα)знак равно1-α

Затем переупорядочить, чтобы решить для , и доверительный интервал будет таким же, как доверительный интервал.μ

Параметры шкалы:

Для параметров масштаба файлы pdf имеют вид . Мы можем взять(Xi|s)Uniform(0,s), что соответствуетf(t)=1. Совместное распределение выборки:п(Икся|s)знак равно1sе(Иксяs)(Икся|s)~UNяеорм(0,s)е(T)знак равно1

п(Икс1,,,ИксN|s)знак равноs-N0<Икс1,,,ИксN<s

ИксмaИкс

пр(ИксмaИкс<Y|s)знак равнопр(Икс1<Y,Икс2<Y,,,ИксN<Y|s)знак равно(Ys)N

Yзнак равноQsQзнак равноs-1ИксмaИкспр(Q<Q)знак равноQNбеTa(N,1)Lα,Uα

пр(Lα<Q<Uα)знак равно1-αзнак равноUαN-LαN

И мы подставляем основное количество:

пр(Lα<s-1ИксмaИкс<Uα)знак равно1-αзнак равнопр(ИксмaИксLα-1>s>ИксмaИксUα-1)

И есть наш доверительный интервал. Для байесовского решения с джефри до этого имеем:

п(s|Икс1,,,ИксN)знак равноs-N-1ИксмaИкср-N-1dрзнак равноN(ИксмaИкс)Ns-N-1
пр(s>T|Икс1,,,ИксN)знак равноN(ИксмaИкс)NTs-N-1dsзнак равно(ИксмaИксT)N

Теперь мы подключаем доверительный интервал и рассчитываем его достоверность.

пр(ИксмaИксLα-1>s>ИксмaИксUα-1|Икс1,,,ИксN)знак равно(ИксмaИксИксмaИксUα-1)N-(ИксмaИксИксмaИксLα-1)N

знак равноUαN-LαNзнак равнопр(Lα<Q<Uα)

1-α


Шедевр, спасибо! Я надеялся, что может быть такой ответ: «при расчете среднего значения выборки из нормального распределения 95% -й доверительный интервал на самом деле также является 95-процентным доверительным интервалом» или что-то простое в этом роде. (Просто составив этот предполагаемый ответ, я не имею ни малейшего понятия о конкретных примерах.)
Уэйн

Я считаю, что частый 95% интервал прогнозирования / толерантности соответствует байесовскому интервалу прогнозирования с регрессией МНК и нормальными ошибками. Похоже, что когда я сравниваю прогнозируемый ответ с имитированным ответом. Это правда?
Уэйн

Yзнак равноα+βИксα,βσ

Большое спасибо! Я пытался объяснить CI для регрессии, которую я сделал, с точки зрения доверительного интервала, и он просто не связан с непрофессионалом, который ожидает Credible Interval. Делает жизнь намного проще для меня ... хотя, возможно, это плохо для всего статистического мира, так как это укрепит непонимание непрофессионала с КИ.
Уэйн

@Wayne - ситуация немного более общая, чем просто семейства масштаба. Обычно CI будет эквивалентен вероятному интервалу, если он основан на «достаточной статистике» (как это было два), где он существует. Если нет достаточной статистики, тогда CI должен обусловить то, что называется «вспомогательной статистикой», чтобы иметь достоверную интервальную интерпретацию.
вероятностная
Используя наш сайт, вы подтверждаете, что прочитали и поняли нашу Политику в отношении файлов cookie и Политику конфиденциальности.
Licensed under cc by-sa 3.0 with attribution required.