Simulates nucleotide data based on parameters under the SELAC model
SelonSimulator(phy, pars, nuc.optim_array, nuc.model, diploid = TRUE,
start.vals_array = NULL)The phylogenetic tree with branch lengths.
A vector of parameters used for the simulation. They are ordered as follows: a0, a1, a2, Ne, base.freqs for A C G, and the nucleotide rates.
A vector of optimal nucleotide for each site to be simulated.
Indicates what type nucleotide model to use. There are three options: "JC", "GTR", or "UNREST".
A logical indicating whether or not the organism is diploid or not.
A vector of nucleotides to be used as the starting nucleotide for each site in the simulation.
Simulates a nucleotide matrix using parameters under the SELON model. Note that the output can be written to a fasta file using the write.dna() function in the ape package.