simulate a sequence assuming substitutions are only accumulated along the branches, using the explicit simulation method (e.g. reverse substitutions are modeled explicitly)
sim_normal_explicit(x, l = 1000, Q = NULL, bf = NULL, rootseq = NULL, rate = 1)list with four items
alignment Phydat object with the resulting alignment
rootseq the rootsequence used
total_branch_substitutions total number of substitutions accumulated on the branches
total_node_substitutions total number of substitutions accumulated at the nodes
a phylogenetic tree tree, i.e. an object of class
phylo or and object of class pml.
length of the sequence to simulate.
the rate matrix.
base frequencies.
a vector of length l containing the root sequence, other root sequence is randomly generated.
mutation rate or scaler for the edge length, a numerical value greater than zero.