Test du chi2

Test du Chi2 :

Sur un vecteur : vect <- c(8, 12, 17); chisq.test(vect) : teste le vecteur contre l'hypothèse nulle d'équiprobabilité.
Sur un vecteur, mais contre une hypothèse autre que celle de l'équiprobabilité : vect <- c(8, 12, 17); chisq.test(vect, p = c(0.8, 0.1, 0.1)) (la somme des probabilités doit faire 1).
Sur une matrice : mat <- matrix(c(70, 100, 150, 350), nrow = 2); chisq.test(mat) : considère la matrice comme un tableau de contingence et renvoie un objet de type htest contenant toutes les informations.
On peut aussi utiliser le test sur une table : tab <- as.table(matrix(c(70, 100, 150, 350), nrow = 2); chisq.test(tab)
Pour une matrice mat, on peut calculer la matrice attendue étant données les fréquences marginales par : mat2 <- rowSums(mat) %o% colSums(mat) / sum(mat) et donc la statistique du Chi2 sans correction est : sum((mat - mat2) ^ 2 / mat2)
Accès aux différents champs de l'objet si tes <- chisq.test(mat) :
- tes$statistic : la statistique du Chi2.
- tes$parameter : le nombre de degrés de libertés.
- tes$p.value : la p-value.
- tes$observed : la matrice observée de départ.
- tes$expected : la matrice attendue sous l'hypothèse nulle d'absence de biais.

programmer en R, tutoriel R, graphes en R