[R] Matrix negative fraction power
Berend Hasselman
bhh at xs4all.nl
Sun Mar 11 18:18:41 CET 2012
On 11-03-2012, at 17:52, Spencer Graves wrote:
> If my memory is correct, the archives of this list contains several discussions of round off error problems associated with different methods for computing things like this. The "Matrix" package (part of the base distribution) contains a function "expm", whose help file says, "The expm package contains newer (partly faster and more accurate) algorithms for expm() and includes logm and sqrtm." This suggests to me that the most numerically stable way to get an arbitrary power p of a matrix M in R might be as follows:
>
>
> expm(p*logm(M))
>
I have just tested and it should be this
library(expm)
expm((-1/2)*logm(A))
expm((-3/2)*logm(A))
Berend
More information about the R-help
mailing list