将转换矩阵拟合到r中的交叉表中

Question

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

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

  1 2 3 4 5 6 
     1 2639 156 69 49 21 42 
     2 150 300 53 11 1 1 
     3 74 20 363 48 8 4 
     4 60 7 46 384 65 12 
     5 42 2 5 75 275 63 
     6 40 0 1 15 52 402 

library(markovchain) 
mcobject<-as(mytableX, "markovchain") 


     1   2   3   4   5   6 
1 0.88676075 0.052419355 0.023185484 0.01646505 0.007056452 0.014112903 
2 0.29069767 0.581395349 0.102713178 0.02131783 0.001937984 0.001937984 
3 0.14313346 0.038684720 0.702127660 0.09284333 0.015473888 0.007736944 
4 0.10452962 0.012195122 0.080139373 0.66898955 0.113240418 0.020905923 
5 0.09090909 0.004329004 0.010822511 0.16233766 0.595238095 0.136363636 
6 0.07843137 0.000000000 0.001960784 0.02941176 0.101960784 0.788235294 

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([email protected], 
           countsTransitionMatrix=ciau.matrix, confidencelevel=.95)