View source: R/create_haplotypes.R
create_haplotypes | R Documentation |
Uses one of multiple methods to create variant haplotypes from a reference genome.
See haps_functions
for the methods available.
create_haplotypes(
reference,
haps_info,
sub = NULL,
ins = NULL,
del = NULL,
epsilon = 0.03,
n_threads = 1,
show_progress = FALSE
)
reference |
A |
haps_info |
Output from one of the |
sub |
Output from one of the |
ins |
Output from the |
del |
Output from the |
epsilon |
Error control parameter for the "tau-leaping" approximation to
the Doob–Gillespie algorithm, as used for the indel portion of the simulations.
Smaller values result in a closer approximation.
Larger values are less exact but faster.
Values must be |
n_threads |
Number of threads to use for parallel processing.
This argument is ignored if OpenMP is not enabled.
Threads are spread across chromosomes, so it
doesn't make sense to supply more threads than chromosomes in the reference genome.
Defaults to |
show_progress |
Boolean for whether to show a progress bar during processing.
Defaults to |
A haplotypes
object.
Cao, Y., D. T. Gillespie, and L. R. Petzold. 2006. Efficient step size selection for the tau-leaping simulation method. The Journal of Chemical Physics 124(4): 044109.
Doob, J. L. 1942. Topics in the theory of markoff chains. Transactions of the American Mathematical Society 52(1): 37–64.
Gillespie, D. T. 1976. A general method for numerically simulating the stochastic time evolution of coupled chemical reactions. Journal of Computational Physics 22(4): 403–434.
Wieder, N., R. H. Fink, and F. von Wegner. 2011. Exact and approximate stochastic simulation of intracellular calcium dynamics. Journal of Biomedicine and Biotechnology 2011: 572492.
r <- create_genome(10, 1000)
v_phylo <- create_haplotypes(r, haps_phylo(ape::rcoal(5)), sub_JC69(0.1))
v_theta <- create_haplotypes(r, haps_theta(0.001, 5), sub_K80(0.1, 0.2))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.