есть ли способ генерировать случайную выборку из цепочки Маркова более высокого порядка?Я использовал пакет clickstream для оценки цепочки Маркова 2-го порядка, и сейчас я пытаюсь сгенерировать образец из него.Я понимаю, как сделать это из матрицы перехода с функцией randomClickstreams, но это сработало бы только для цепочки Маркова 1-го порядка.
Вот воспроизводимый пример, в котором мы генерируем образец из матрицы перехода и затем подгоняемпорядок цепей Маркова на выборке:
trans_mat <- matrix(c(0, 0.2, 0.7, 0, 0.1,
0.2, 0, 0.5, 0, 0.3,
0.1, 0.1, 0.1, 0.7, 0,
0, 0.4, 0.2, 0.1, 0.3,
0, 0 , 0 , 0, 1), nrow = 5)
cls <- randomClickstreams(states = c("P1", "P2", "P3", "P4", "end"),
startProbabilities = c(0.5, 0.5, 0, 0, 0),
transitionMatrix = trans_mat,
meanLength = 20, n = 1000)
# fit 2nd order markov chain:
mc <- fitMarkovChain(clickstreamList = cls, order = 2,
control = list(optimizer = "quadratic"))
Это состоит из 2-х матриц перехода и 2-х лямбда-параметров:
Какмогу ли я использовать эти элементы для создания случайной выборки, скажем, 10000 поездок?