从马尔可夫链对象导出转移矩阵

Exporting the Transition Matrix from a markovchain object

我将一系列状态组织成一个数据框,如下所示:

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

我使用了 markovchain 的 markovchainListfit 函数,因为我想为这个事件序列构造 t运行 位置矩阵,所以我 运行:

myModel <- markovchainListFit(data = myData) 

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

代码似乎有效,但我找回了一个马尔可夫链对象,我无法仅将 t运行sition 矩阵导出到数据框或矩阵,以便我可以分析它。我知道这非常简单,但出于某种原因我没有看到它。

我试过以下方法:

  1. 将整个列表导出到一个 txt 文件,但没有成功:

    capture.output(summary(myModel), file = "myModel.txt")
    
  2. 像这样访问列表的元素,但我得到一个空列表:

    tranMatrix <- myModel$estimate@markovchains 
    
  3. 我还尝试强制转换为我知道如何使用的不同对象,例如数据框,使用:

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

    但是我得到一个错误

    (Error in as(myModel, from = "markovchain", to = "data.frame"): unused arguments (from = "markovchain", to = "data.frame")

谢谢!

为了提取转换矩阵,我使用了:

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