Discussion :

[jack-no]

Bonsoir cher(e)s ami(e)s,
Depuis un moment maintenant je cherche à calculer le mode (la fréquence la plus élevée) d'un vecteur.
Je n'ai pas trouvé de fonction toute faite ( ce qui est dommage pour un programme aussi réputé en stat).
Bref, je l'ai codé moi-même (ce n'était pas très compliqué).
Le problème est que lorsque je veux utiliser le résultat pour faire autre chose avec j'ai une erreur.
Voici mon code:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
mod <- function(x)
{
  m<- unique(x)
  L<- length(m)
  w<- NULL
 
  for(i in 1:L)
  w<-c(w,length(x==m[i]))
 
  h<- sort(w, decreasing =TRUE)[1]
  res<- x[grep(h,w)][1]
 
  return(res)
}

Et voici l'utilisation:

Code :

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

c <-mod( vecteur )

J'obtiens:

Code :

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

1
2
3
c
[1] 2
Levels: 1 2 3

Dans mon cas c=2 (Je ne comprends pas pourquoi il m'affiche des "levels");
Je fais ensuite:

Code :

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

comb <- Prono(c,comb,S,2,1)

J'ai une erreur; Alors qu'en remplaçant c par 2 directement comme ceci:

Code :

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

comb <- Prono(2,comb,S,2,1)

Ça fonctionne! Franchement R est en train de me tourner en bourrique!

[jack-no]

Alléluia! Le problème est résolu!
Dans ma fonction Prono(), je faisais un test du type:

Code :

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

if(!x)

Que j'ai remplacé par:

Code :

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

if(x !=0)

Et comme par miracle plus de soucis.

[juliatheric]

Bonsoir,
Même si le problème est résolu, j'aimerais clarifier deux points.

Bonsoir cher(e)s ami(e)s,
Depuis un moment maintenant je cherche à calculer le mode (la fréquence la plus élevée) d'un vecteur.
Je n'ai pas trouvé de fonction toute faite ( ce qui est dommage pour un programme aussi réputé en stat). [...]

R propose bel et bien des outils pour ça. Les fréquences du contenu d'un objet donné peuvent être calculées avec la fonction table. Après, il suffit d'extraire le maximum du résultat.

Code :

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

> x <- sample(15, replace=TRUE); x

 [1] 10  4 12 15  7  6  6  7 14 15  4  2  7  5 13

Code :

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

1
2
> y <- table(x)
> y; y[which.max(y)]

x
 2  4  5  6  7 10 12 13 14 15 
 1  2  1  2  3  1  1  1  1  2 
7 
3

[...]
Dans mon cas c=2 (Je ne comprends pas pourquoi il m'affiche des "levels");
[...]

Les "levels" correspondent aux modalités (en Statistique). C'est un type d'objets au même titre que numeric, character, ...