Чи існує сполучник до розподілу Лапласа?


13

Чи існує сполучник до розподілу Лапласа ? Якщо ні, то чи відомий вираз закритої форми, який наближає до заднього для параметрів розподілу Лапласа?

Я гуляв досить багато, не маючи успіху, тому моє теперішнє здогадка "ні" на вищезазначені питання ...


1
"Середні дисперсії середніх дисперсій серед полісонів та скоттів" - це дасть вам приблизний рівень використання MAP за допомогою алгоритму em.
ймовірністьлогічний

Відповіді:


12

Давайте спочатку розглянемо їх по черзі (приймаючи інше як дане).

За посиланням (з модифікацією відповідно до конвенції про використання грецьких символів для параметрів):

f(x|μ,τ)=12τexp(|xμ|τ)

- параметр масштабу :

L(τ)τk1eSτ

при певних значеннях і . Тобто вірогідність зворотного гамма-форми.kS

Таким чином, параметр шкали має кон'югат попереднього - шляхом огляду попереднього сполучення є зворотною гаммою.

- параметр розташування

Це, дійсно, складніше, тому щоне спрощується в щось зручне в ; Я не думаю, що існує спосіб "зібрати терміни" (ну і таким чином, але нам все одно не потрібно).i|xiμ|μ

Рівномірний пріоритет просто уріже задню частину, що не так вже й погано працювати, якщо це здається правдоподібним як попереднє.

Цікавою можливістю, яка інколи може бути корисною, є досить легко включити попередній Лаплас (той, що має ту саму шкалу, що і дані), використовуючи псевдоспостереження. Можна також наблизити якусь іншу (більш жорстку) попередню за допомогою декількох псевдоспостережень)

Насправді, щоб узагальнити це, якби я працював з Лапласом, я б спокусився просто узагальнити від постійної ваги-постійної ваги до роботи із зваженою версією спостереження Лапласа (що рівно потенційно відрізняється шкалою для кожна точка даних) - імовірність журналу все ще є лише безперервною кусочно-лінійною функцією, але нахил може змінюватися на нецілі суми в точках приєднання. Тоді існує зручне попереднє "сполучення" - просто ще один "зважений" Лаплас або, справді, будь-яка форма абоexp ( - j w j | μ - θ j | )exp(j|μθj|/ϕj)exp(jwj|μθj|)(хоча його потрібно буде відповідним чином масштабувати, щоб зробити фактичну щільність) - дуже гнучка сім'я розподілів, яка, очевидно, призводить до задньої "тієї ж форми", що ймовірність зваженого спостереження, і щось легко працювати з малювати; насправді навіть псевдоспостереження все ще працює.

Він також досить гнучкий, щоб його можна було використовувати для наближення інших пріорів.

(Більш загальне все-таки, можна було б працювати в масштабі журналу і використовувати безперервний, кусково-лінійний логаритмічний увігненний до цього, а задній також буде такої форми; це має включати асиметричний Лаплас як особливий випадок)

Приклад

Тільки щоб показати, що з цим досить легко впоратися - нижче - пріоритетний (пунктирний сірий), вірогідність (пунктирний, чорний) та задній (суцільний, червоний) для параметра розташування зваженого Лапласа (... це було з відомими масштабами ).

введіть тут опис зображення

Зважений підхід Лапласа добре би працював у MCMC, я думаю.

-

Цікаво, чи отриманий задній режим є середньозваженою?

- насправді (щоб відповісти на моє власне запитання), схоже, що відповідь на це «так». Це робить досить приємно працювати.

-

Спільний пріоритет

Очевидним підходом було б написати : було б відносно легко мати у тому ж вигляді, що і вище - де може бути коефіцієнтом масштабування на попередній, тому попередній буде вказаний відносно - а потім обернена гамма до , безумовно.μ | τ τ τ τf(μ,τ)=f(μ|τ)f(τ)μ|ττττ

Безсумнівно, щось більш загальне для спільного початку є цілком можливим, але я не думаю, що я буду продовжувати спільну справу далі, ніж тут.

-

Я ніколи не бачив і не чув цього попереднього підходу із зваженим рівнем, але придумати його було досить просто, так що, мабуть, це вже було зроблено. (Довідники вітаються, якщо хтось знає про що-небудь.)

Якщо взагалі ніхто не знає жодних посилань, можливо, я повинен щось написати, але це було б дивним.


Вау, чудова відповідь. Я впевнений, не знаю жодних посилань на щось подібне. Якщо ви щось знайдете або щось написати, будь ласка, повідомте мене!
Rasmus Bååth

1
Один з можливих способів отримати параметр розташування - це використовувати звичайне представлення дисперсійної суміші лапласа. Це умовно сполучається попереднє, хоча ...
ймовірність відомості

@probabilityislogic, що цікаво. У попередніх редакціях я вказував рядок, вказуючи на те, що Лаплас - це експоненціальна суміш нормалів, тому що я замислювався, чи можна з цим щось зробити, але як редагував відповідь далі, він більше нікуди не підходив, і я взяв це знову. З вашого корисного коментаря звучить так, що його можна використовувати таким чином; це, ймовірно, стане в нагоді.
Glen_b -Встановити Моніку
Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.