Чи визначають два квантування бета-розподілу його параметри?


9

Якщо я даю два та їх відповідні місця (кожен) у відкритому інтервалі , чи завжди я можу знайти параметри бета-розподілу, які мають ці кванти у зазначених місцях?(q1,q2)(l1,l2)(0,1)


1
Ні, основний контрприклад (q1, q2) = (0,1) і (l1, l2) = (0,1) незалежно від параметрів.
Тім

1
@Tim Я думаю, що я бачу вашу думку, але ваш контрприклад не відповідає умовам, які я вказав (наприклад, що розташування знаходяться у відкритому інтервалі (0,1)).
Бота

1
Я думаю, ви можете це зробити чисельно (і щоб було унікальне рішення), але це зажадає трохи зусиль.
Glen_b -Встановіть Моніку

1
Я теж думаю - чисельне вирішення не складно, але знайти аргумент для унікальності непросто.
Елвіс

1
@Elvis насправді, я підозрюю, що може бути спосіб це зробити, переглянувши логи обох змінних (ОП l і q).
Glen_b -Встановити Моніку

Відповіді:


9

Відповідь "так", якщо дані відповідають очевидним вимогам узгодженості. Аргумент простий, заснований на простому побудові, але він потребує певного налаштування. Це зводиться до інтуїтивно привабливого факту: збільшення параметраa у бета-версії(a,b) розподіл збільшує значення його щільності (PDF) більше для більшого x ніж менший x; і збільшуєтьсяb робить навпаки: менший x є, чим більше збільшується значення PDF.

Деталі випливають.


Нехай бажане q1 квантил бути x1 і бажане q2 квантил бути x2 з 1>q2>q1>0 і таким чином) 1>x2>x1>0. Тоді є унікальніa і b для яких бета(a,b) розподіл має ці кванти.

Складність з демонстрацією цього полягає в тому, що бета-розподіл передбачає непокірну постійну нормалізацію. Нагадаємо визначення: дляa>0 і b>0, бета-версія(a,b) розподіл має функцію щільності (PDF)

f(x;a,b)=1B(a,b)xa1(1x)b1.

Константа нормалізації - функція Бета

B(a,b)=01xa1(1x)b1dx=Γ(a)Γ(b)Γ(a+b).

Все стає безладним, якщо ми спробуємо розмежуватись f(x;a,b) безпосередньо стосовно a і b, що було б найгрубішим способом спроби демонстрації.

Один із способів уникнути необхідності аналізу бета-функції - зазначити, що квантили відносні області. Це є,

qi=F(xi;a,b)=0xif(x;a,b)dx01f(x;a,b)dx

для i=1,2. Ось, наприклад, функція PDF та накопичувальної дистрибуції (CDF)F бета-версії(1.15,0.57) розповсюдження для якого x1=1/3 і q1=1/6.

Фігура 1

Функція щільності xf(x;a,b) зображено зліва. q1- площа під кривою ліворуч відx1, показана червоним кольором відносно загальної площі під кривою. q2 - область зліва від x2, що дорівнює сумі червоної та синьої областей, знову ж таки відносно загальної площі . CDF праворуч показує, як(x1,q1) і (x2,q2) позначте на ній дві чіткі точки.

На цій фігурі (x1,q1) було зафіксовано в (1/3,1/6), a було обрано 1.15, а потім значення b було знайдено за що (x1,q1) лежить на бета-версії(a,b) CDF.

Лема : Такаb завжди можна знайти.

Щоб бути конкретним, нехай (x1,q1)фіксуватися раз і назавжди. (Вони залишаються однаковими на наведених нижче ілюстраціях: у всіх трьох випадках відносна область зліва відx1 дорівнює q1.) Для будь-якого a>0, лема стверджує, що існує унікальне значення b, написана b(a), для котрого x1 є q1 квантил бета(a,b(a)) розповсюдження.

Щоб зрозуміти, чому, зверніть увагу, що спочатку b наближається до нуля, вся ймовірність накопичується поблизу значень 0, звідки F(x1;a,b) підходи 1. Якb наближається до нескінченності, вся ймовірність накопичується поблизу значень 1, звідки F(x1;a,b) підходи 0. Між тим функція bF(x1;a,b) суворо зростає в Росії b.

Це твердження є геометрично очевидним: воно означає, що якщо ми подивимось на ділянку зліва під кривою xxa1(1x)b1 відносно загальної площі під кривою і порівняйте її з відносною площею під кривоюxxa1(1x)b1 для b>b, то остання площа порівняно більша. Співвідношення цих двох функцій є(1x)bb. Це функція, рівна1 коли x=0, неухильно опускаючись до 0 коли x=1. Тому висоти функції xf(x;a,b)є відносно більше , ніж висотиxf(x;a,b) для x зліва від x1 ніж вони за x праворуч від x1. Отже, область ліворуч відx1в першому повинен бути порівняно більшим за площу праворуч відx1. (Це просто перекласти в суворий аргумент, використовуючи, наприклад, суму Рімана.)

Ми бачили, що функція bf(x1;a,b) строго монотонно зростає з обмежуючими значеннями при 0 і 1 як b0 і b,відповідно. Він також (чітко) безперервний. Отже, існує числоb(a) де f(x1;a,b(a))=q1 і це число унікальне, що підтверджує лему.

Цей же аргумент показує, що як b збільшується, площа зліва від x2збільшується. Отже, значенняf(x2;a,b(a)) діапазон через деякий інтервал чисел як a прогресує майже 0 майже . Межа f(x2;a,b(a)) як a0 є q1.

Ось приклад, де a близький до 0 (це дорівнює 0.1). Зx1=1/3 і q1=1/6 (як на попередньому малюнку), b(a)0.02. Між ними майже немає ділянки x1 і x2:

Малюнок 2

CDF практично рівна між ними x1 і x2, звідки q2 практично на вершині q1. У ліміті як a0, q2q1.

З іншого боку, досить великі значення a призводить до F(x2;a,b(a)) довільно близький до 1. Ось приклад з (x1,q1) як і раніше.

Малюнок 3

Ось a=8 і b(a) майже 10. Тепер F(x2;a,b(a)) є по суті 1: Праворуч від району майже немає x2.

Отже, ви можете вибрати будь-яку q2 між q1 і 1 і відрегулювати a до F(x2;a,a(b))=q2. Як і раніше, це aповинен бути унікальним, QED .


Робочий Rкод для пошуку рішень розміщений у розділі Визначення параметрів бета-розподілуα і βз двох довільних точок (квантилів) .


Ця відповідь показує, що якщо ми обираємо фіксовану a або bми знайдемо унікальне відповідне значення. Можна було б побудувати функції, які мають фіксовану область у[0,x1], [x1,x2] і [x2,1]. Я не відразу бачу, чому це гарантувало б, що набірα і βє унікальним. Чи бажаєте ви розробити і просвітити мене?
Ян

@Jan Не могли б пояснити, що ви маєте на увазі під набором " α і β"? Ці символи ніде не з’являються в цій нитці.
whuber
Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.