Converts a dosage matrix into a matrix of centered dominance deviations. This sets up the "classical" (aka "Statistical") partition of additive dominance in terms of breeding values and dom. deviations. See function dose2domDevGenotypic() to set-up the "genotypic" (aka "biological") partition in terms of genotypic effects. See Vitezica et al. 2013. Also see Varona et al. 2018. https://doi.org/10.3389/fgene.2018.00078
dose2domDev(M)
| M | dosage matrix. Assumes SNPs in M coded 0, 1, 2 (requires rounding dosages to integers). M is Nind x Mrow, numeric matrix, with row/columanes to indicate SNP/ind ID |
|---|
domDev<-dose2domDev(M)#> Error in is.data.frame(x): object 'M' not found