Last chance! 50% off unlimited learning
Sale ends in
Computes the power of a positive semi-definite symmetric matrix.
powerofpsdmatrix( M , p,...)
U%*% D^p %*% t(P)
where the UDP
is the singular value
decomposition of M.
a matrix, preferably symmetric
a single number giving the power
further arguments to the singular value decomposition
K.Gerald v.d. Boogaart http://www.stat.boogaart.de
for a symmetric matrix the computed result can actually be considered as
a version of the given power of the matrix fullfilling the relation:
data(SimulatedAmounts)
d <- ilr(sa.lognormals)
var( d %*% powerofpsdmatrix(var(d),-1/2)) # Unit matrix
Run the code above in your browser using DataLab