estimate_pt {cfda} | R Documentation |
Estimate probabilities to be in each state
estimate_pt(data, NAafterTmax = FALSE)
data |
data.frame containing |
NAafterTmax |
if TRUE, return NA if t > Tmax otherwise return the state associated with Tmax (useful when individuals has different lengths) |
A list of two elements:
t: vector of time
pt: a matrix with K (= number of states) rows and with length(t)
columns containing the probabilities to be in each state at each time.
Cristian Preda, 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) d_JK2 <- cut_data(d_JK, 10) # estimate probabilities estimate_pt(d_JK2)