Mesures appariées
Un article de IMSP - Formation continue.
On trouve typiquement des mesures appariées lors de mesures séquentielles chez les mêmes individus dans différentes situations, ou encore dans les études cas-contrôle. Les calculs et les critères sont similaires à ceux effectués avec un échantillon unique, à la différence près qu'ils sont effectués sur la différence des deux mesures. Les valeurs sont ainsi compensées des variations intra-personnelles, et l'erreur standard sera inférieure à celle qu'on obtiendrait si les mesures n'étaient pas appariées.
Sommaire |
Intervalles de confiance
le calcul de la moyenne
et de l'écart-type s s'effectue sur la différence des n paires, l'erreur standard e.s. est
.
- Pour les grands échantillons (60 paires et plus) on utilise le z-score:
- Pour les petits échantillons (moins de 60 paires) on utilise la valeur de t avec n − 1 degrés de liberté:
Test t de Student pour mesures appariées
Selon la taille de l'échantillon on utilise soit un test z apparié ou un test t apparié.
étant la moyenne des différences et s l'écart-type des différences:
grand effectif
effectif réduit
R
Un exercice donné par l'Université d'Alberta:
mod <- c(16.85, 16.40, 17.21, 16.35, 16.52, 17.04, 16.96, 17.15, 16.59, 16.57) ; unmod <- c(17.50, 17.63, 18.25, 18.00, 17.86, 17.75, 18.22, 17.90, 17.96, 18.15) ; differences <- mod-unmod
differences summary(differences) ;
effectuons un test t de Student apparié:
> t.test(mod, unmod,paired=TRUE) ;
Paired t-test
data: mod and unmod
t = -10.2311, df = 9, p-value = 2.958e-06
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
-1.4140405 -0.9019595
sample estimates:
mean of the differences
-1.158
comparons avec un test t de Student simple sur les différences:
> t.test(differences) ;
One Sample t-test
data: differences
t = -10.2311, df = 9, p-value = 2.958e-06
alternative hypothesis: true mean is not equal to 0
95 percent confidence interval:
-1.4140405 -0.9019595
sample estimates:
mean of x
-1.158
enfin comparons avec un test non apparié sur les mêmes données:
> t.test(mod, unmod,paired=FALSE)
Welch Two Sample t-test
data: mod and unmod
t = -9.1094, df = 17.025, p-value = 5.894e-08
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
-1.4261741 -0.8898259
sample estimates:
mean of x mean of y
16.764 17.922
STATA
en construction
Pour continuer
|
|
