R稀疏矩阵电源

Question

我正在使用Matrix软件包来创建大量零（〜14000x14000）稀疏矩阵。有谁知道计算这个矩阵的功率的最好方法吗？

我试过A_pow2 = A％^％2但我得到错误：A％^％2中的错误：不是矩阵。下面是返回相同的错误一个简单的例子：

A = matrix(3,2,2) 
A = Matrix(A,sparse=TRUE) 
Apow2 = A%^%2 

Answer 1

（编辑感谢@罗兰的评论）

自定义功能或许能够解决您的问题。每?expm::`%^%`

Compute the k-th power of a matrix. Whereas x^k computes element wise powers, x %^% k corresponds to k - 1 matrix multiplications, x %*% x %*% ... %*% x.

文档，我们可以写一个新的管道符执行乘法K-1倍。不知道它将如何扩展，但它在更小的例子中起作用。

> library(Matrix) 
> library(expm) 
> A = matrix(3,2,2) 
> B = Matrix(A,sparse=TRUE) 
> 
> # changed lapply to rep list 
> `%^^%` = function(x, k) Reduce(`%*%`, rep(list(x), k)) 
> # per Roland for loop approach will be better on memory 
> `%^^%` = function(x, k) {for (i in 1:(k - 1)) x <- x %*% x; x} 
> 
> as.matrix(B%^^%2) 
    [,1] [,2] 
[1,] 18 18 
[2,] 18 18 
> A%^%2 
    [,1] [,2] 
[1,] 18 18 
[2,] 18 18