Wie man eine dünn besetzte Matrix in eine Matrix aus Index und Wert eines Nicht-Null-Elements umwandelt

Question

Wir können eine dünn besetzte Matrix aus einem Index und einem Wert eines Nicht-Null-Elements mit sparseMatrix oder spMatrix konstruieren. Gibt es irgendeine Funktion, die eine dünn besetzte Matrix zurück in einen Index und einen Wert aller Nicht-Null-Elemente umwandelt? Zum Beispiel

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)) 

Gibt es eine Funktion, um eine ähnliche Arbeit wie sparseToVector tun?

Answer 1

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 

die Sie können leicht zu as.data.frame oder as.matrix gehen:

---

sparseToVector <- function(x)as.matrix(summary(x)) 
B <- sparseToVector(A) 
## test case: 
identical(B,cbind(i,j,x)) 
# [1] TRUE 

Answer 2

Verwendung which mit 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 

Answer 3

Ihre Matrix A ist in der spärlichen komprimierten Format (Klasse dgCMatrix). Sie können es in den nicht komprimierten spärlichen Format von

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

coerce Oder Sie giveCsparse = TRUE im sparseMatrix Aufruf angegeben haben könnte.

Die Triplett-Form dgTMatrix grundsätzlich alle für Sie suchen enthält in Schlitze i, j und x, nur i und j Indizierung mit 0 basierten Offsets getan:

> 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