Vinv {SEAGLE}

Function for applying V inverse in Algorithm 1

Description

This function applies V inverse via the Woodbury matrix identity

Usage

Vinv(G, qrM, tau_over_sigma, sigma, RHS)

Arguments

G

Matrix of genotype markers (size n x L)

qrM

Pre-computation for LxL linear system solve

tau_over_sigma

Tau over sigma from precomputation

sigma

Variance component from model noise epsilon

RHS

Matrix or vector on right-hand side of V inverse

Value

Matrix or vector resulting from left multiplication of Vinv with input RHS