blup | R Documentation |
Estimate the best linear unbiased prediction (BLUP) for various effects in the model.
blup(object)
object |
An object from |
fixed |
BLUP for fixed effects. |
R, etc. |
BLUP for random effects. |
estVC
and aicVC
.
data(miscEx)
## Not run:
# only consider additive genetic variance component
pheno<- pdatF8[!is.na(pdatF8$bwt) & !is.na(pdatF8$sex),]
ii<- match(rownames(pheno), rownames(gmF8$AA))
v<- list(A=gmF8$AA[ii,ii],D=gmF8$DD[ii,ii])
vc<- estVC(y=pheno$bwt, x=pheno$sex, v=v)
b<- blup(vc)
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.