Function for applying R inverse to AtG in REML EM algorithm

Source: R/em_reml.R

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