Марковская цепь подходит из множества отдельных цепочек (в R)

Question

С использованием markovchain пакет, я работаю с набором данных, состоящим из шести ежемесячных наблюдений для каждого из 23 тыс. человек.

Когда я иду, чтобы соответствовать DTMC, используя markovchainFit функция, функция, кажется, хочет принять то, что было бы только одним из моих 23k людей. Как я могу установить DTMC на полную совокупность последовательностей по 23 к 6 периодов?

(К вашему сведению - я в состоянии рассчитать / построить / описать / и т. Д. MC просто отлично - я просто хочу иметь возможность генерировать некоторые точные прогнозы и использовать преимущества остальной части функциональности пакета, и для этого мне кажется, что мне нужно приспособленный объект MC?)

Итак: как я могу соответствовать объекту MC, используя markovchain или другой пакет, использующий кучу наблюдений той же последовательности из 6 периодов, и который позволит мне затем сгенерировать некоторые прогнозы для будущих шагов?

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