QTLRel (version 0.1)

blup: Best Linear Unbiased Prediction

Description

Estimate the best linear unbiased prediction (BLUP) for various effects in the model.

Usage

blup(object)

Arguments

object

an object from estVC or aicVC.

Value

• fixed
  BLUP for fixed effects.
• AA,DD,...
  BLUP for random (genetic) variance components "AA", "DD", ...
• EE
  BLUP for residual effect.

See Also

estVC and aicVC.

Examples

y<- rnorm(100)
x<- matrix(1,nrow=100,ncol=1)
v<- cov(matrix(rnorm(10^6),ncol=100))

o<- estVC(y, x, v = list(AA=v,DD=NULL,HH=NULL,AD=NULL,
   MH=NULL,EE=diag(100)))
b<- blup(o)