Learn R Programming

statgenQTLxT (version 1.0.3)

updatePrec: Helper function for updating precision matrix.

Description

Helper function for updating the precision matrices in the EMFA algorithm.

Usage

updatePrec(m, nc, omega, w, p, wNew, pNew, cNew, het, maxDiag)

Arguments

m

An integer, the order of the model.

nc

An integer, the number of traits or genotypes.

omega

A computed matrix for the current step in the algoritm.

w

A model matrix for the current step in the algorithm.

p

A model matrix for the current step in the algorithm.

wNew

A pointer to the updated model matrix for w.

pNew

A pointer to the updated model matrix for p.

cNew

A pointer to the updated matrix c.

het

Should an extra diagonal part be added in the model for the precision matrix.

maxDiag

A numerical value for the maximum value of sigma2.