aml (version 0.1-1)

aml.estimate: Estimate Genetic Values Using AML Fit

Description

This function calculate the genetic values using the result from amltest.

Usage

aml.estimate(fit, marker, response, kin)

Arguments

fit

An object generated by amltest to be used to calculate genetic values.

marker

A matrix or data frame for the markers (or genetic effects). It should be the same one used by amltest to generate fit.

response

A numerical vector of trait (phenotype) values. It should be the same one used by amltest to generate fit.

kin

The kinship matrix representing relationships between lines. It should be the same one used by amltest to generate fit.

Value

A numeric vector containing estimated genetic values for lines analyzed by amltest to generate fit.

Details

This function is used to genetic values, i.e., fitted phenotypic values using genetic marker information. It requires an adaptive mixed LASSO model has been fitted to the lines using amltest and the result is given in fit. Thus this function will only calculate genetic values for lines with observed phenotypic values. To make prediction for lines with only genetic information but no observed phenotypic values, use the function aml.pred.outside .

References

Wang, D., Eskridge, K.M. and Crossa, J. (2011) Identifying QTLs and Epistasis in Structured Plant Populations Using Adaptive Mixed LASSO. Journal of Agricultural, Biological, and Environmental Statistics, 16:170-184.

Wang, D., et al. (2012) Prediction of genetic values of quantitative traits with epistatic effects in plant breeding populations. Heredity, 109: 313-319.

See Also

amltest, aml.pred.outside.

Examples

     ## estimate genetic values for lines in the wheat data set
     data("wheat")
     clmarker<- cleanclust(wheat$marker, nafrac=0.2, mafb=0.1, corbnd=0.5, method="complete")
     intermat <- epigen(wheat$y, clmarker$newmarker, wheat$A, numkeep=100, selectvar=30, 
                          corbnd=0.5, mafb=0.04)
     resepi <- amltest(wheat$y, intermat$effects, wheat$A, numkeep=80, selectvar=40)
     predall<-aml.estimate(resepi, intermat$effects, wheat$y, wheat$A)