# NOT RUN {
# Construct a population with additive and epistatic effects
pop <- Population(
popSize = 200, map = map100snp, QTL = 20,
alleleFrequencies = runif(100),
broadH2 = 0.9, narrowh2 = 0.6, traitVar = 40
)
pop <- addEffects(pop)
pop <- attachEpiNet(pop)
# Run the simulator
pop2 <- runSim(pop, generations = 150)
# Retrieve the population pedigree from the simulation
ped <- getPedigree(pop2)
# Re-run the simulation using the same pedigree
pop3 <- runSim(pop, ped)
# }
Run the code above in your browser using DataLab