Discussion :

[ricoben]

je cherche une manière élégante de trouver le couplet ligne/colonne d'un tableau de contingence qui correspond au maximum de la valeur du tableau,
sans jouer sur des boucles emboitées.

je pensais à utiliser une syntaxe type apply...
quelqu'un a t il déjà été confronté au pb?
Merci

[pitipoisson]

Bonjour,

Si tu veux utiliser des apply, toute la difficulté c'est de traiter les cas d'égalités... ça oblige à de petites contorsions :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
> data <- data.frame(a=letters[1:5][round(runif(30, 0.2, 1) * 5)], B=LETTERS[1:4][round(runif(30, 0.25, 1) * 4)])
> tab <- table(data$a, data$B)
> tab
 
    A B C D
  a 1 1 3 0
  b 0 2 3 2
  c 2 3 1 1
  d 0 2 2 0
  e 3 2 2 0
 
> x <- apply(tab, 1, function(col) c(which.max(col), max(col))); colnames(x) <- NULL
 
> y <- apply(tab, 2, function(row) c(which.max(row), max(row))); colnames(y) <- NULL
 
> c(row=which.max(x[2, ]), col=x[1, which.max(x[2, ])])   ## Priorité aux colonnes.
row col 
  1   3 
 
> c(row=y[1, which.max(y[2, ])], col=which.max(y[2, ]))   ## Priorité aux lignes.
row col 
  5   1

Ça ne te donne que la position de la première occurrence... avec priorité aux lignes ou aux colonnes, au choix.

Tu peux également faire plus simple sachant que which.max() fonctionne également sur les matrices :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
> c(row = ifelse(which.max(tab) %% nrow(tab), which.max(tab) %% nrow(tab), nrow(tab)), col = ceiling(which.max(tab) / nrow(tab)))
row col 
  5   1

(équivalent priorité aux lignes pour les égalités).
Mais attention, je pense que ça doit donner des résultats erronés si la matrice a été créée avec l'option byrow=TRUE !

[ricoben]

Merci ca va m'aider;