У меня проблемы с реализацией функции rmarkovchain из пакета markovchain. Я хотел бы смоделировать случайную последовательность с помощью функции rmarkovchain ie:
library(markovchain)
mcWeather <- new("markovchain", states = c("sunny", "cloudy", "rain"),
transitionMatrix = matrix(data = c(0.70, 0.2, 0.1,
0.3, 0.4, 0.3,
0, 0, 1), byrow = byRow, nrow = 3),
name = "Weather")
rmarkovchain(n=100, object = mcWeather, t0 = "sunny")
Однако одним из указанных состояний является состояние поглощения (дождь). Поэтому вместо определения n = "some number" я хотел бы смоделировать столько, сколько n, пока не произойдет состояние поглощения, а затем закончить моделирование.
Спасибо за любые предложения.