qui.max des liens de méthode dans la R

which.max et which.min sera de retour le plus petit indice de la valeur max ou min si il y a des liens.

Est-il un moyen de contourner cela afin que le plus grand indice est retourné sans affectent l'efficacité de la fonction?

max.col a cette fonctionnalité, mais j'ai affaire à un vecteur non d'une matrice.

OriginalL'auteur by0 | 2013-03-26

Vous pouviez faire comme ceci:

x<-c(1,2,1,4,3,4)
#identical to which.max, except returns all indices with max
which(x==max(x)) 
[1] 4 6
z<-which(x==max(x))
z[length(z)]
[1] 6
#or with tail
tail(which(x==max(x)),1)
[1] 6

edit:

Ou, vous pouvez également utiliser max.col fonction des vecteurs comme ceci:

max.col(t(x),"last")
[1] 6
#or
max.col(matrix(x,nrow=1),"last")
[1] 6

edit: Certains d'analyse comparative:

x<-sample(1:1000,size=10000,replace=TRUE)
library(microbenchmark)
microbenchmark(which.max(x),{z<-which(x==max(x));z[length(z)]}, 
     tail(which(x==max(x)),1),max.col(matrix(x,nrow=1),"last"),
     max.col(t(x),"last"),which.max(rev(x)),times=1000)
Unit: microseconds
                                             expr     min      lq  median      uq       max neval
                                     which.max(x)  29.390  30.323  30.323  31.256 17550.276  1000
 {     z <- which(x == max(x))     z[length(z)] }  40.586  42.452  42.919  44.318   631.178  1000
                      tail(which(x == max(x)), 1)  57.380  60.646  61.579  64.844   596.657  1000
             max.col(matrix(x, nrow = 1), "last") 134.353 138.085 139.485 144.383   710.949  1000
                            max.col(t(x), "last") 116.159 119.425 121.291 125.956   729.610  1000
                                which.max(rev(x))  89.569  91.435  92.368  96.566   746.404  1000

Donc toutes les méthodes semblent être plus lent que l'original (qui donne le mauvais résultat), mais z <- which(x == max(x));z[length(z)] semble être l'option la plus rapide de ces.

L'analyse comparative est magnifique!!! Merci, je pense que je vais aller avec l'option 2! ` z <- qui(x == max(x)) z[longueur(z)]`

OriginalL'auteur Jouni Helske

4

Vous pouvez inverser x
```
which.max(rev(x))
which.min(rev(x))
```
OriginalL'auteur mnel
4

La which la fonction a un 'arr.ind paramètre normalement défini à FALSE, mais utilement la valeur TRUE dans ce cas:
```
x <- sample(1:20, 50, repl=TRUE)

> which(x==max(x), arr.ind=TRUE)
[1] 11 23
> tail(which(x==max(x), arr.ind=TRUE) , 1)
[1] 23
```
À l'aide de l'arr.ind argument est particulièrement utile avec la matrice ou un tableau de structures, mais il ne fonctionne pas avec le atomique vecteurs.

DWin, est-il une différence dans la mise en arr.ind=TRUE pour les vecteurs?
Je ne le pense pas. Je devrais probablement juste utilisé which. Parfois, j'ai trop réfléchir.

OriginalL'auteur 42-

De s'étendre sur les Jouni de réponse, vous pouvez à la place utiliser max sur le résultat de which:

x <- c(1, 2, 1, 4, 3, 4)
which(x == max(x)) 
[1] 4 6
max(which(x == max(x)))
[1] 6

Benchmarking:

x <- sample(1:1000, size = 10000, replace = TRUE)
library(microbenchmark)
microbenchmark(which.max(x), {z <- which(x == max(x)); z[length(z)]}, 
tail(which(x == max(x)), 1), max.col(matrix(x, nrow = 1), "last"),
max.col(t(x), "last"), which.max(rev(x)), max(which(x == max(x))), times = 1000)
Unit: microseconds
expr     min      lq       mean  median      uq      max neval
which.max(x)   6.322   6.717   7.171838   7.112   7.112   40.297  1000
{     z <- which(x == max(x))     z[length(z)] }  27.260  28.445  37.126964  28.840  29.630 2276.346  1000
tail(which(x == max(x)), 1)  35.952  37.927  45.198484  38.718  40.298 1005.038  1000
max.col(matrix(x, nrow = 1), "last") 160.791 162.766 181.698171 163.557 169.087 1688.494  1000
max.col(t(x), "last")  84.149  86.124 100.249921  86.915  89.680 1230.618  1000
which.max(rev(x))  53.729  55.310  69.442985  56.100  57.680 1076.149  1000
max(which(x == max(x)))  26.865  27.655  35.552256  28.050  28.841 1029.137  1000

OriginalL'auteur BobbyBean

Vous devez vous connecter pour publier un commentaire.