R：映射稀疏矩陣中所有條目的方法

Question

我有一個從R的Matrix包創建的稀疏矩陣。我想迭代矩陣中的每個條目並執行一個操作，將結果保存在與原始矩陣具有相同索引的另一個稀疏矩陣中。

例如，假設我有稀疏矩陣A：

1 . 1 
2 . . 
. . 4 

ColSums會是什麼樣子：

3 . 5 

RowSums會是什麼樣子：

2 
2 
4 

我想迭代A並執行此操作

(1,1) > 3*2 
(2,1) > 2*3 
(1,3) > 2*5 
(3,3) > 4*5 

創建B：

6 . 10 
6 . . 
. . 20 

我怎麼會去一個量化的方式這樣做呢？

我認爲函數foo看起來像：

B=fooMap(A,fun) 

和樂趣會是什麼樣子：

fun(row,col) = RowSums(row) * ColSums(col) 

什麼fooMap？

編輯：

我用flodel的解決方案。它使用匯總將稀疏矩陣轉換爲i，j，x數據幀，然後與&朋友一起使用來對該幀執行操作，然後將結果轉換回稀疏矩陣。使用這種技術，with/within運算符是fooMap;稀疏矩陣必須首先轉換爲i，j，x數據幀，以便可以使用/ within。

這是解決這個特殊問題的單線程。

B = with(summary(A), sparseMatrix(i=i, j=j, x = rowSums(A)[i] * colSums(A)[j])) 

Answer 1

每當我有稀疏矩陣元素的運算，我來回走了矩陣本身及其summary表示之間：

summ.B <- summary(A) 
summ.B <- within(summ.B, x <- rowSums(A)[i]*colSums(A)[j]) 
B <- sparseMatrix(i = summ.B$i, j = summ.B$j, x = summ.B$x) 
B 
# 3 x 3 sparse Matrix of class "dgCMatrix" 
#    
# [1,] 6 . 10 
# [2,] 6 . . 
# [3,] . . 20 

Answer 2

下面是一個在每一步都使用稀疏矩陣的方法。

## Load library and create example sparse matrix 
library(Matrix) 
m <- sparseMatrix(i = c(1,2,1,3), j = c(1,1,3,3), x = c(1,2,1,4)) 

## Multiply each cell by its respective row and column sums. 
Diagonal(x = rowSums(m)) %*% (1*(m!=0)) %*% Diagonal(x = colSums(m)) 
# 3 x 3 sparse Matrix of class "dgCMatrix" 
#    
# [1,] 6 . 10 
# [2,] 6 . . 
# [3,] . . 20 

（在(1*(m!=0))的1*被用來強制類「lgCMatrix」由m!=0產生回數字矩陣類「dgCMatrix」的邏輯矩陣。較長嗦（但也許更清楚）替代方案將到在它的地方使用as(m!=0, "dgCMatrix")）