Question

使用markovchain套餐，我正在使用一个数据集，该数据集包含每月23个人的六个月观察结果。

当我使用markovchainFit功能去适应DTMC时，该功能似乎想要接收我的23k个人中的一个。如何在完整的23k 6周期序列中安装DTMC？

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

那么：我如何使用markovchain或其他包使用相同的6个周期序列的一堆观察来拟合MC对象，然后让我为未来的步骤生成一些预测？

Answer 1

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

马尔可夫链适合许多单独的链条（在R中）

1 个答案: