14

sparseMatrixorを使用して、非ゼロ要素のインデックスと値からスパース行列を作成できます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か？

4

3 に答える 3

9

マトリックスAはスパース圧縮形式 (クラスdgCMatrix) です。圧縮されていないスパース形式に強制することができます

A.nc <- as (A, "dgTMatrix")

または、呼び出しで指定giveCsparse = TRUEすることもできます。sparseMatrix

のトリプレット形式には基本的に、スロット、、およびでdgTMatrix探しているすべてが含まれており、インデックスは 0 ベースのオフセットで行われます。ijxij

> 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
于 2013-04-06T12:13:04.493 に答える
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
于 2013-04-06T12:41:53.620 に答える
3

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
于 2013-04-06T10:21:45.453 に答える