拟合过渡矩阵到r中的可交换性

Question

我在安装转换矩阵时遇到了一些麻烦。 在分析面板数据时，我需要将转换矩阵（具有标准误差）拟合到下面显示的交叉矩阵或基础数据。它显示参与者在前两波中的收入状态之间转换。

通过使用markovchain包将其转换为markovchain对象，我只获得转换百分比，但没有标准错误。我怎么得到那些？

$s

Answer 1

您可以使用 multinomialConfidenceIntervals 功能。

library(markovchain)
#get your data
ciau.matrix<-matrix(c(2639, 156  ,69,   49  , 21,   42,
150,  300  ,53 ,  11  , 1  ,  1,
74,   20  , 363 , 48   ,8   , 4,
60,   7  ,  46   ,384 , 65   ,12,
42,   2 ,   5    ,75 ,  275 , 63,
40,   0,    1    ,15,   52  , 402),nrow = 6,byrow = TRUE)

ciau.table<-as.table(ciau.matrix)

ciau.Markovchain<-as(ciau.table,"markovchain")
#use the multinomialConfidenceInterval to obtain confidence intervals
multinomialConfidenceIntervals(transitionMatrix=ciau.Markovchain@transitionMatrix, 
                               countsTransitionMatrix=ciau.matrix, confidencelevel=.95)