Я пытаюсь смоделировать траекторию выборки, используя непрерывную цепочку времени Маркова.Это то, что я сделал:
set.seed(183427)
require(ECctmc)
# rates
r1 <- 1 # 1->2
r2 <- 0.75 # 2->3
r3 <- 0.5 # 3->1
r4 <- 0.5 # 3-> 2
Q <- matrix(c(-r1, r1, 0, 0, -r2, r2, r3, r4, -(r3+r4)),nrow = 3, byrow = TRUE)
# sample path
path <- sample_path(a=1, b = 3,t0=0, t1=5, Q=Q)
plot(stepfun(x=path[1:(nrow(path)-1),"time"], y = path[,"state"]), xlim = c(0,5), xlab = "Time", ylab = "State", main = "Sample path")
Это отлично работает, и я получаю график в качестве вывода, где у меня есть состояние по оси Y и время по оси X.Однако в этом случае я должен указать свое конечное состояние b.Я хочу изменить этот код так, чтобы мне не нужно было указывать конечное состояние и просто указывать в качестве входных данных матрицу скоростей, начальное распределение и длительность временного интервала T.Но я не знаю, как это сделать.Помощь будет оценена.