从马尔可夫链对象导出转移矩阵
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 矩阵导出到数据框或矩阵,以便我可以分析它。我知道这非常简单,但出于某种原因我没有看到它。
我试过以下方法:
将整个列表导出到一个 txt 文件,但没有成功:
capture.output(summary(myModel), file = "myModel.txt")
像这样访问列表的元素,但我得到一个空列表:
tranMatrix <- myModel$estimate@markovchains
我还尝试强制转换为我知道如何使用的不同对象,例如数据框,使用:
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
我将一系列状态组织成一个数据框,如下所示:
Year1 Year2 Year3 ...
1 2 5 ...
3 9 4 ...
我使用了 markovchain 的 markovchainListfit
函数,因为我想为这个事件序列构造 t运行 位置矩阵,所以我 运行:
myModel <- markovchainListFit(data = myData)
其中 myData 是上面显示的数据。
代码似乎有效,但我找回了一个马尔可夫链对象,我无法仅将 t运行sition 矩阵导出到数据框或矩阵,以便我可以分析它。我知道这非常简单,但出于某种原因我没有看到它。
我试过以下方法:
将整个列表导出到一个 txt 文件,但没有成功:
capture.output(summary(myModel), file = "myModel.txt")
像这样访问列表的元素,但我得到一个空列表:
tranMatrix <- myModel$estimate@markovchains
我还尝试强制转换为我知道如何使用的不同对象,例如数据框,使用:
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