powerofpsdmatrix {compositions}R Documentation

power transform of a matrix

Description

Computes a power of a positive semidefinite symmetric matrix.

Usage

          powerofpsdmatrix( M , p,...)
          

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 a singular value decomposition of M.

Author(s)

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

References

See Also

Examples

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

[Package compositions version 0.9-11 Index]