gl.sim.Neconst: Simulate a population with constant mutation rate

In dartR.sim: Computer Simulations of 'SNP' Data

Simulate a population with constant mutation rate

Description

This function simulates a population with a constant mutation rate using a beta distribution for allele frequencies.

Usage

gl.sim.Neconst(ninds, nlocs, mutation_rate = 1e-08, verbose = 0)

Arguments

ninds	Number of individuals in the population [required].
nlocs	Number of loci in the population [required].
mutation_rate	Mutation rate per generation (default is 1e-8) [default 1e-8].
verbose	Verbosity level (default is 0).

Details

The function generates a genlight object with the specified number of individuals and loci, simulating allele frequencies based on a beta distribution. The mutation rate is used to calculate theta, which in turn is the parameter in the beta function.

Value

A genlight object representing the simulated population.

Author(s)

Bernd Gruber (Post to https://groups.google.com/d/forum/dartr)

Examples

# Simulate a population with 50 individuals and 4000 loci
gg <- gl.sim.Neconst(ninds = 50, nlocs = 4000, mutation_rate = 1e-8, verbose = 0)
dartR.popgen::gl.sfs(gg)

dartR.sim documentation built on March 17, 2026, 5:06 p.m.