Solutions des TPs
653
table(sexe) # Effectifs.
# Fr´ equences en pourcentages:
round(table(sexe)/length(sexe) * 100,2)
names(which.max(table(sexe))) # Mode de la variable.
barplot(sort(table(sexe),TRUE),col=c("pink","blue"),
pareto=TRUE)
9.5- yeuxcheveux <- table(yeux,cheveux) # Table de contingence.
9.6- fchev <- margin.table(yeuxcheveux,2)/592 # Profils colonnes.
fchev
9.7- # Nombre d’individus ayant les yeux bleux:
nbleus <- margin.table(yeuxcheveux,1)[1]
nbleus
9.8- round(fchev * nbleus)
9.9- margeX <- margin.table(yeuxcheveux,1)
tab.ind1 <- margeX% * %t(fchev)
round(tab.ind1)
9.10- fyeux <- margin.table(yeuxcheveux,1)/592 # Profils lignes.
fyeux
9.11- # Nombre d’individus ayant les cheveux blonds:
nblonds <- margin.table(yeuxcheveux,2)[1]
nblonds
9.12- round(fyeux * nblonds)
9.13- margeY <- margin.table(yeuxcheveux,2)
tab.ind2 <- fyeux% * %t(margeY)
round(tab.ind2)
9.14- all.equal(tab.ind1,tab.ind2)
Les deux tableaux sont identiques, ce qui est rassurant puisque s’int´ eresser
` a l’ind´ ependance des yeux et des cheveux est ´ equivalent `
a s’int´ eresser `
a
l’ind´ ependance des cheveux et des yeux.
9.15- (yeuxcheveux-tab.ind1)^2
9.16- tab.contr <- (yeuxcheveux-tab.ind1)^2/tab.ind1
tab.contr
9.17- khi2 <- sum(tab.contr)
khi2
Phi2 <- khi2/sum(yeuxcheveux)
Phi2