Forum de Statistiques
Vous souhaitez réagir à ce message ? Créez un compte en quelques clics ou connectez-vous pour continuer.
Les posteurs les plus actifs de la semaine
Aucun utilisateur

-40%
Le deal à ne pas rater :
-40% sur le Pack Gaming Mario PDP Manette filaire + Casque filaire ...
29.99 € 49.99 €
Voir le deal

Test de wilcox sour R

2 participants

Aller en bas

Test de wilcox sour R Empty Test de wilcox sour R

Message par Komine Ven 2 Juil 2010 - 12:28

Bonjour,
j'ai fait un test de wilcox en utilisant cette syntaxe suivante et ça a marché.
Code:
wilcox.test(x=tab$FMC[tab$Classe_ech==2],y=tab$FMC[tab$Classe_ech==1],alternative="greater")
Cependant je veux ajouter un deuxième facteur (tab$Site==] en plus du premier facteur qui est [tab$Classe_ech==]

Donc j'ai écrit la syntaxe de cette façon, mais cela ne marche pas:
Code:
wilcox.test(x=tab$FMC[tab$Classe_ech==2,[tab$Site==1]],y=tab$FMC[tab$Classe_ech==1,[tab$Site==1]],alternative="greater")

Je sollicite votre aide pour écrire cette syntaxe.
Merci d'avance

Komine

Nombre de messages : 78
Date d'inscription : 02/06/2010

Revenir en haut Aller en bas

Test de wilcox sour R Empty Re: Test de wilcox sour R

Message par Nik Ven 2 Juil 2010 - 12:38

Comme ça ça devrait marcher :

Code:
wilcox.test(x=tab$FMC[tab$Classe_ech==2 & tab$Site==1],y=tab$FMC[tab$Classe_ech==1 & tab$Site==1],alternative="greater")

nik

Nik

Nombre de messages : 1606
Date d'inscription : 23/05/2008

Revenir en haut Aller en bas

Test de wilcox sour R Empty Test de wilcox sour R

Message par Komine Ven 2 Juil 2010 - 14:06

Bonjour,
J'ai apporté des changement de variables dans ma table mais l'esprit de la question reste le même. Cependant l'opération ne parche toujours pas. Je vous montre un échantillon de mes donnée et le message d'erreur.
Merci encore
Code:
  Site Jour FMC_ligneux FMC_herb Diff_FMC EWT._ligneux EWT._herb Diff_EWT
1    1    1      366.67    23.06  343.61      244.44    56.22  188.22
2    1    2      135.56    39.04    96.52      118.52    57.33    61.19
3    1    3      110.03    53.96    56.07      103.70    59.89    43.81
4    1    4      89.67    38.10    51.57        83.33    43.56    39.77
5    1    5      71.79    80.28    -8.49      107.41    74.33    33.08
6    2    1      87.34    46.18    41.16      149.36    69.50    79.86
7    2    2      88.42    73.16    15.26      121.15    66.50    78

 
 wilcox.test(x=tab$FMC_ligneux[tab$Jour==(1:5)& tab$Site==1],(y=tab$FMC_herb[tab$jour==(1:5)& tab$Site==1],alternative="greater")
Erreur : ',' inattendu(e) dans "wilcox.test(x=tab$FMC_ligneux[tab$Jour==(1:5)& tab$Site==1],(y=tab$FMC_herb[tab$jour==(1:5)& tab$Site==1],"

Komine

Nombre de messages : 78
Date d'inscription : 02/06/2010

Revenir en haut Aller en bas

Test de wilcox sour R Empty Test de wilcox sour R

Message par Komine Ven 2 Juil 2010 - 14:31

Bonjour,
Finalement ça marche. Il y avait une parenthèse de plus. Mais dans le résultat, il ya un message d'avis, est ce un problème aussi? voir résultat:
Code:
 Wilcoxon rank sum test

data:  tab$FMC_ligneux[tab$Jour == (1:5) & tab$Site == 1] and tab$FMC_herb[tab$Jour == (1:5) & tab$Site == 1]
W = 24, p-value = 0.007937
alternative hypothesis: true location shift is greater than 0
Messages d'avis :
1: In tab$Jour == (1:5) :
  la taille d'un objet plus long n'est pas multiple de la taille d'un objet plus court
2: In tab$Jour == (1:5) :
  la taille d'un objet plus long n'est pas multiple de la taille d'un objet plus court
Merci encore

Komine

Nombre de messages : 78
Date d'inscription : 02/06/2010

Revenir en haut Aller en bas

Test de wilcox sour R Empty Re: Test de wilcox sour R

Message par Komine Dim 4 Juil 2010 - 11:32

Bonjour,
j'ai un problème pour faire le test de wilcoxon.
je montre d'abord un échantillon des mes données:
Code:
  Site Jour FMC_ligneux FMC_herb Diff_FMC EWT._ligneux EWT._herb Diff_EWT
1    1    1      366.67    23.06  343.61      244.44    56.22  188.22
2    1    2      135.56    39.04    96.52      118.52    57.33    61.19
3    1    3      110.03    53.96    56.07      103.70    59.89    43.81
4    1    4      89.67    38.10    51.57        83.33    43.56    39.77
5    1    5      71.79    80.28    -8.49      107.41    74.33    33.08
6    2    1      87.34    46.18    41.16      149.36    69.50    79.86
7    2    2      88.42    73.16    15.26      121.15    66.50    54.65
8    2    3      96.83    76.23    20.60      139.10    73.50    65.60
9    2    4      100.64    57.51    43.13      109.62    49.75    59.87
10    2    5      105.38    57.93    47.45      112.82    48.75    64.07
11    2    6      95.66    42.74    52.92      110.90    41.75    69.15
12    2    7      87.51    46.41    41.10        92.95    38.75    54.20
13    2    8      75.47    36.69    38.78        76.92    28.50    48.42
j'utilise le code suivant pour faire le test de wilcoxon:

Code:
wilcox.test(x=tab$FMC_ligneux[tab$Jour==(1:5)& tab$Site==1],y=tab$FMC_herb[tab$Jour==(1:5)& tab$Site==1],alternative="greater")
Le résultat est le suivant et je ne comprends pas le message d'avis?

Code:
        Wilcoxon rank sum test

data:  tab$FMC_ligneux[tab$Jour == (1:5) & tab$Site == 1] and tab$FMC_herb[tab$Jour == (1:5) & tab$Site == 1]
W = 24, p-value = 0.007937
alternative hypothesis: true location shift is greater than 0
Messages d'avis :
1: In tab$Jour == (1:5) :
  la taille d'un objet plus long n'est pas multiple de la taille d'un objet plus court
2: In tab$Jour == (1:5) :
  la taille d'un objet plus long n'est pas multiple de la taille d'un objet plus court

Je mets un 2ème code qui lui traite la 2ème partie de mes données:
Code:
 wilcox.test(x=tab$FMC_ligneux[tab$Jour==(6:13)& tab$Site==2],y=tab$FMC_herb[tab$Jour==(6:13)& tab$Site==2],alternative="greater")

Le résultat n'est pas bon:
Code:
Erreur dans wilcox.test.default(x = tab$FMC_ligneux[tab$Jour == (6:13) &  :
  pas assez d'observations 'x' (définies)
De plus : Messages d'avis :
1: In tab$Jour == (6:13) :
  la taille d'un objet plus long n'est pas multiple de la taille d'un objet plus court
2: In tab$Jour == (6:13) :
  la taille d'un objet plus long n'est pas multiple de la taille d'un objet plus court

Pourriez vous me dire où se trouve le problème sur sur le dernier résultat?
Merci d'avance

Komine

Nombre de messages : 78
Date d'inscription : 02/06/2010

Revenir en haut Aller en bas

Test de wilcox sour R Empty Re: Test de wilcox sour R

Message par Contenu sponsorisé


Contenu sponsorisé


Revenir en haut Aller en bas

Revenir en haut

- Sujets similaires

 
Permission de ce forum:
Vous ne pouvez pas répondre aux sujets dans ce forum