synbreed (version 0.9-2)

predict.gpMod: Prediction for genomic prediction models.

Description

S3 predict method for objects of class gpMod. A genomic prediction model is used to predict the genetic performance for e.g. unphenotyped individuals using an object of class gpMod estimated by a training set.

Usage

## S3 method for class 'gpMod':
predict(object,newdata,...)

Arguments

object

object of class gpMod which is the model used for the prediction.

newdata

for model="BL" and "BRR" an object of class gpData with the marker data of the unphenotyped individuals. For model="BLUP" a character vector with the names of the individuals to predict.

...

not used

Value

• a named vector with the predicted genetic values for all individuals in newdata.

Details

For models, model="RR" and "BL", the prediction for the unphenotyped individuals is given by $$\hat{g}'=\hat{\mu} + W'\hat{m}$$ with the estimates taken from the gpMod object. For the prediction using model="BLUP", the full relationship matrix including individuals of the training set and the prediction set must be specified in the gpMod. This model is used to predict the unphenotyped individuals of the prediction set by solving the corresponding mixed model equations using the variance components of the fit in gpMod.

References

Henderson C (1977) Best linear unbiased prediction of breeding values not in the model for records. Journal of Dairy Science 60:783-787 Henderson CR (1984). Applications of linear models in animal breeding. University of Guelph.

See Also

gpMod

Examples

# Example from Henderson (1977)
dat <- data.frame(y=c(132,147,156,172),time=c(1,2,1,2),animal=c(1,2,3,4))
ped <- create.pedigree(ID=c(6,5,1,2,3,4),Par1=c(0,0,5,5,1,6),Par2=c(0,0,0,0,6,2))
gp <- create.gpData(pheno=dat,pedigree=ped)
A <- kin(gp,ret="add")

# assuming h2=sigma2u/(sigma2u+sigma2)=0.5
# no REML fit possible due to the limited number of observations
y <- c(132,147,156,172)
names(y) <- paste(1:4)
mod1 <- list(fit=list(sigma=c(1,1)),kin=A,model="BLUP",y=y,m=NULL)
class(mod1) <- "gpMod"
predict(mod1,c("5","6"))

# prediction 'by hand'
X <- matrix(1,ncol=1,nrow=4)
Z <- diag(6)[-c(1,2),]
AI <- solve(A)
RI <- diag(4)

res <- MME(X,Z,AI,RI,y)
res$b + res$u[1:2]