r - R中的矩阵幂

Question

尝试计算R中矩阵的幂次，我发现expm包实现了运算符%^%。

所以x％^％k计算矩阵的第k次幂。

> A<-matrix(c(1,3,0,2,8,4,1,1,1),nrow=3)

> A %^% 5
      [,1]  [,2] [,3]
[1,]  6469 18038 2929
[2,] 21837 60902 9889
[3,] 10440 29116 4729

但是，令我惊讶的是:

> A
     [,1] [,2] [,3]
[1,]  691 1926  312
[2,] 2331 6502 1056
[3,] 1116 3108  505

不知何故，初始矩阵A已更改为A％^％4 !!!

您如何执行矩阵幂运算？

Answer 1

我已经修复了R-forge来源(“expm”程序包中)的错误，
svn版本53.-> expm R-forge page
由于某些原因，该网页仍显示修订版52，因此可能尚未显示以下内容
解决您的问题(但应在24小时之内):

 install.packages("expm", repos="http://R-Forge.R-project.org")

否则，直接获取svn版本，并自行安装:

 svn checkout svn://svn.r-forge.r-project.org/svnroot/expm

感谢“gd047”，他通过电子邮件提醒我该问题。
请注意，R-forge还具有自己的错误跟踪功能。
马丁特