sim_pop: Simulate a population

In neyhartj/pbsim: Simulations of Plant Breeding Experiments

Simulate a population

Description

Simulates a random population of given population size. SNPs are simulated to be in linkage equilibrium.

Usage

sim_pop(genome, n.ind, ignore.gen.model = FALSE)

Arguments

genome	A genome object.
n.ind	The number of individual in the population.
ignore.gen.model	Logical - should the gene model be ignored?

Value

A pop object.

Examples

# Simulate a genome
n.mar  <- c(505, 505, 505)
len <- c(120, 130, 140)

genome <- sim_genome(len, n.mar)

# Simulate a quantitative trait influenced by 50 QTL
qtl.model <- matrix(NA, 50, 4)
genome <- sim_gen_model(genome = genome, qtl.model = qtl.model, 
                        add.dist = "geometric", max.qtl = 50)
                        
# Simulate the population
pop <- sim_pop(genome = genome, n.ind = 100)

neyhartj/pbsim documentation built on Nov. 11, 2023, 4:07 p.m.