blup: Function to calculate breeding values using an animal model

In pedigree: Pedigree Functions

Function to calculate breeding values using an animal model

Description

Fit an animal model to data, use a given variance ratio (α = \frac{σ^2_e}{σ^ 2_a}). Calculate inverse of the additive genetic relationship matrix using function makeInv() of this package.

Usage

blup(formula, ped, alpha,trim = FALSE)

Arguments

formula	formula of the model, do not include the random effect due to animal (generally ID).
ped	data.frame with columns corresponding to ID, SIRE, DAM and the columns in the formula.
alpha	Variance ratio (\frac{σ^2_e}{σ^ 2_a}).
trim	If TRUE, trims the pedigree using the available phenotype data using function trimPed.

Value

Vector of solutions to the model, including random animal effects.

See Also

SamplePedigree, gblup, makeAinv,blup

Examples

example(gblup)
sol <- blup(P~1,ped = ped,alpha = 1/h2 - 1)

pedigree documentation built on Aug. 14, 2022, 1:06 a.m.