0

使用该markovchain软件包，我正在处理一个数据集，该数据集由 23k 个人每个月的六个月观察组成。

当我使用该功能安装 DTMC 时markovchainFit，该功能似乎想要容纳我的 23k 个人中的一个。如何在 23k 6 周期序列的全部人口中拟合 DTMC？

（仅供参考 - 我能够计算/绘图/描述/等。MC就好了 - 我只想能够生成一些干净的预测并利用包的其余功能，为此我似乎需要一个合适的 MC 对象？）

那么：如何使用markovchain相同 6 周期序列的一组观察值来拟合 MC 对象或另一个包，这将让我对未来的步骤产生一些预测？

4

1 回答 1

0

该markovchainFit函数可以处理matrix,data.frame和list作为输入：

library(markovchain)
#getting from matrix / data.frame
data(holson) #load a data.frame (also it works with matrices)
head(holson) #load a matrix of pop * time observations
singleMc<-markovchainFit(data=holson[,2:7],name="holson") #fit the MC

#getting from list
myList<-list()
for (i in 1:100) {
  myList[[i]]<-sample(x = c("a","b"),size = 6,replace = TRUE,prob = c(0.5,0.5))
}
singleMcFromList<-markovchainFit(data=myList,name="holson") #fit the MC
于 2017-09-11T14:46:05.530 回答