Оценка параметра равномерного распределения: неправильный априор?


10

У нас есть N выборок из равномерного распределения где неизвестно. Оцените из данных.Xiθ θ[0,θ]θθ

Итак, правило Байеса ...

f(θ|Xi)=f(Xi|θ)f(θ)f(Xi)

и вероятность:

0Xiθif(Xi|θ)=i=1N1θ (редактировать: когда для всех и 0 в противном случае - спасибо whuber)0Xiθi

но без какой-либо другой информации о , кажется, что априор должен быть пропорционален (то есть равномерно) или (приор Джеффриса?) на но тогда мои интегралы не не сходятся, и я не знаю, как поступить. Любые идеи?1 1θ1 [0,]1L[0,]


2
Ваша вероятность неверна: она будет равна нулю всякий раз, когда меньше наибольшего . X iθXi
whuber

Можете ли вы показать, какие интегралы вы принимаете?

Да, так что, я думаю, я просто не знаю, как справиться с неподобающим предшественником. Например, я хочу написатьf[Xi]=Θf(Xi|θ)f(θ)dθ
будет

1
Для неправильного априора = = и для предшествующего вы аналогично получаетеПоскольку почти наверняка, есть уверенность, что интегралы сойдутся. f[Xi]=Θf(Xi|θ)f(θ)dθmax(Xi)θNdθmax(Xi)1N/(N1)f(θ)1/θmax(Xi)N/N.maxXi>0
whuber

1
Бернардо ссылается на задний план Парето - см. Каталог неинформативных априорных значений .
Стефан Лоран,

Ответы:


4

Это вызвало некоторые интересные дебаты, но обратите внимание, что это действительно не имеет большого значения для вопроса интереса. Лично я думаю, что, поскольку является параметром масштаба, аргумент группы преобразования является подходящим, приводя кθ

p(θ|I)=θ1log(UL)θ1L<θ<U

Это распределение имеет ту же форму при масштабировании задачи (вероятность также остается «инвариантной» при масштабировании). Ядро этого априора может быть получено путем решения функционального уравнения . Значения зависят от проблемы и действительно имеют значение, только если размер выборки очень мал (например, 1 или 2). Задняя часть представляет собой усеченное парето, определяемое как:f(y)=y1af(ay)=f(y)L,U

p(θ|DI)=NθN1(L)NUNL<θ<UwhereL=max(L,X(N))
Где - N-е число статистика заказа или максимальное значение выборки. Мы получаем апостериорное среднее Если мы установите и чтобы получить более простое выражение .X(N)
E(θ|DI)=N((L)1NU1N)(N1)((L)NUN)=NN1L(1[LU]N11[LU]N)
UL0E(θ|DI)=NN1X(N)

Но теперь предположим, что мы используем более общий априор, заданный (обратите внимание, что мы сохраняем пределы чтобы гарантировать, что все правильно - никаких особых математических выражений тогда ). Апостериор тогда такой же, как и выше, но с заменой на - при условии, что . Повторяя приведенные выше расчеты, мы упростили среднееp(θ|cI)θc1L,UNc+Nc+N0

E(θ|DI)=N+cN+c1X(N)

Таким образом, единообразный априор ( ) даст оценку при условии, что (среднее значение для бесконечно ). Это показывает, что спор здесь немного похож на то, использовать или нет или в качестве делителя в оценке дисперсии.c=1N1N2X(N)N2N=2NN1

Один аргумент против использования неподходящей униформы априора в этом случае состоит в том, что апостериор неправильный, когда , так как он пропорционален . Но это имеет значение только если или очень мало.N=1θ1N=1


1

Поскольку цель здесь, по-видимому, заключается в получении некоторой достоверной и полезной оценки , предварительное распределение должно соответствовать спецификации распределения населения, из которого происходит выборка. Это никоим образом не означает, что мы «вычислим» предшествующее использование самого образца - это сведет на нет действительность всей процедуры. Мы действительно знаем, что совокупность, из которой происходит выборка, представляет собой совокупность одинаковых случайных величин, каждая из которых находится в диапазоне . Это предположение поддерживается и является частью предшествующей информации, которой мы обладаем (и она не имеет ничего общего с выборкой , то есть с конкретной реализацией подмножества этих случайных переменных).θ[0,θ]

Теперь предположим, что эта совокупность состоит из случайных величин (в то время как наша выборка состоит из реализаций случайных величин). Поддерживаемое предположение говорит нам, что mn<mn

maxi=1,...,n{Xi}maxj=1,...,m{Xj}θ

Обозначим для компактности . Тогда у нас есть который также может быть записан как maxi=1,...,n{Xi}XθX

θ=cXc1

Функция плотности из из IID Равномерное с.в. в диапазоне является maxN[0,θ]

fX(x)=N(x)N1θN

для поддержки и ноль в другом месте. Затем, используя и применяя формулу изменения переменной, мы получаем предварительное распределение для , которое согласуется с поддерживаемым предположением: [0,θ]θ=cXθ

fp(θ)=N(θc)N1θN1c=NcNθ1θ[x,]

что может быть неправильно, если мы не укажем константу соответствующим образом. Но наш интерес заключается в том, чтобы иметь правильную апостериор для , а также мы не хотим ограничивать возможные значения (за исключением ограничения, подразумеваемого поддерживаемым предположением). Таким образом , мы выходим неопределенными. Затем, написав апостериорныйcθθc
X={x1,..,xn}

f(θX)θNNcNθ1f(θX)=ANcNθ(N+1)

для некоторой нормализующей константы A. Мы хотим

Sθf(θX)dθ=1xANcNθ(N+1)dθ=1

ANcN1NθN|x=1A=(cx)N

Вставка в заднюю часть

f(θX)=(cx)NNcNθ(N+1)=N(x)Nθ(N+1)

Обратите внимание, что неопределенная константа предыдущего распределения удобно удаляется.c

Апостериор суммирует всю информацию, которую конкретный образец может дать нам относительно значения . Если мы хотим получить конкретное значение для мы можем легко рассчитать ожидаемое значение апостериора, θθ

E(θX)=xθN(x)Nθ(N+1)dθ=NN1(x)NθN+1|x=NN1x

Есть ли интуиция в этом результате? Что ж, по мере того, как число увеличивается, более вероятно, что максимальная реализация среди них будет все ближе и ближе к их верхней границе, - что в точности соответствует последнему среднему значению : если, скажем, , , но если , Это показывает, что наша тактика в отношении выбора предыдущего была разумной и соответствовала рассматриваемой проблеме, но не обязательно «оптимальной» в некотором смысле.XθθN=2E(θX)=2xN=10E(θX)=109x


1
Основываясь на предварительных данных звучит подозрительно для меня. Как вы обосновываете такой подход?
whuber

2
Я ничего не имею против того, что ваш предшественник не "лучший". Где я сказал что-то подобное? Я просто пытаюсь понять твой подход. Я еще не понимаю это равенство. Если постоянна в равенстве , значит ли это, что и неслучайны? Кстати, вы не используете тот факт, что при выводе предыдущего, не так ли? (cc @whuber)cθ=cXXθc1
Стефан Лоран

1
А поддержка твоего приора зависит от данных? ( )θ[x,[
Стефан Лоран,

3
Априорная зависимость (даже если это только через поддержку) от данных звучит неправильно: вы не можете знать максимум сэмпла до того, как сэмпл был сгенерирован . Более того, вы утверждаете, что является почти верным равенством, причем как и random (таким образом, существует корреляция ). Но это подразумевает, что апостериорное распределение (которое является условным распределением данного образца) является массой Дирака в . И это противоречит вашему выводу апостериорного распределения. ... (персонажей не осталось ...)θ=cXθX1θθcx
Стефан Лоран,

1
Последующее распределение - это Дирак в точке означает, что - это . Теорема Байеса не является причиной. Вы уничтожаете все, предполагая . Это означает, что , поэтому условное распределение заданное является массой Дирака в , тогда как исходное предположение состоит в том, что это распределение является равномерным распределением на . θcxθ cxθ=cXX=θ/cXθθ/c(0,θ)
Стефан Лоран

0

Теорема о равномерном априорном распределении (интервальный случай):

«Если вся информация о внешняя по отношению к данным , захвачена одним предложением тогда Ваша единственная возможная логически-внутренне согласованная предварительная спецификация - это θD

B={{Possible values for θ}={the interval (a,b)},a<b}
f(θ)=Uniform(a,b)

Таким образом, ваша предварительная спецификация должна соответствовать предварительной версии Джеффри, если вы действительно верите в приведенную выше теорему ".

Не является частью теоремы о равномерном предыдущем распределении:

В качестве альтернативы вы можете указать свое предыдущее распределение как распределение Парето, которое является сопряженным распределением для униформы, зная, что ваше последующее распределение должно быть другим равномерным распределением по сопряженности. Однако, если вы используете дистрибутив Pareto, вам нужно будет каким-то образом указать параметры распределения Pareto.f(θ)


4
Сначала вы говорите, что «единственно возможным логически внутренне согласованным» ответом является равномерное распределение, а затем вы приступаете к предложению альтернативы. Это звучит нелогично и противоречиво для меня :-).
whuber

2
Я не могу согласиться Например, также является множествомКогда PDF для равен для . Но согласно "теореме", , pdf которого равен в этом интервале. Короче говоря, хотя утверждение не зависит от того, как проблема параметризована, вывод «теоремы» зависит от параметризации, поэтому она неоднозначна. B{θ|θ3(a3,b3)}.ΘUniform(a,b),Ψ=Θ31/(3ψ2/3(ba))a3<ψ<b3ΨUniform(a3,b3)1/(b3a3)
whuber

2
BabakP: Как можно сказать, что это теорема ? Теорема - это математическое утверждение с математическим доказательством. Эту «теорему» более целесообразно назвать «принципом», но она не имеет смысла, потому что она противоречива, как показывает @whuber.
Стефан Лоран

2
Спасибо за ссылку BabakP. Я хотел бы отметить, что «эскиз доказательства» является поддельным. Дрейпер делит интервал на конечное число одинаково расположенных значений и «переходит к пределу». Любой может разделить интервал на значения, отстоящие друг от друга для аппроксимации любой плотности, которую они любят, и аналогичным образом переходят к пределу, создавая совершенно произвольные «только возможные логически-внутренне согласованные предыдущие спецификации». Подобные вещи - а именно, использование плохой математики в попытке показать, что небайесовцы нелогичны - дают байесовскому анализу (незаслуженно) дурную славу. (Копия @ Стефан.)
whuber

1
@ Стефан Пожалуйста, прости мою нечувствительность ( insensibilité ) - я восхищаюсь твоим умением общаться здесь на другом языке и не использую заведомо неясные термины! Bogus - это прилагательное, которое происходит от 200-летнего американского сленгового термина, обозначающего машину для подделки денег. В данном случае это математическая машина для подделок теорем :-).
whuber
Используя наш сайт, вы подтверждаете, что прочитали и поняли нашу Политику в отношении файлов cookie и Политику конфиденциальности.
Licensed under cc by-sa 3.0 with attribution required.