Cálculo da probabilidade logarítmica para determinado MLE (cadeias de Markov)

9

Atualmente, estou trabalhando com cadeias de Markov e calculei a Estimativa de Máxima Verossimilhança usando probabilidades de transição, conforme sugerido por várias fontes (ou seja, número de transições de a para b dividido pelo número de transições gerais de a para outros nós).

Agora eu quero calcular a probabilidade de log do MLE.

fsociety
fonte
Você já calculou a estimativa de probabilidade máxima das probabilidades de transição e agora deseja calcular a probabilidade de log do que exatamente?
Nick

Respostas:

11

Seja um caminho da cadeia de markov e seja a probabilidade de observar o caminho quando é o verdadeiro valor do parâmetro (também conhecido como a função de probabilidade para ). Usando a definição de probabilidade condicional, sabemos P θ ( X 1 , . . . , X T ) θ θ{Xi}i=1TPθ(X1,...,XT)θθ

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

Como essa é uma cadeia de markov, sabemos que , então isso simplifica isso paraPθ(XT|XT1,...,X1)=Pθ(XT|XT1)

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

Agora, se você repetir essa mesma lógica vezes, obtémT

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

onde deve ser interpretado como o estado inicial do processo. Os termos do lado direito são apenas elementos da matriz de transição. Como foi a probabilidade de log que você solicitou, a resposta final é:X0

eu(θ)=Eu=1 1Tregistro(Pθ(XEu|XEu-1 1))

Essa é a probabilidade de uma única cadeia de markov - se o seu conjunto de dados incluir várias cadeias de markov (independentes), a probabilidade total será uma soma dos termos deste formulário.

Macro
fonte
Uau, muito obrigado pela resposta. Nesse caso, é a probabilidade de "transição" obtida do MLE, certo? Pθ
Fsociety
@ ph_singer, você é muito bem-vindo. é a probabilidade de passar do estado para , considerando o valor do parâmetro . Se você não impôs nenhuma estrutura à matriz de transição (que é o que parece), então apenas indica o vetor de probabilidades de transição (e os MLEs são apenas as proporções da amostra, conforme indicado corretamente na declaração da pergunta), então sim : será apenas a proporção da amostra de movimentos do estado que terminaram no estado . X i - 1 X i θ θ P θ M G E ( X i | X I - 1 ) X i - 1 X iPθ(XEu|XEu-1 1)XEu-1 1XEuθθPθ^MeuE(XEu|XEu-1 1)XEu-1 1XEu
Macro
Obrigado novamente! Apenas mais uma pergunta: se eu usar outro pedido (por exemplo, k = 2), como esse processo funcionaria?
fsociety
Você pode esclarecer o que você quer dizer com "ordem"?
Macro
(+1) O OP provavelmente significa para denotar um MC de segunda ordem , ou seja, dependendo dos dois estados anteriores e não apenas o mais recente . X i - 1 , X i - 2 X i - 1k=2XEu-1 1,XEu-2XEu-1 1
cardeal