estimate_Markov {cfda} | R Documentation |
Calculates crude initial values for transition intensities by assuming that the data represent the exact transition times of the Markov process.
estimate_Markov(data)
data |
data.frame containing |
list of two elements: Q
, the estimated transition matrix, and lambda
, the estimated time spent in each state
Cristian Preda
# 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 = 100, K = K, P = PJK, lambda = lambda_PJK, Tmax = 10) # estimation mark <- estimate_Markov(d_JK) mark$P mark$lambda