randTraceHutchinson {randPedPCA}R Documentation

Trace estimation for sparse L inverse matrices

Description

Using Hutchinson's method

Usage

randTraceHutchinson(L, numVectors)

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.

Value

a scalar


[Package randPedPCA version 1.1.3 Index]