statetable {cfda} | R Documentation |
Calculates a frequency table counting the number of times each pair of states were observed in successive observation times.
statetable(data, removeDiagonal = FALSE)
data |
data.frame containing |
removeDiagonal |
if TRUE, does not count transition from a state i to i |
a vector of length K
containing the total time spent in each state
Quentin Grimonprez
# Simulate the Jukes-Cantor model of nucleotide replacement K <- 4 PJK <- matrix(1/3, nrow = K, ncol = K) - diag(rep(1/3, K)) lambda_PJK <- c(1, 1, 1, 1) d_JK <- generate_Markov(n = 10, K = K, P = PJK, lambda = lambda_PJK, Tmax = 10) # table of transitions statetable(d_JK)