web-dev-qa-db-fra.com

`which ()` fonction pour les indices matriciels

Disons que j'ai une matrice, par exemple:

> m = matrix(rep(c(0, 0, 1), 4), nrow = 4)
> m
     [,1] [,2] [,3]
[1,]    0    0    1
[2,]    0    1    0
[3,]    1    0    0
[4,]    0    0    1

Si je lance which, j'obtiens la liste des indices normaux:

> which(m == 1)
[1]  3  6  9 12

Je veux obtenir quelque chose comme des indices matriciels - chaque index contenant le numéro de ligne et de colonne:

     [,1] [,2]
[1,]    3    1
[2,]    2    2
[3,]    1    3
[4,]    4    3

Existe-t-il une fonction simple pour ce faire? De plus, elle devrait en quelque sorte contenir les noms des lignes et des colonnes:

> rownames(m) = letters[1:4]
> colnames(m) = letters[5:7]
> m
  e f g
a 0 0 1
b 0 1 0
c 1 0 0
d 0 0 1

mais je ne sais pas comment, peut-être

     [,1] [,2] [,3] [,4]
[1,]    3    1    c    e
[2,]    2    2    b    f
[3,]    1    3    a    g
[4,]    4    3    d    g

ou, peut-être retourner 2 vecteurs (pour les lignes et les colonnes), comme

c b a d
3 2 1 4

e f g g
1 2 3 3
rmatrixindexing
31
TMS

Pour votre première question, vous devez également passer arr.ind= TRUE à which:

> which(m == 1, arr.ind = TRUE)
     row col
[1,]   3   1
[2,]   2   2
[3,]   1   3
[4,]   4   3
51
David Alber

Vous ne pouvez pas mélanger numérique et alpha dans une matrice, mais vous pouvez dans un data.frame:

> indices <- data.frame(ind= which(m == 1, arr.ind=TRUE))
> indices$rnm <- rownames(m)[indices$ind.row]
> indices$cnm <- colnames(m)[indices$ind.col]
> indices
  ind.row ind.col rnm cnm
c       3       1   c   e
b       2       2   b   f
a       1       3   a   g
d       4       3   d   g
6
42-