从markovchain对象导出过渡矩阵

Question

我将状态序列组织成一个数据框，如下所示：

Year1 Year2 Year3 ...
1       2     5   ...
3       9     4   ...

我使用markovchain的markovchainListfit函数是因为我想为此事件序列构造转换矩阵，所以我运行了：

myModel <- markovchainListFit(data = myData) 

其中 myData 是上面显示的数据。

代码似乎可以正常工作，但是我得到了markovchain对象，但在将过渡矩阵仅导出到数据帧或矩阵时遇到了麻烦，因此我可以对其进行分析。我知道这很简单，但是由于某种原因我没有看到它。

我尝试了以下方法：

1. 将整个列表导出到无法正常工作的txt文件：

   capture.output(summary(myModel), file = "myModel.txt")
   

2. 像这样访问列表的元素，但是我得到了一个空列表：

   tranMatrix <- myModel$estimate@markovchains 
   

3. 我还尝试过强制转换为一个我知道如何使用的对象，例如使用以下数据框：

   as(myModel, from = "markovchain", to = "data.frame") 
   

   但是我遇到了错误

     

   （as（myModel，from =“ markovchain”，to =“ data.frame”中的错误：未使用的参数（from =“ markovchain”，to =“ data.frame”）

谢谢！

Answer 1

我使用以下方法提取过渡矩阵：

myModel$estimate@markovchains[[1]]@transitionMatrix