Description Usage Arguments Details Value References See Also Examples
This function calculate the genetic values using the result from amltest
.
1 2 | aml.estimate(fit, marker, response, kin)
|
fit |
An object generated by |
marker |
A matrix or data frame for the markers (or genetic effects). It should be the same one used by |
response |
A numerical vector of trait (phenotype) values. It should be the same one used by |
kin |
The kinship matrix representing relationships between lines. It should be the same one used by |
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
.
A numeric vector containing estimated genetic values for lines analyzed by amltest
to generate fit
.
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.
1 2 3 4 5 6 7 8 9 | ## 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)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.