Bonjour tout le monde,
Je sollicite votre aide, car je suis confronté à une erreur que je ne comprends pas.
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
read.table("dureeft.txt",header=TRUE)
tous<-read.table("dureeft.txt",header=TRUE)
Ymvta<-tous$Ymvta
compresse<-tous$compresse
a1<-aov(Ymvta~compresse)
shapiro.test(a1$res)
friedman.test(Ymvta~compresse | nom_bebe, data=tous)
friedman.test.with.post.hoc(Ymvta~compresse | nom_bebe, data=tous)
Le code ci dessus fonctionne correctement.
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
23
24
25
26
27
 
library("pgirmess", lib.loc="~/R/win-library/3.3")
read.table("dureeft.txt",header=TRUE)
tous<-read.table("dureeft.txt",header=TRUE)
compresse=tous$compresse
banane=0
for (cpt in 3:ncol(tous)) {
  comportement=unlist(tous[cpt])
  ao<-aov(comportement~compresse)
  if (shapiro.test(ao$res)$p.value<0.05){
    if (friedman.test(comportement~compresse | nom_bebe, data=tous)$p.value<0.05){
      source("https://www.r-statistics.com/wp-content/uploads/2010/02/Friedman-Test-with-Post-Hoc.r.txt")
      banane=banane+1
      print(cpt)
      friedman.test.with.post.hoc(comportement~compresse | nom_bebe, data=tous)
    }
  }else{
    if (bartlett.test(comportement~compresse, data=tous)$p.value>0.05){
      truc=summary(ao)[[1]]
      if(truc[1,"Pr(>F)"]<0.05){
        print(summary(ao))
      }
    }
  }
 
}
print(banane)
En revanche celui ci ne fonctionne pas, je me retrouve avec l'erreur suivante: "Error in `[.data.frame`(data, , c(Y.name, X.name, block.name)) :
undefined columns selected"
L'erreur provient de la fonction post hoc cependant, c'est la même syntaxe que pour le test de friedman. Or le test de friedman fonctionne parfaitement mais pas le test post hoc.
Si l'un d'entre vous à des suggestions, je vous remercie par avance.

Mathieu