SDselectiongain | R Documentation |
This function is used to calculate the standard deviation of sel gain acording to longin 2015
SDselectiongain(Ob, maseff, VGCAandE, VSCA, VLine, years, Genotypes)
Ob |
matrix object produced by the function multistageoptimum.search or multistageoptiumum.grid |
maseff |
is the efficiency of marker-assisted selection (MAS). The default value is NA, which means there is no MAS. If a value between 0 and 1 is assigned to |
VGCAandE |
is the vector of variance components of genetic effect, genotype \times location interaction, genotype \times year interaction, genotype \times location \times year interaction and the plot error. When |
VSCA |
is the vector of variance components for specific combining ability. The default value is 0,0,0,0. |
VLine |
is the vector of variance components for line per se. The default value is 0,0,0,0,0. |
years |
Duration of the breeding scheme in years, it is used only to compute the anual selection gain |
Genotypes |
character vector to indicate the function which variance components we are using. Pssible values are "Hybrids" if we are using GCA and SCA variance components or "Lines" if we are using line perse variance components |
for the new added to parameters "alpha.nursery" and "cost.nursery" since v2.0.47:
After producing new DH lines, breeders do NOT go directly for a selection stage in the field, neither for genomic selection. Most of the times, they prefer to make a small field experiment (called "nursery") in which all DH lines are observed and discarded for other traits as disease resistance. That means, all DH lines with poor resistance will be discarded. At the end of the nursery stage only certain amount of DH lines (alpha) advance to the first selection stage (phenotypic or genomic). Specially in maize that makes sense, because in experience around 90 percent of the new DH lines are very weak in terms of per se performance what make them not suitable as new hybrid parents. Then, budget should not be used to make genotyping on or testcrossing with them. Only the alpha fraction should be used for entering the stage 1 of the multistageoptimum.search function.
More details are available in the Crop Science and Computational Statistics papers.
The output is equivalent to the matrix object produced by the functions multistageoptimum.search or multistageoptimum.grid but with two columns added, one for the values of the anual selection gain and the second for the standard deviation of selection gain
no further comment
Jose Marulanda
C. Longin, X. Mi and T. Wuerschum. Genomic selection in wheat: optimum allocation of test resources and comparison of breeding strategies for line and hybrid breeding. Theoretical and Applied Genetics 128: 1297-1306. 2015.
C. Longin, H.F. Utz., J. Reif, T. Wegenast, W. Schipprack and A.E. Melchinger. Hybrid maize breeding with doubled haploids: III. Efficiency of early testing prior to doubled haploid production in two-stage selection for testcross performance. Theor. Appl. Genet. 115: 519-527, 2007.
E.L. Heffner, A.J. Lorenz, J.L. Jannink, and M.E. Sorrells. Plant breeding with genomic selection: gain per unit time and cost. Crop Sci. 50: 1681-1690, 2010.
selectiongain()
CostProd =c(0.5,1,1) CostTest = c(0.5,1,1) Budget=1021 # Budget is very small here to save time in package checking # for the example in Heffner's paper, please change it to Budget=10021 VCGCAandError=c(0.4,0.2,0.2,0.4,2) VCSCA=c(0.2,0.1,0.1,0.2) Nf=10 maseff=0.4 years=7 # this breeding scheme takes 7 years from the initial cross to the final field testing. # See references for more details Ob<-multistageoptimum.search (maseff=maseff, VGCAandE=VCGCAandError, VSCA=VCSCA, CostProd = CostProd, CostTest = CostTest, Nf = Nf, Budget = Budget, N2grid = c(11, 1211, 30), N3grid = c(11, 211, 5), L2grid=c(1,1,1), L3grid=c(6,6,1), T2grid=c(1,2,1), T3grid=c(3,5,1), R2=1, R3=1, alg = Miwa(), detail=TRUE, fig=FALSE, t2free=TRUE) SDselectiongain(Ob=Ob,maseff=maseff,VGCAandE=VCGCAandError,VSCA=VCSCA, years=years,Genotypes="Hybrids")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.