Verilen MLE (Markov Zincirleri) için günlük olasılığını hesaplama


9

Şu anda Markov zincirleri ile çalışıyorum ve birkaç kaynak tarafından önerilen geçiş olasılıklarını kullanarak Maksimum Olabilirlik Tahmini'ni hesapladım (yani, a'dan b'ye geçişlerin sayısının a'dan diğer düğümlere toplam geçiş sayısına bölünmesiyle).

Şimdi MLE'nin log olasılığını hesaplamak istiyorum.


Geçiş olasılıklarının maksimum olabilirlik tahminini zaten hesapladınız ve şimdi tam olarak neyin log olabilirliğini hesaplamak istiyorsunuz?
Nick

Yanıtlar:


11

Let markov zincirinin bir yol olabilir ve izin yolu gözlemleme ihtimali olduğu zaman isimli gerçek parametre değeri (aka için olasılık işlevi ). Koşullu olasılık tanımını kullanarak biliyoruz{Xi}i=1TPθ(X1,...,XT)θθ

Pθ(X1,...,XT)=Pθ(XT|XT1,...,X1)Pθ(X1,...,XT1)

Bu bir Markov zinciri olduğu için, bunu biliyoruz , bu nedenle bu bunu basitleştirirPθ(XT|XT1,...,X1)=Pθ(XT|XT1)

Pθ(X1,...,XT)=Pθ(XT|XT1)Pθ(X1,...,XT1)

Şimdi aynı mantığı kez tekrarlarsanız ,T

Pθ(X1,...,XT)=i=1TPθ(Xi|Xi1)

burada , sürecin başlangıç ​​durumu olarak yorumlanacaktır. Sağ taraftaki terimler sadece geçiş matrisinin öğeleridir. İstediğiniz günlük olasılığı olduğundan, son yanıt:X0

L(θ)=i=1Tlog(Pθ(Xi|Xi1))

Bu tek bir markov zincirinin olasılığıdır - veri kümeniz birkaç (bağımsız) markov zinciri içeriyorsa, tam olasılık bu formun terimlerinin toplamı olacaktır.


Vay be, cevap için çok teşekkürler. Bu durumda MLE'den alınan "geçiş" olasılığı, değil mi? Pθ
fsociety

@ ph_singer, çok hoş geldiniz. , parametre değeri göz önüne alındığında durumundan durumuna geçme olasılığıdır . Geçiş matrisine hiçbir yapı koymadıysanız (ki bu kulağa nasıl geliyor) o zaman sadece geçiş olasılıklarının vektörünü gösterir (ve MLE'ler soru ifadenizde doğru bir şekilde belirttiğiniz gibi sadece örnek oranlardır), yani, evet : sadece durumundan hamle örnek oranı olacak durum içinde sona erdi . Pθ(Xi|Xi1)Xi1XiθθPθ^MLE(Xi|Xi1)Xi1Xi
Makro

Tekrar teşekkürler! Sadece bir soru daha: Başka bir sipariş kullanırsam (örneğin, k = 2), bu süreç nasıl çalışır?
fsociety

"Sipariş" ile ne demek istediğinizi açıklığa kavuşturabilir misiniz?
Makro

(+1) OP muhtemelen ikinci dereceden bir MC'yi belirtmek için anlamına gelir , yani önceki iki duruma bağlı olarak sadece en sonuncusu olan . k=2Xi1,Xi2Xi1
kardinal
Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.