SoyNAM (version 1.0)

BLUP: Best Linear Unbias Predictor

Description

Genetic values for a given trait computed by REML.

Usage

BLUP(trait="yield",family="all",env="all",use.check=TRUE,clean.rep=TRUE)

Arguments

trait

Character. Trait of interest. The options are: "yield" (grain yield in Kg/ha), "maturity" (days to maturity), "height" (plant height in cm), "lodging" (lodging score from 1 to 5), "protein" (protein percentage in the grain), "oil" (oil percentage in the g

family

Numberic vector or "all". Which SoyNAM families to use.

env

Numberic vector or "all". Which environments to use. The environments are coded as follows: 1 (IA_2012), 2 (IA_2013), 3 (IL_2011), 4 (IL_2012), 5 (IL_2013), 6 (IN_2012), 7(IN_2013), 8 (KS_2012), 9 (KS_2013), 10 (MI_2012), 11 (MO_2012), 12 (MO_2013), 13 (N

use.check

Logical. If TRUE, it includes a control term as fixed effect in the model.

clean.rep

Logical. If TRUE, it removes genotypes that appear two or more times and replace by the expectation (mean value). Genotypes are treated as identical when the genotypes are more than 95 percent identical.

Value

• This function returns a list with four objects. A numeric vector with the BLUP solution of the phenotyes ("Phen"); the corresponding genotypes ("Gen"); a vector with the respective family ("Fam"); and a numeric vector with the number of SNPs per cromosome ("Chrom"). The output of this fuction has the exact input format for the NAM package to perform genome-wide association analysis.

Details

The algorithm start from selecting the chosen families and environment that will be used for the best linear unbias predictor (BLUP). The BLUP values are calculates based on the following model: ($Trait = Control + Environment + Genotype$). Where control is a covariate set as fixed effect based on the checks of each set (microenvironment); Environment is a random effect that represents the combination of location and year; and Genotype is the random effect associated to the lines. The BLUP values are the regression coefficients corresponding to the Genotype effect. The BLUP is calculated using the R package lme4 (Bates 2010) using REML. If checks are used as covariate (use.check=TRUE), then the best linear unbias estimator (BLUE) of the check effects is assigned to each set as a micro-environmental control. Each set had between one and five controls, including the SoyNAM parents and five other cultivars. These genotypes are normalized by environment and the BLUE of each set is calculated. All genotypes in a same set will have the same check effect.

References

Bates, D. M. (2010). lme4: Mixed-effects modeling with R. URL http://lme4.r-forge.r-project.org/book.

Examples

Test=BLUP(trait="yield",family=2:5,env=2:5)