我正在尝试了解为带有r包“ msm”的面板数据确定马尔可夫链的q矩阵的方法。
（我是msm模型的新手）
在r-package示例中，确定为：
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项为零的非对角位置上包含零。所有其他非对角位置均包含对应的跃迁强度的初始值”
但是我不明白为什么选择这个值。
此外，我有一个经典的面板数据，该数据由与状态相关的整个时间的ID决定。看起来像这样：
（这是示例，不是我的真实数据。假设有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。
感谢您的帮助