Dominance deviations

Returns dominance deviations for all traits

Usage

dd(pop, simParam = NULL)

Arguments

pop

an object of Pop-class

simParam

an object of SimParam

Examples

#Create founder haplotypes
founderPop = quickHaplo(nInd=10, nChr=1, segSites=10)

#Set simulation parameters
SP = SimParam$new(founderPop)
SP$addTraitAD(10, meanDD=0.5)
SP$setVarE(h2=0.5)
SP$nThreads = 1L

#Create population
pop = newPop(founderPop, simParam=SP)
dd(pop, simParam=SP)
#>            Trait1
#>  [1,]  0.47307923
#>  [2,]  0.07417245
#>  [3,] -0.17579194
#>  [4,]  0.08077610
#>  [5,] -0.27650589
#>  [6,] -0.17123583
#>  [7,]  0.42400041
#>  [8,] -0.41609159
#>  [9,]  0.06537566
#> [10,] -0.07777860