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