Function for applying R inverse to u in REML EM algorithm

Source: R/em_reml.R

Function for applying R inverse to u in REML EM algorithm

Rinv.u(G, AtG, GtAAtG, GtAu, u, tau, sigma)

Arguments

G

Matrix of genotype markers (size n x L)
AtG

AtG from precomputation
GtAAtG

GtAAtG from precomputation
GtAu

GtAu from precomputation
u

u=Aty from REML EM
tau

Variance component from G main effect
sigma

Variance component from model noise epsilon

Value

Vector resulting from left multiplication of Rinv with input vector u