Jak przekształcić rzadką macierz w macierz indeksu i wartość niezerowego elementu

Question

Możemy zbudować rzadką macierz z indeksu i wartość niezerowego elementu z sparseMatrix lub spMatrix. Czy istnieje jakakolwiek funkcja przekształca rzadką macierz z powrotem w indeks i wartość wszystkich elementów niezerowych? Na przykład:

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

Czy istnieje jakaś funkcja wykonująca podobną pracę, jak sparseToVector?

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 

który możesz z łatwością przejść pod numer as.data.frame lub as.matrix:

---

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

Answer 2

Zastosowanie which z 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

Twój matryca A jest w rzadkiej skompresowanym formacie (klasa dgCMatrix). Można zmusić go do non-skompresowanym formacie rozrzedzony przez

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

Albo mogłeś określony giveCsparse = TRUE w zaproszeniu sparseMatrix.

Formularz triplet dgTMatrix zasadniczo zawiera wszystko, czego szukasz w gniazdach i, j, a x tylko i i j indeksowanie odbywa się z offsetu 0 opartymi:

> 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