使用markovchain
软件包,我正在使用一个数据集,该数据集由23K个人中的每个人每月观察组成。
当我使用markovchainFit
函数拟合DTMC时,该函数似乎想接收我的23K个人之一。我如何在23k 6 period序列的全部人群上安装DTMC?
(fyi-我能够计算/绘图/描述/等。MC很好 - 我只想能够生成一些干净的预测并利用包装的其余功能,为此似乎我需要一个安装的MC对象?(
so:如何使用markovchain
或其他软件包使用一堆观察到相同的6个周期序列拟合MC对象,然后让我为将来的步骤生成一些预测?
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