ラプラス分布の前に共役が存在しますか?


13

ラプラス分布の共役事前分布はありますか?そうでない場合、ラプラス分布のパラメーターの事後を近似する既知の閉形式表現はありますか?

私は多くのことをグーグルで検索して成功しなかったので、現在の推測は上記の質問に「いいえ」です...


1
Googleの「polsonとscottの正規分散の平均混合」-これにより、emアルゴリズムを介してMAPを使用して、おおよそのベイが得られます。
確率論的

回答:


12

最初にそれらを1つずつ見ていきましょう(もう一方を与えられたとおりに取ります)。

リンクから(パラメーターにギリシャ語の記号を使用する規則に従うように変更して):

fバツ|μτ=12τexp|バツμ|τ

- スケールパラメーター

Lττk1eSτ

およびSの特定の値に対して。つまり、尤度は逆ガンマ形式です。kS

したがって、スケールパラメーターには共役事前分布があります。検査により、共役事前分布は逆ガンマになります。

- 位置パラメータ

|バツμ|μ

均一な事前分布は、単純に事後を切り捨てます。これは、事前分布として妥当であると思われる場合には、それほど悪いことではありません。

時折役に立つかもしれない興味深い可能性の1つは、擬似観測を使用することにより、ラプラスの事前データ(データと同じスケールを持つもの)を含めるのがかなり簡単なことです。また、いくつかの疑似観測を介して、事前に他の(より厳密な)近似を行う場合もあります)

expj|μθj|/ϕjexpjwj|μθj|

また、他の事前分布を近似するために使用できるほど十分に柔軟です。

(より一般的には、対数スケールで作業し、連続的で区分的に線形の対数凹を使用することができ、前も後もその形式になります。これには、特殊なケースとして非対称ラプラスが含まれます)

対処が非常に簡単であることを示すために、以下は加重ラプラスの位置パラメーターの事前(点線の灰色)、尤度(破線、黒)および事後(実線、赤)です(...これは既知のスケールでした)。

ここに画像の説明を入力してください

重み付きラプラスアプローチはMCMCでうまく機能すると思います。

-

結果の事後モードは加重中央値であるのだろうか?

-実際には(自分の質問に答えるために)、それに対する答えは「はい」のように見えます。それで作業するのはかなり良いです。

-

共同事前

fμτ=fμ|τfτμ|ττττ

ジョイントの事前のより一般的なことは間違いなく可能ですが、ここよりもジョイントケースを追求するつもりはありません。

-

以前にこの重み付きラプラスのアプローチを目にしたことも聞いたこともありませんでしたが、考え出すのはかなり簡単だったので、おそらくすでに行われています。(誰かが知っているなら、参照は大歓迎です。)

誰も参照をまったく知らない場合は、おそらく何かを書く必要がありますが、それは驚くべきことです。


うわー、素晴らしい答え。私は確かに同様のものへの参照を知らない。何かを見つけたり、何かを書いたりしたら、私に知らせてください!
ラスマスバース14

1
locationパラメーターを取得する1つの方法は、laplaceの通常の分散混合表現を使用することです。ただし、これは事前に条件付き共役です...
potentialislogic

@probabilityislogicそれは興味深いです。以前の編集では、ラプラスは法線の指数スケールの混合物であったことを指摘する行を入れましたが、それで何かできるかどうか疑問に思っていましたが、答えをさらに編集すると、どこにも収まらず、再びそれ。有益なコメントから、そのように使用できるように思えます。それは便利だろう。
-Glen_b-モニカの復活14
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.