powerofmatrix: power transform of a matrix

Description Usage Arguments Details Value Author(s) Examples

Description

Computes the power of a positive semi-definite symmetric matrix.

Usage

1
2

Arguments

M

a matrix, preferably symmetric

p

a single number giving the power

...

further arguments to the singular value decomposition

Details

for a symmetric matrix the computed result can actually be considered as a version of the given power of the matrix fullfilling the relation:

M^pM^q=M^{p+q}

The symmetry of the matrix is not checked.

Value

U%*% D^p %*% t(P) where the UDP is the singular value decomposition of M.

Author(s)

K.Gerald v.d. Boogaart http://www.stat.boogaart.de

Examples

1
2
3
data(SimulatedAmounts)
d <- ilr(sa.lognormals)
var( d %*% powerofpsdmatrix(var(d),-1/2)) # Unit matrix


Search within the compositions package
Search all R packages, documentation and source code

Questions? Problems? Suggestions? or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.