马尔可夫链适合许多单独的链(在R中)



使用markovchain软件包,我正在使用一个数据集,该数据集由23K个人中的每个人每月观察组成。

当我使用markovchainFit函数拟合DTMC时,该函数似乎想接收我的23K个人之一。我如何在23k 6 period序列的全部人群上安装DTMC?

(fyi-我能够计算/绘图/描述/等。MC很好 - 我只想能够生成一些干净的预测并利用包装的其余功能,为此似乎我需要一个安装的MC对象?(

so:如何使用markovchain或其他软件包使用一堆观察到相同的6个周期序列拟合MC对象,然后让我为将来的步骤生成一些预测?

markovchainFit功能可以处理matrixdata.framelist作为输入:

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

最新更新