我们可以使用sparseMatrix
或spMatrix
从索引和非零元素的值构造稀疏矩阵。有没有函数将稀疏矩阵转换回索引和所有非零元素的值?例如
i <- c(1,3,5); j <- c(1,3,4); x <- 1:3
A <- sparseMatrix(i, j, x = x)
B <- sparseToVector(A)
## test case:
identical(B,cbind(i,j,x))
是否有任何功能与sparseToVector
做类似的工作?
答案 0 :(得分:9)
您的矩阵A
采用稀疏压缩格式(类dgCMatrix
)。您可以通过
A.nc <- as (A, "dgTMatrix")
或者,您可以在giveCsparse = TRUE
来电中指定sparseMatrix
。
dgTMatrix
的三元组形式基本上包含您在广告资源i
,j
和x
中寻找的所有内容,只有i
和使用基于0的偏移完成j
索引:
> str (A.nc)
Formal class 'dgTMatrix' [package "Matrix"] with 6 slots
..@ i : int [1:3] 0 2 4
..@ j : int [1:3] 0 2 3
..@ Dim : int [1:2] 5 4
..@ Dimnames:List of 2
.. ..$ : NULL
.. ..$ : NULL
..@ x : num [1:3] 1 2 3
..@ factors : list()
> cbind (i = A.nc@i + 1, j = A.nc@j + 1, x = A.nc@x)
i j x
[1,] 1 1 1
[2,] 3 3 2
[3,] 5 4 3
> all (cbind (i = A.nc@i + 1, j = A.nc@j + 1, x = A.nc@x) == cbind (i, j, x))
[1] TRUE
答案 1 :(得分:8)
summary(A)
# 5 x 4 sparse Matrix of class "dgCMatrix", with 3 entries
# i j x
# 1 1 1 1
# 2 3 3 2
# 3 5 4 3
您可以轻松传递到as.data.frame
或as.matrix
:
sparseToVector <- function(x)as.matrix(summary(x))
B <- sparseToVector(A)
## test case:
identical(B,cbind(i,j,x))
# [1] TRUE
答案 2 :(得分:2)
将which
与arr.ind
:
idx <- which(A != 0, arr.ind=TRUE)
cbind(idx, A[idx])
# [,1] [,2] [,3]
# [1,] 1 1 1
# [2,] 3 3 2
# [3,] 5 4 3