hierfstat (version 0.5-11)

sim.genot: 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

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
# }