给定的MLE(马尔可夫链)计算对数似然


9

我目前正在与马尔可夫链一起工作,并使用多个来源建议的转移概率(即,从a到b的转移数除以从a到其他节点的总体转移数)来计算最大似然估计。

我现在要计算MLE的对数似然率。


您已经计算了转移概率的最大似然估计,现在想计算出什么是对数似然?
尼克,

Answers:


11

让是马尔可夫链的路径,并让是观察路径的概率时就是真实参数值(也称为的似然函数)。使用条件概率的定义,我们知道{Xi}i=1TPθ(X1,...,XT)θθ

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

由于这是一个马尔可夫链,我们知道,因此简化为Pθ(XT|XT1,...,X1)=Pθ(XT|XT1)

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

现在,如果您重复相同的逻辑次,您将得到T

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

其中将被解释为过程的初始状态。右侧的术语只是过渡矩阵的元素。由于这是您要求的对数可能性,因此最终答案是:X0

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

这是单个马尔可夫链的可能性-如果您的数据集包含多个(独立的)马尔可夫链,则全部可能性将是这种形式的项的总和。


哇,非常感谢您的回答。在这种情况下,是从MLE中获取的“过渡”概率,对吗?Pθ
fsociety 2013年

@ph_singer,非常欢迎您。是给定参数值从状态到的概率。如果您没有在过渡矩阵上施加任何结构(听起来像这样),那么只是表示过渡概率的向量(MLE只是样本比例,正如您在问题陈述中正确指出的那样),所以,是的: 将只是从状态移动样本比例这结束了在状态。Pθ(Xi|Xi1)Xi1XiθθPθ^MLE(Xi|Xi1)Xi1Xi
2013年

再次感谢!只是一个问题:如果我使用其他顺序(例如k = 2),那么此过程将如何工作?
fsociety

您能否阐明“订单”的含义?
2013年

(+1)OP可能意味着表示一个二阶 MC,即取决于前两个状态而不仅仅是最近的一个。k=2Xi1,Xi2Xi1
红衣主教
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.