Skip to contents

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)

#Create population
pop = newPop(founderPop, simParam=SP)
dd(pop, simParam=SP)
#>             Trait1
#>  [1,] -0.538699475
#>  [2,]  0.009939649
#>  [3,]  0.498418920
#>  [4,]  0.357675877
#>  [5,] -0.330793890
#>  [6,]  0.334681096
#>  [7,] -0.411143211
#>  [8,] -0.225089870
#>  [9,]  0.289978568
#> [10,]  0.015032337