emreml: Estimate variance components via REML using EM

In timflutre/rutilstimflutre: Timothee Flutre's personal R code

Estimate variance components via REML using EM

Description

Given the MME components, estimate variance components by the EM algorithm. References are:

• Suzuki, M. 2007. Applied Animal Breeding & Genetics. Course Notes. Obihiro University of Agriculture and Veterinary Medicine.

• Mrode, R.A. 2005. Linear Models for the Prediction of Animal Breeding Values. CAB International, Oxon, UK.

Usage

emreml(Ainv, y, X, Z, initE, initU, verbose = 0)

Arguments

Ainv	inverse of additive relationship matrix
y	vector of responses
X	incidence matrix for fixed effects
Z	incidence matrix for random effects
initE	initial value for the residual variance
initU	initial value for the additive genetic variance
verbose	verbosity level (0/1/2)

Value

vector

Note

Unknown parents should be coded as zero. Last modified by Morota on April 8, 2010.

Author(s)

Goto Morota [aut], Timothee Flutre [ctb]

See Also

makeMmeElements, makeMmeRhs, makeMmeLhs

timflutre/rutilstimflutre documentation built on Feb. 7, 2024, 8:17 a.m.