Simulate allele frequencies

Simulate allele frequencies

Usage

simulate_allele_frequencies(alpha, num_loci)

Arguments

alpha

vector parameter controlling the Dirichlet distribution

num_loci

total number of loci to draw

Value

Numeric matrix with length(alpha) rows and num_loci columns; each column is an allele frequency vector for one locus.

Details

Simulate allele frequency vectors as a draw from a Dirichlet distribution

Examples

simulate_allele_frequencies(alpha = c(1, 1, 1, 1), num_loci = 3)
#>            [,1]      [,2]      [,3]
#> [1,] 0.44764475 0.1051512 0.3560853
#> [2,] 0.05500502 0.4197373 0.3549046
#> [3,] 0.11079711 0.1933853 0.1738189
#> [4,] 0.38655312 0.2817262 0.1151912