simSeq from the phangorn package.Simulate sequences for a given evolutionary tree, using a standard model
of sequence evolution along the branches. Code for this function was heavily
inspired by the function simSeq from the phangorn package.
sim_normal(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
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
Klaus Schliep klaus.schliep@gmail.com