The `simufreqD`

function simulate single population allele frequencies for independent loci.
Allele frequencies are generated as random deviates from a Dirichlet distribution, whose parameters control
the mean and the variance of the simulated allele frequencies.

1 | ```
simufreqD(nloc = 1, nal = 2, alpha = 1)
``` |

`nloc` |
the number of loci to simulate |

`nal` |
the numbers of alleles per locus. Either an integer, if the loci have the same number of alleles, or an integer vector, if the number of alleles differ between loci |

`alpha` |
the parameter used to simulate allele frequencies from the Dirichlet distribution. If the
When the number of alleles differ between loci, |

Allele frequencies for independent loci are simulated using a Dirichlet distribution with parameter
`alpha`

. At a given locus L with n alleles, the allele frequencies are modeled as a vector of random
variables
p=(p1, ..., pn), following a Dirichlet distribution with parameters:

alpha = (alpha1, ..., alphan) where p1+...+pn=1 and alpha1,..., alphan > 0.

A matrix containing the simulated allele frequencies. The data is presented in the format of the Journal of Forensic Sciences for genetic data: allele names are given in the first column, and frequencies for a given allele are read in rows for the different markers in columns. When an allele is not observed for a given locus, the value is coded NA (instead of "-" in the original format).

The code used here for the generation of random Dirichlet deviates was previously implemented in the gtools library.

Hinda Haned contact@hindahaned.info

Johnson NL, Kotz S, Balakrishnan N. Continuous Univariate Distributions, vol 2. John Wiley & Sons, 1995.

Wright S. The genetical structure of populations. Ann Eugen 1951;15:323-354.

1 2 3 |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.