特定のMLE(マルコフ連鎖)の対数尤度の計算


9

私は現在、マルコフチェーンを使用しており、いくつかのソース(つまり、aからbへの遷移の数をaから他のノードへの全体的な遷移の数で割ったもの)によって提案された遷移確率を使用して最尤推定を計算しました。

次に、MLEの対数尤度を計算します。


遷移確率の最尤推定値を既に計算しているので、何の対数尤度を正確に計算しますか?
Nick

回答:


11

ましょうマルコフ連鎖の経路であるとletパスを観察する確率である場合 IS真のパラメーター値(別名尤度関数)条件付き確率の定義を使用すると、 P θX 1X Tθ θ{バツ}=1TPθバツ1バツTθθ

Pθバツ1バツT=PθバツT|バツT1バツ1Pθバツ1バツT1

これはマルコフ連鎖であるため、であることがこれを簡素化しますPθバツT|バツT1バツ1=PθバツT|バツT1

Pθバツ1バツT=PθバツT|バツT1Pθバツ1バツT1

これと同じロジックを回繰り返すと、T

Pθバツ1バツT=Π=1TPθバツ|バツ1

ここで、はプロセスの初期状態として解釈されます。右側の項は、遷移行列の単なる要素です。あなたが要求した対数尤度だったので、最終的な答えは次のとおりです:バツ0

Lθ=Σ=1TログPθバツ|バツ1

これは単一のマルコフチェーンの可能性です。データセットに複数の(独立した)マルコフチェーンが含まれている場合、完全な可能性はこの形式の項の合計になります。


うわー、答えてくれてありがとう。この場合、はMLEから取得された「遷移」確率ですよね?Pθ
fsociety 2013年

@ph_singer、どういたしまして。は、パラメーター値与えられた場合に、状態からに移動する確率です。遷移行列に構造を課さない場合(これはそのように聞こえます)、は遷移確率のベクトルを示すだけです(そしてMLEは、質問文で正しく示したように、単なるサンプルの比率です)。 : は、状態で終了した状態からの移動のサンプル比率になります。X I - 1 X I θ θ P θ M L EX I | X I - 1X I - 1 X IPθバツ|バツ1バツ1バツθθPθ^MLEバツ|バツ1バツ1バツ
Macro

再度、感謝します!もう1つの質問:別の次数(たとえば、k = 2)を使用する場合、このプロセスはどのように機能しますか?
学会2013年

「注文」の意味を教えてください。
マクロ

(+1)OPはおそらく、を意味し、2次 MC を示します。つまり、直前の2つの状態に応じて、最新の状態。X i 1X i 2 X i 1k=2バツ1バツ2バツ1
枢機卿
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.