Learn R Programming

randPedPCA (version 1.0.1)

randTraceHutchinson: Trace estimation for sparse L inverse matrices

Description

Using Hutchinson's method

Usage

randTraceHutchinson(L, numVectors)

Value

a scalar

Arguments

L

A pedigree's L inverse matrix

numVectors,

an integer specifying how many random vectors to use

If you do not have a good reason to do otherwise, use the function hutchpp instead.

The higher numVectors, the higher the accuracy and the longer the running time. Accuracy can be estimated with the function getNumVectorsHutchinson.