cut_data {cfda} | R Documentation |
Cut data to a maximal given time
cut_data(data, Tmax)
data |
data.frame containing |
Tmax |
max time considered |
a data.frame with the same format as data
where each individual has Tmax
as last time entry.
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 = 10, K = K, P = PJK, lambda = lambda_PJK, Tmax = 10) tail(d_JK) # cut at Tmax = 8 d_JK2 <- cut_data(d_JK, Tmax = 8) tail(d_JK2)