determin.qualiquali<- function(x,y,tabcontingence=FALSE) { # détermination du coeffcient de Cramer V, de la taille d'effet w et de la probabilité critique # pour le croisement de deux variable qualitatives # ********************************************************* # determin.qualiquali<- function(x,y) : # * Variables d'entrées : # * x,y : les deux vecteurs des données qualitatives # * tabcontingence : booléen faux par défaut ; si vrai, sortie de la table de contingence # * Variables de sortie : # * table.cont : éventuellement table de contingence # * chi2 : le chi au carré # * V : le coefficient de Cramer # * w : la taille d'effet # * pc: probabilité critique res <- chisq.test(x, y, correct = FALSE) chi2 <- as.numeric(res$statistic) n <- length(x) p <- length(levels(x)) q <- length(levels(y)) m <- min(p - 1, q - 1) V <- sqrt(chi2/(n * m)) w <- sqrt(chi2/n) if (tabcontingence) return(list(table.cont=res$observed,chi2=chi2,V=V,w=w,pc=res$p.value)) else return(list(chi2=chi2,V=V,w=w,pc=res$p.value)) }