cpowop: Square matrix power operator

In cpgen: Parallelized Genomic Prediction and GWAS

Description

This operator computes an arbitrary power of a positive definite square matrix using an Eigen-decomposition: \mathbf{X}^p = \mathbf{UD}^{p}\mathbf{U}'

Usage

X %**% power

Arguments

X	Positive definite square matrix
power	numeric scalar - desired power of X

Value

Matrix X to the power p

Examples

## Not run: 
# Inverse Square Root of a positive definite square matrix
X <- matrix(rnorm(100*5000),100,1000)

XX <- ccross(X)

XX_InvSqrt <- XX %**% -0.5

# check result: ((XX')^-0.5 (XX')^-0.5)^-1 = XX'
table(round(csolve(XX_InvSqrt %c% XX_InvSqrt),digits=2) == round(XX,digits=2) )

## End(Not run)

cpgen documentation built on May 2, 2019, 8:15 a.m.