analyse croisée de toutes les variables d'un tableau

Version imprimable

17/02/2012, 12h17
CharlesHUS

analyse croisée de toutes les variables d'un tableau
Bonjour
J'ai un data.frame à 32 variables (en colonnes) et 125 individus en ligne.
J'aimerais effectuer un test croisé (du chi² en l'occurence) de toutes les variables entre elles.
J'arrive a croiser une variable avec toutes les autres avec cette instruction :

Code:

round((sapply(tableauanalysechideux[,2:32],function(x)chisq.test(tableauanalysechideux$SEXE,x)$p.value)),3)

ou celle-ci :
Code:

1 2 3 4 chisq.pvalue <- function(x) { chisq.test(tableauanalysechideux$SEXE,x)$p.value } sapply(tableauanalysechideux[,2:32],chisq.pvalue)
mais pas de toutes les variables entre-elles.
J'ai essayé ceci mais sans succès :
Code:

1 2 3 4 analysechideux <- function(x) chisq.test(table(ex$x,ex$y)) x<- tableauanalysechideux[,2:32] y<- tableauanalysechideux[,2:32] analysechideux(x,y)
Je pense que la solution est assez simple mais je stagne depuis un bout de temps.

merci d'avance de vos suggestions

Bonjour,

Code:

1
2
3
4
5
6
for (i in 1:(ncol(USJudgeRatings) - 1)) {
 for (j in (i+1):ncol(USJudgeRatings)) {
  print(paste("Chisq.test between : ", colnames(USJudgeRatings)[i], " and ", colnames(USJudgeRatings)[j], sep=""))
  print(chisq.test(USJudgeRatings[, i], USJudgeRatings[, j])$p.value)
 }
}

Un truc dans ce genre là, non?

HTH

Vincent

Merci pour ta réponse.
Je vais l'essayer pour voir.

Pour n'obtenir que la valeur p d'un test du chi² :

Code:

1
2
3
4
5
6
7
8
9
resultatschi²<-matrix(0,nrow=dim(variablesqualitatives)[2],ncol=dim(variablesqualitatives)[2])
colnames(resultatschi²)<-names(variablesqualitatives)
rownames(resultatschi²)<-names(variablesqualitatives)
for(i in 1:(dim(variablesqualitatives)[2])){
   for(j in i:(dim(variablesqualitatives)[2])){
       resultatschi²[i,j]<-chisq.test(variablesqualitatives[,i],variablesqualitatives[,j])$p.value
       }
   }
round((resultatschi²),4)

Pour n'obtenir que la valeur de p pour une régression linéaire :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
#construction du tableau
resultatsreglin<-matrix(0,nrow=dim(variablesquantitatives)[2],ncol=dim(variablesquantitatives)[2])
colnames(resultatsreglin)<-names(variablesquantitatives)
rownames(resultatsreglin)<-names(variablesquantitatives)
 
# boucle de programmation
for(i in 1:(dim(variablesquantitatives)[2])){
   for(j in i:(dim(variablesquantitatives)[2])){
       resultatsreglin[i,j]<-round((summary(lm(variablesquantitatives[,i]~variablesquantitatives[,j], data=variablesquantitatives))$coef[,"Pr(>|t|)",drop=F][2,]),4)
       }
   }
resultatsreglin

Je galère maintenant pour ajouter automatiquement la droite de régression sur l'ensemble de mes graphiques "plot".

Merci.

[EDIT]

Pour les graphiques j'ai trouvé ceci :

Code:

1
2
3
4
5
6
7
par(ask=T,mfrow=c(5,5))
for(i in 1:(dim(variablesquantitatives)[2])){
   for(j in i:(dim(variablesquantitatives)[2])){
       plot(variablesquantitatives[,i],variablesquantitatives[,j], ylab=names(variablesquantitatives)[j],xlab=names(variablesquantitatives)[i],col="blue")
       abline(lm(variablesquantitatives[,i]~variablesquantitatives[,j], data=variablesquantitatives),col="red")
}
   }