Fast analytical REML estimation for single variance component

For a single variance component model, the REML estimate of λ = σe²/σu² has a closed-form solution that can be computed efficiently.

Arguments

omega

Transformed response vector Q'y

theta

Transformed eigenvalues

tol

Convergence tolerance for Newton iterations

max_iter

Maximum Newton iterations

Value

Estimated variance ratio λ