我试图找到一种有效的方法来检索 R 中稀疏矩阵的非零上三角元素的列表/向量/数组。
例如:

    library(igraph)
    Gmini <- as.directed(graph.lattice(c(3,5)))
    GminiMat <- sparseMatrix(i=get.edgelist(Gmini)[,1],j=get.edgelist(Gmini)[,2],x=1:length(E(Gmini)))
    GminiMat

    15 x 15 sparse Matrix of class "dgCMatrix"

    [1,]  .  1  .  2  .  .  .  .  .  .  .  .  .  .  .
    [2,] 23  .  3  .  4  .  .  .  .  .  .  .  .  .  .
    [3,]  . 25  .  .  .  5  .  .  .  .  .  .  .  .  .
    [4,] 24  .  .  .  6  .  7  .  .  .  .  .  .  .  .
    [5,]  . 26  . 28  .  8  .  9  .  .  .  .  .  .  .
    [6,]  .  . 27  . 30  .  .  . 10  .  .  .  .  .  .
    [7,]  .  .  . 29  .  .  . 11  . 12  .  .  .  .  .
    [8,]  .  .  .  . 31  . 33  . 13  . 14  .  .  .  .
    [9,]  .  .  .  .  . 32  . 35  .  .  . 15  .  .  .
    [10,]  .  .  .  .  .  . 34  .  .  . 16  . 17  .  .
    [11,]  .  .  .  .  .  .  . 36  . 38  . 18  . 19  .
    [12,]  .  .  .  .  .  .  .  . 37  . 40  .  .  . 20
    [13,]  .  .  .  .  .  .  .  .  . 39  .  .  . 21  .
    [14,]  .  .  .  .  .  .  .  .  .  . 41  . 43  . 22
    [15,]  .  .  .  .  .  .  .  .  .  .  . 42  . 44  .

所以理想情况下我想创建一个函数 getUpper(mat) 使得
getUpper(GminiMat)
将产生 1:22 的向量(GminiMat 的上三角非零项)

理想情况下,我需要一种内存和速度相当高效的方法,因为我可能需要将其应用于大型系统(例如,矩阵可能来自每个维度有数百个节点的多维点阵)。

最佳答案

Matrix 包中有一个 triu 方法,它将在保留稀疏性的同时返回上三角形:

triu(GminiMat)
15 x 15 sparse Matrix of class "dtCMatrix"

 [1,] . 1 . 2 . . .  .  .  .  .  .  .  .  .
 [2,] . . 3 . 4 . .  .  .  .  .  .  .  .  .
 [3,] . . . . . 5 .  .  .  .  .  .  .  .  .
 [4,] . . . . 6 . 7  .  .  .  .  .  .  .  .
 [5,] . . . . . 8 .  9  .  .  .  .  .  .  .
 [6,] . . . . . . .  . 10  .  .  .  .  .  .
 [7,] . . . . . . . 11  . 12  .  .  .  .  .
 [8,] . . . . . . .  . 13  . 14  .  .  .  .
 [9,] . . . . . . .  .  .  .  . 15  .  .  .
[10,] . . . . . . .  .  .  . 16  . 17  .  .
[11,] . . . . . . .  .  .  .  . 18  . 19  .
[12,] . . . . . . .  .  .  .  .  .  .  . 20
[13,] . . . . . . .  .  .  .  .  .  . 21  .
[14,] . . . . . . .  .  .  .  .  .  .  . 22
[15,] . . . . . . .  .  .  .  .  .  .  .  .

关于r - 获取 R 中稀疏矩阵的上三角非零元素,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/23395487/

10-12 17:22