Enoncé :
Source : exemple 10-items.sta fourni avec Statistica 7.1
On a fait passer un questionnaire comportant 10 items à un échantillon de 10 sujets. On souhaite mesurer la cohérence du questionnaire à l'aide du coefficient alpha de Cronbach.
Les données se présentent sous la forme suivante :
|
|
Q1
|
Q2
|
Q3
|
Q4
|
Q5
|
...
|
|
S1
|
4
|
5
|
4
|
6
|
4
|
...
|
|
S2
|
4
|
3
|
4
|
5
|
5
|
...
|
|
S3
|
4
|
4
|
3
|
5
|
7
|
...
|
|
S4
|
4
|
3
|
6
|
5
|
4
|
...
|
|
S5
|
4
|
6
|
7
|
5
|
5
|
...
|
|
S6
|
5
|
6
|
6
|
8
|
4
|
...
|
|
...
|
...
|
...
|
...
|
...
|
...
|
...
|
Ces données sont enregistrées au format .csv d'Excel dans le fichier : dix-items.csv, avec la structure indiquée ci-dessus.
On peut les charger dans R à l'aide de la commande :
> dix_items <- read.csv2("dix-items.csv")
N.B. read.csv2 permet de lire des données au format "table", avec les spécifications par défaut suivantes : séparateur = ";", séparateur décimal = ",".
On utilise la commande cronbach, qui se trouve dans le package psy.
> library(psy)
> dix_items
X Q1 Q2 Q3 Q4 Q5 Q6 Q7 Q8 Q9 Q10
1 S1 4 5 4 6 4 7 4 6 6 3
2 S2 4 3 4 5 5 4 4 4 2 4
3 S3 4 4 3 5 7 5 3 5 5 3
...
> cronbach(dix_items[,2:11])
$sample.size
[1] 100
$number.of.items
[1] 10
$alpha
[1] 0.794313
Les valeurs retournées par la commande cronbach :
> unlist(cronbach(dix_items[,2:11]))
sample.size number.of.items alpha
100.000000 10.000000 0.794313
Pour calculer le coefficient en se limitant aux questions 1 à 4 et 7 à 10 :
> cronbach(dix_items[,c(2:5,8:11)])
$sample.size
[1] 100
$number.of.items
[1] 8
$alpha
[1] 0.854978