用过渡矩阵模拟马尔科夫链
#transition matrix
P = matrix(c(.4, .3, .5, 0, .5, 0, 0, .3, .4), nrow = 3, ncol = 5,
byrow = TRUE)
我不清楚你到底在找什么。如果你是想获得以下信息后的过渡矩阵 10
和 100
次,你可以跟踪进化进度,通过。Reduce
,其中选项 accumulate
应设置为 TRUE
P10 <- Reduce(`%*%`,replicate(10,P,simplify = FALSE),accumulate = TRUE)
P100 <- Reduce(`%*%`,replicate(100,P,simplify = FALSE),accumulate = TRUE)