Function for applying R inverse to AtG in REML EM algorithm

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

Arguments

G

Matrix of genotype markers (size n x L)

AtG

AtG from pre-computation

GtAAtG

GtAAtG from pre-computation

tau

Variance component from G main effect

sigma

Variance component from model noise epsilon

Value

Matrix resulting from left multiplication of Rinv with input matrix AtG