R语言矩阵的乘方如何编程

stephen7 发表于 2012-11-27 02:12
rmutil包里有一个 %^%函数，可以做方阵的幂，但是这个包在R里不能直接安装，要到作者的网站去下载，然后解 ...

只有这个函数用不了啊，说是不存在

http://cos.name/cn/topic/5619#post-5619
统计之都的回复，写的不错，可以查看。
PS：统计之都谈论的R问题比这边深刻很多。

一般方针阵先对角化，然后特征值矩阵自己的幂方，再把特征向量乘回去，估计会比较快。

循环或者连城都可以的

matpower=function(da,n){
  dat1=diag(dim(dat)[1])
  for(i in 1:n){
    dat1=dat1%*%da
  }
  return(dat1)
}

上官雨凌 发表于 2014-10-4 09:52
只有这个函数用不了啊，说是不存在

请问前辈，这个问题你解决了吗？加载了rmutil包，没有发现 %^% 函数