轉換矩陣中的R

三個定義的列鑑於m：轉換矩陣中的R

m <- structure(c(5, 1, 3, 2, 1, 4, 5, 2, 5, 1, 1, 5, 1, 4, 0, 4, 5, 
5, 3, 2, 0, 0, 3, 0, 3, 2, 3, 0, 0, 0, 0, 0, 0, 0, 0), .Dim = c(7L, 
5L)) 

    # [,1] [,2] [,3] [,4] [,5] 
# [1,] 5 2 0 0 0 
# [2,] 1 5 4 3 0 
# [3,] 3 1 5 0 0 
# [4,] 2 1 5 3 0 
# [5,] 1 5 3 2 0 
# [6,] 4 1 2 3 0 
# [7,] 5 4 0 0 0

考慮元件1，它出現在5行(2, 3, 4, ,5, 6)和各列方向指數(1, 2, 2, 1, 2)。我想有以下幾點：

再舉一個例子，考慮元素2，它出現在4行(1, 4, 5, 6)和相應逐列指數(2, 1, 4, 3)我們有：

我想要的是所有1-5的n*3矩陣。 優選地，在基礎R

來源

2016-07-15 989

我們可以用與whicharr.ind=TRUE

cbind(val= 1, which(m==1, arr.ind=TRUE)) 
#  val row col 
#[1,] 1 2 1 
#[2,] 1 5 1 
#[3,] 1 3 2 
#[4,] 1 4 2 
#[5,] 1 6 2

對於多情況下，如@RHertel提到

for(i in 1:5) print(cbind(i,which(m==i, arr.ind=TRUE)))

或用lapply

do.call(rbind, lapply(1:2, function(i) { 
     m1 <-cbind(val=i,which(m==i, arr.ind=TRUE)) 
     m1[order(m1[,2]),]})) 
#  val row col 
#[1,] 1 2 1 
#[2,] 1 3 2 
#[3,] 1 4 2 
#[4,] 1 5 1 
#[5,] 1 6 2 
#[6,] 2 1 2 
#[7,] 2 4 1 
#[8,] 2 5 4 
#[9,] 2 6 3

由於OP提到約base R解決方案，上述將有所幫助。但是，在情況下，如果有人想一個緊湊的解決方案，

library(reshape2) 
melt(m)

，然後子集的感興趣的值。

來源

2016-07-15 15:11:58 akrun

一個方便的方式來改造它是使用sparseMatrix從Matrix庫，由於所需輸出非常接近稀疏矩陣的表示：

library(Matrix) 
summary(Matrix(m, sparse = T)) 
# 7 x 5 sparse Matrix of class "dgCMatrix", with 23 entries 
# i j x 
# 1 1 1 5 
# 2 2 1 1 
# 3 3 1 3 
# 4 4 1 2 
# 5 5 1 1 
# 6 6 1 4 
# 7 7 1 5 
# 8 1 2 2 
# 9 2 2 5 
# 10 3 2 1 
# 11 4 2 1 
# 12 5 2 5 
# 13 6 2 1 
# 14 7 2 4 
# 15 2 3 4 
# 16 3 3 5 
# 17 4 3 5 
# 18 5 3 3 
# 19 6 3 2 
# 20 2 4 3 
# 21 4 4 3 
# 22 5 4 2 
# 23 6 4 3

爲了看得更清楚：

summary(Matrix(m, sparse = T)) %>% dplyr::arrange(x) 
# i j x 
# 1 2 1 1 
# 2 5 1 1 
# 3 3 2 1 
# 4 4 2 1 
# 5 6 2 1 
# 6 4 1 2 
# 7 1 2 2 
# 8 6 3 2 
# 9 5 4 2 
# 10 3 1 3 
# 11 5 3 3 
# 12 2 4 3 
# 13 4 4 3 
# 14 6 4 3 
# 15 6 1 4 
# 16 7 2 4 
# 17 2 3 4 
# 18 1 1 5 
# 19 7 1 5 
# 20 2 2 5 
# 21 5 2 5 
# 22 3 3 5 
# 23 4 3 5

來源

2016-07-15 15:18:03 Psidom

只需使用row和col即可。

> data.frame(m=as.vector(m), row=as.vector(row(m)), col=as.vector(col(m))) 
    m row col 
1 5 1 1 
2 1 2 1 
3 3 3 1 
4 2 4 1 
5 1 5 1 
...

子集，排序，並且如期望的打印。

> tmp <- out[order(out$m, out$row), ] 
> print(subset(tmp, m==1), row.names=FALSE) 
m row col 
1 2 1 
1 3 2 
1 4 2 
1 5 1 
1 6 2

來源

2016-07-15 17:32:42 Aaron

轉換矩陣中的R

回答

相關問題