Я пытаюсь понять, как определить матрицу q для цепочки Маркова с данными панели с пакетом r "msm".
(я новичок в моделях MSM)
В примере r-пакета оно определяется как:
twoway4.q <- rbind(c(0, 0.25, 0, 0.25), c(0.166, 0, 0.166, 0.166),
c(0, 0.25, 0, 0.25), c(0, 0, 0, 0))
пример ссылки: https://www.jstatsoft.org/article/view/v038i08/v38i08.pdf
Но я не понимаю, как определить эти значения (0, 0,25, 0, 0,25). По словам автора, «содержащего нули в недиагональных позициях, где элементы Q равны нулю. Все остальные недиагональные позиции содержат начальное значение для соответствующей интенсивности перехода»
Но я не могу понять, почему выбирают эти значения.
Кроме того, у меня есть классическая панель данных, определяемая идентификатором по времени, связанному с состоянием. Выглядит так:
(Это пример, а не мои реальные данные. Предположим, 4 состояния)
id year month state
1 1960 1 2
1 1960 2 1
1 1960 3 4
1 1960 4 3
. . . .
. . . .
. . . .
1 2016 12 4
2 1960 1 3
. . . .
. . . .
. . . .
n 2016 12 state_n
Мне нужно создать цепочку Маркова с данными, похожими на пример выше, но я не знаю, как определить qmatrix.
Спасибо за вашу помощь