sim.genot: Simulates genotypes in an island model at equilibrium

In jgx65/hierfstat: Estimation and Tests of Hierarchical F-Statistics

Simulates genotypes in an island model at equilibrium

Description

Simulates genotypes from several individuals in several populations at several loci in an island model at equilibrium. The islands may differ in size and inbreeding coeeficients.

Usage

sim.genot(size=50,nbal=4,nbloc=5,nbpop=3,N=1000,mig=0.001,mut=0.0001,f=0)

Arguments

size	The number of individuals to sample per population
nbal	The maximum number of alleles present at a locus
nbloc	The number of loci to simulate
nbpop	The number of populations to simulate
N	The population sizes for each island
mig	the proportion of migration among islands
mut	The loci mutation rate
f	the inbreeding coefficient for each island

Value

a data frame with nbpop*size lines and nbloc+1 columns. Individuals are in rows and genotypes in columns, the first column being the population identifier

Author(s)

Jerome Goudet [e-mail]

Examples

## Not run: 
dat<-sim.genot(nbpop=4,nbal=20,nbloc=10,mig=0.001,mut=0.0001,N=c(100,100,1000,1000),f=0)
betas(dat)$betaiovl

## End(Not run)

jgx65/hierfstat documentation built on April 20, 2023, 8:34 a.m.