Modérateur : Groupe des modérateurs
Code : Tout sélectionner
#l'exemple de la fiche br3 je pense qu'il y a une autre fiche plus détaillée
pomme<-matrix(0,nrow=5,ncol=4)
pomme[1,]<-c(527,604,606,533)
pomme[2,]<-c(633,600,650,567)
pomme[3,]<-c(642,708,662,504)
pomme[4,]<-c(623,550,562,667)
pomme[5,]<-c(377,408,500,333)
colnames(pomme) <- paste("bloc.",1:4,sep="")
rownames(pomme) <- paste("traitement.",1:5,sep="")
pomme
bloc.1 bloc.2 bloc.3 bloc.4
traitement.1 527 604 606 533
traitement.2 633 600 650 567
traitement.3 642 708 662 504
traitement.4 623 550 562 667
traitement.5 377 408 500 333
pomme.vec<-as.vector(pomme)
traitement<-as.factor(row(pomme))
traitement
[1] 1 2 3 4 5 1 2 3 4 5 1 2 3 4 5 1 2 3 4 5
Levels: 1 2 3 4 5
bloc<-as.factor(col(pomme))
bloc
[1] 1 1 1 1 1 2 2 2 2 2 3 3 3 3 3 4 4 4 4 4
Levels: 1 2 3 4
friedman.test(pomme.vec,traitement,bloc)
Friedman rank sum test
data: pomme.vec, traitement and bloc
Friedman chi-squared = 9.8, df = 4, p-value = 0.04393
# ce qui équivaut à mettre en ligne les blocs et en colonne les traitements :
pomme <- t(pomme)
pomme
traitement.1 traitement.2 traitement.3 traitement.4 traitement.5
bloc.1 527 633 642 623 377
bloc.2 604 600 708 550 408
bloc.3 606 650 662 562 500
bloc.4 533 567 504 667 333
friedman.test(pomme)
Friedman rank sum test
data: pomme
Friedman chi-squared = 9.8, df = 4, p-value = 0.04393
# pour ton cas à toi, il faut commencer par vérifier que pour une période donnée tu n'as qu'une seule valeur pour les aliments :
aliments <- read.table("friedman.txt",dec=",",h=T)
aliments
age periode dia conso
1 308 J15-J18 d2.5 0.13
2 308 J18-J22 d2.5 0.95
3 308 J15-J18 d3.5 0.27
4 308 J18-J22 d3.5 0.90
table(aliments$periode,aliments$dia)
d2.5 d3.5
J15-J18 1 1
J18-J22 1 1
# jusque la ça va, il te faut ensuite mettre en forme le tableau :
aliments <- reshape(aliments,direction="wide",idvar="periode",timevar="dia",drop="age")
rownames(aliments) <- aliments[,1]
aliments <- aliments[,-1]
aliments <- as.matrix(aliments)
aliments
conso.d2.5 conso.d3.5
J15-J18 0.13 0.27
J18-J22 0.95 0.90
# après tu n'as plus qu'a faire le test
friedman.test(aliments)
Code : Tout sélectionner
library(lattice)
bwplot(dia ~ conso | periode, data = aliment, layout = c(1, 8))
Retourner vers « Questions en cours »
Utilisateurs parcourant ce forum : Aucun utilisateur enregistré et 1 invité