# NOT RUN {
# Construct a new population with additive effects
pop <- Population(
popSize = 20, map = map100snp, QTL = 20,
broadH2 = 0.4, narrowh2 = 0.4, traitVar = 40,
alleleFrequencies = runif(100, 0.05, 0.5)
)
pop <- addEffects(pop)
# Find the additive contribution to the individuals' phenotypes
hap <- getHaplo(pop)
hap <- (hap[[1]] + hap[[2]])[, getQTL(pop)$Index]
(hap %*% getAddCoefs(pop))[, 1] + getAddOffset(pop)
# Compare with additive component from getComponents()
getComponents(pop)$Additive
# }
Run the code above in your browser using DataLab