Initiation ` a l'analyse statistique des donn' ees exp ...

Description

I
n
i
t
i
a
t
i
o
n
`
a
l
’
a
n
a
l
y
s
e
s
t
a
t
i
s
t
i
q
u
e
d
e
s
d
o
n
n
´
e
e
s
e
x
p
´
e
r
i
m
e
n
t
a
l
e
s
(
C
o
m
p
a
r
a
i
s
o
n
s
d
e
m
o
y
e
n
n
e
s
)
C
h
r
i
s
t
o
p
h
e
P
a
l
l
i
e
r
e
t
C
h
r
i
s
t
o
p
h
e
L
a
l
a
n
n
e
c
h
r
i
s
t
o
p
h
e
.
p
a
l
l
i
e
r
@
m
4
x
.
o
r
g
/
c
h
r
i
s
t
o
p
h
e
.
l
a
l
a
n
n
e
@
g
m
x
.
n
e
t
M
a
s
t
`
e
r
e
d
e
S
c
i
e
n
c
e
s
C
o
g
n
i
t
i
v
e
s
,
E
H
E
S
S
–
P
a
r
i
s
5
–
E
N
S
S
t
a
t
s
a
v
e
c
R
–
p
.
1R
a
p
p
e
l
s
b
a
s
e
d
e
l
’
i
n
f
´
e
r
e
n
c
e
p
o
p
u
l
a
t
i
o
n
,
´
e
c
h
a
n
t
i
l
l
o
n
a
l
´
e
a
t
o
i
r
e
,
d
i
s
t
r
i
b
u
t
i
o
n
d
’
´
e
c
h
a
n
t
i
l
l
o
n
n
a
g
e
t
e
s
t
s
d
’
h
y
p
o
t
h
`
e
s
e
s
u
r
u
n
e
s
t
a
t
i
s
t
i
q
u
e
,
s
t
a
t
i
s
t
i
q
u
e
d
e
t
e
s
t
f
o
r
m
u
l
a
t
i
o
n
d
e
H
0
r
i
s
q
u
e
s
d
’
e
r
r
e
u
r
i
n
t
e
r
v
a
l
l
e
s
d
e
c
o
n
fi
a
n
c
e
H
v
r
a
i
e
H
f
a
u
s
s
e
0
0
R
e
j
e
t
d
e
H
α
1
−
α
0
A
c
c
e
p
t
a
t
i
o
n
d
e
H
1
−
β
β
0
c
o
m
p
a
r
a
i
s
o
n
s
d
e
m
o
y
e
n
n
e
s
:
t
e
s
t
t
e
t
A
N
O
V
A
`
a
u
n
f
a
c
t
e
u
r
´
e
c
h
a
n
t
i
l
l
o
n
s
i
n
d
´
e
p
e
n
d
a
n
t
s
o
u
a
p
p
a
r
i
´
e
s
o
b
s
e
r
v
a
t
i
o
n
u
n
i
q
u
e
o
u
r
´
e
p
l
i
q
u
e
s
S
t
a
t
s
a
v
e
c
R
–
p
.
2P
u
i
s
s
a
n
c
e
d
’
u
n
t
e
s
t
S
t
a
t
s
a
v
e
c
R
–
p
.
3C
o
m
p
a
r
a
i
s
o
n
`
a
u
n
e
m
o
y
e
n
n
e
d
e
r
´
e
f
´
e
r
e
n
c
e
1
´
e
c
h
a
n
t
i
l
l
o
n
a
l
´
e
a
t
o
i
r
e
d
e
t
a
i
l
l
e
n
H
:
μ
=
μ
0
0
e
x
e
m
p
l
e
:
x
<
-
r
n
o
r
m
(
2
0
,
m
e
a
n
=
0
.
7
5
)
p
a
r
(
m
f
c
o
l
=
c
(
1
...

Sujets

• Ittoqqortoormiit
• Apparitions mariales de Me?ugorje
• Banque d'Estonie
• GLUT
• Rite de passage
• Design transport

Informations

• Publié par: Sewyag
• Nombre de lectures: 97
• Langue: Español

Extrait

Initiation`al’analysestatistiquedesdonn´ees expe´rimentales (Comparaisons de moyennes)

Christophe Pallier et Christophe Lalanne [e-mail] / christophe.lalanne@gmx.n et

Mast`eredeSciencesCognitives,EHESS–Paris5–ENS

Stats avec R – p. 1

Rappels

basedel’infe´rence population,e´chantillonal´eatoire,distributiond’e´chantillonnage testsd’hypothe`sesurunestatistique,statistiquedetest formulation deH 0 risques d’erreur intervalles de confiance

Rejet deH 0 Acceptation deH 0

H 0 vraie α 1−β

comparaisons de moyennes : testtetANOVA`aunfacteur e´chantillonsinde´pendantsouapparie´s observationuniqueoure´pliques

H 0 fausse 1−α β

Stats avec R – p. 2

Puissance d’un test

Stats avec R – p. 3

Comparaisona`unemoyennedere´fe´rence

1e´chantillonal´eatoiredetaillen H 0 :µ=µ 0 exemple : x < rnorm(20,mean=0.75) par(mfcol=c(1,2)) stripchart(x,method="jitter",vert=T,pch=19) abline(h=0,lty=3) abline(h=mean(x),col="red") boxplot(x) abline(h=0,lty=3) t.test(x,mu=0,alt="two.sided")

comparableautestdetypicalit´e:notiond’limaet´exˆetrd’une statistique parrapport`aunedistributionder´ef´erence(ici,N(0; 1))

Stats avec R – p. 4

Comparaison de deux moyennes (ind.) 1/2

2e´chantillonsale´atoiresi.i.d.detaillen H 0 :µ 1 =µ 2 (≡µ 1 −µ 2 = 0) exemple : nsuj < 60 effectsize < 0.5 subject < gl(nsuj,1) group < gl(2,nsuj/2) scores < 1+rnorm(nsuj)+((0:(nsuj1))%/%(nsuj/2))*effectsize par(mfcol=c(1,3)) stripchart(scores~group,vert=T,method=’jitter’) plot.design(scores~subject+group) boxplot(scores~group) by(scores,group,mean) by(scores,group,summary) t.test(scores~group)

Stats avec R – p. 5

Comparaison de deux moyennes (ind.) 2/2

alternativenonparame´trique:testdeWilcoxon(rangs) wilcox.test(scores~group)

N.B. ve´rificationdeshypoth`eses:(i)normalite´desre´sidus(e.g.qqnorm(), ks.test()it´e()(ietoh)icsomade´citsbartlett.test()) pard´efaut,Reffectueletestmodifie´deWelch(aussipuissantqueletest classiquedeStudent),sionnesp´ecifiepasl’optionvar.equal=T

Stats avec R – p. 6

Comparaison de deux moyennes (app.) 1/2

2e´chantillonsal´eatoiresapparie´sdetaillen meˆmehypoth`esenulle exemple : nsuj < 20 subject < gl(nsuj,2) cond < gl(2,1,2*nsuj) effectsize < 2 scores < (as.numeric(cond)1)*effectsize+rnorm(2*nsuj) data < data.frame(subject,cond,scores) par(mfcol=c(1,2)) plot.design(data) interaction.plot(cond,subject,scores) t.test(scores~cond,paired=T)

Stats avec R – p. 7

Comparaison de deux moyennes (app.) 2/2

on tient compte de la structure d’appariement dans l’estimation de la variance alternativenonparam´etrique:testdeMannWhitneyWilcoxon(rangs) wilcox.test()avec l’optionpaired=T

Stats avec R – p. 8

Comparaisondekenoym3/1sen

klitnahce´resiatoial´elonseialld.tei.d.n k H 0 :µ 1 =µ 2 =. . .=µ k vs.H 1 :au moins une pairedmeyoneendsffi`ere exemple : ns < c(15,28,10,20,35) n < length(ns) group < factor(rep(1:n,ns),labels=paste("g",1:n,sep="")) data < rnorm(length(group),mean = 100 + (as.numeric(group)2)^2) par(mfrow=c(2,2)) plot(data~group) stripchart(data~group,method="jitter",vertical=T) plot.design(data~group) l < aov(data~group) summary(l)

Stats avec R – p. 9

Comparaisondeknney/2sem3o

poursuite de l’analyse : quelles sont les paires de moyennes significativementdiff´erentes(comparaisonsmultiples,posthoc) (hsd < TukeyHSD(l, "group", ordered = FALSE)) plot(hsd)

taille des effets : model.tables(l,se=T) summary.lm(l)

Stats avec R – p. 10

Comparaisondekoymneen/3s3

alternativenonparam´etrique:ANOVAdeKruskalWallis(rangs),avec la fonctionkruskal.test() N.B. ve´rificationdeshypothe`ses:(i)normalit´edesr´esidus(e.g.qqnorm(), ks.test()(e´ticitsade´cosom)hiit()ebartlett.test())

Stats avec R – p. 11