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