examples/sim_markov_forward.R

end_points <- list(inc_times = c(0,10),
                   inc_behav = c(3))
param <- list(lambda = c(0.1,0.2,0.3),
              q = matrix(c(NA,0.5,0.3,0.35,NA,0.7,0.65,0.5,NA),nrow=3))

set.seed(123)
sim_markov_forward(3, end_points, param)
# $times
# [1]  0.000000  2.811524 10.000000
#
# $states
# [1] 3 1 1

set.seed(2712)
sim_markov_forward(3, end_points, param)
# $times
# [1]  0.000000  1.898085  3.603405  4.275325 10.000000
#
# $states
# [1] 3 2 1 3 3
a-parton/CTStepTurn documentation built on May 14, 2019, 4:17 p.m.