mpow(P, 18)
をベクトル形式&行列形式で書こうとしています。誰もそれで私を助けることができますか?遷移確率行列が与えられているマルコフ連鎖の定常分布を得るには
また、私は各状態の定常分布を見つけようとしています。ここで
Pi_0 = ?
Pi_1 = ?
Pi_2 = ?
...
Pi_5 = ?
は、私が書いたコードです:あなたの質問に
P <- matrix(c(0, 0, 0, 0.5, 0, 0.5, 0.1, 0.1, 0, 0.4, 0, 0.4, 0, 0.2, 0.2, 0.3, 0, 0.3, 0, 0, 0.3, 0.5, 0, 0.2, 0, 0, 0, 0.4, 0.6, 0, 0, 0, 0, 0, 0.4, 0.6), nrow = 6, ncol = 6, byrow = TRUE)
mpow <- function(P, n) {
if (n == 0) diag(nrow(P))
else if (n == 1) P
else P %*% mpow(P, n - 1)
}
mpow(P, 18)
固定分布は、それぞれの状態で終わる確率を表す行ベクトルです。行の合計を 'rowSums(mp(P、18))'でチェックすると、それらは合計で1になります。 – toni057
あなたの定常的な分布です。 – toni057
私はそれが正しい答えだとは思わない。 Pi_0 ..... Pi_5は、おそらく1未満の小数部であるべきです。 – PeterNiklas