Draw random Uniform or Beta ancestral allele frequencies
This is simply a wrapper around \link[stats]{runif} or \link[stats]{rbeta} (depending on parameters) with different defaults and additional validations.
draw_p_anc(m_loci, p_min = 0.01, p_max = 0.5, beta = NA)
m_loci |
Number of loci to draw. |
p_min |
Minimum allele frequency to draw (Uniform case only). |
p_max |
Maximum allele frequency to draw (Uniform case only). |
beta |
Shape parameter for a symmetric Beta.
If |
A length-m vector of random ancestral allele frequencies
# Default is uniform with range between 0.01 and 0.5 p_anc <- draw_p_anc(m_loci = 10) # Use of `beta` triggers a symmetric Beta distribution. # This parameter has increased density for rare minor allele frequencies, # resembling the 1000 Genomes allele frequency distribution p_anc <- draw_p_anc(m_loci = 10, beta = 0.03)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.