RRBLUP_D {AlphaSimR} | R Documentation |
Fits an RR-BLUP model for genomic predictions that includes dominance effects.
RRBLUP_D(pop, traits = 1, use = "pheno", snpChip = 1, useQtl = FALSE, maxIter = 40L, simParam = NULL, ...)
pop |
a |
traits |
an integer indicating the trait to model, or a function of the traits returning a single value. |
use |
train model using phenotypes "pheno", genetic values "gv", estimated breeding values "ebv", breeding values "bv", or randomly "rand" |
snpChip |
an integer indicating which SNP chip genotype to use |
useQtl |
should QTL genotypes be used instead of a SNP chip. If TRUE, snpChip specifies which trait's QTL to use, and thus these QTL may not match the QTL underlying the phenotype supplied in traits. |
maxIter |
maximum number of iterations. Only used when number of traits is greater than 1. |
simParam |
an object of |
... |
additional arguments if using a function for traits |
#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$addSnpChip(10) #Create population pop = newPop(founderPop, simParam=SP) #Run GS model and set EBV ans = RRBLUP_D(pop, simParam=SP) pop = setEBV(pop, ans, simParam=SP) #Evaluate accuracy cor(gv(pop), ebv(pop))