Lezione 7

14 Pages·2011·0.4 MB·Italian

LaboratoriodiST1-Lezione7
AntoniettadiSalvatore
DipartimentodiMatematica
Universita´ degliStudiRomaTre
Outline
(cid:73) Lacostruzionedelteststatistico
(cid:73) errorediprimaspecieαesecondaspecieβ
(cid:73) Esercizisuitest
Relazionetraαeβ
Aparita´ din(numerosita´ campionaria),larelazionechesussistetraαeβe´
generalmenteesprimibilemedianteunafunzionedecrescente.
Consideriamodiaverestrattouncamionedin=50osservazionidaunaX(cid:118)N(µ,36)
evogliamoverificareilseguentesistemadiipotesi
H0:µ=0
H1:µ=3
v=sqrt(36/50)
curve(dnorm(x,0,v),-5,7, axes = TRUE, ylim =c(0,.5))
lines(c(0,0),c(0.45,0),lty=2)
curve(dnorm(x,3,v),-5,7, ylim =c(0,.5), add=T)
lines(c(3,3),c(0.45,0),lty=2)
disegnamoα
z=qnorm(0.95)
zv=z*v
x1=seq(zv,4,0.05)
y1=dnorm(x1,0,v)
polygon(c(x1,zv),c(y1,0), col=2)
eoradisegnamoβ
x2=seq(-1,zv,0.05)
y2=dnorm(x2,3,v)
polygon(c(x2,zv),c(y2,0), col=3)
legend(’topleft’, c(’alpha’,’beta’), col=c(2,3), pch=19)
5
0. alpha
beta
4
0.
v) 3
0, 0.
x,
m(
or 2
dn 0.
1
0.
0
0.
-4 -2 0 2 4 6
x
alpha=1-pnorm(x1,0,v)
beta=pnorm(x1,3,v)
plot(alpha,beta,type=’l’)
8
0.
6
0.
a
et
b 4
0.
2
0.
0
0.
0.00 0.01 0.02 0.03 0.04 0.05
alpha
Potenzadeltestγ
Visualizziamolapotenzadeltestperl’esercizioprecedente
curve(dnorm(x,0,v),-5,7, axes = TRUE, ylim =c(0,.5))
lines(c(0,0),c(0.45,0),lty=2)
curve(dnorm(x,3,v),-5,7, ylim =c(0,.5), add=T)
lines(c(3,3),c(0.45,0),lty=2)
z=qnorm(0.95)
Sappiamochelapotenzadeltesttestγe´ lapobabilita´ dirifiutareH0quandoquestae´
falsa,quindisihache
γ=1−β
disegnamoβ
x2=seq(-1,zv,0.05)
y2=dnorm(x2,3,v)
polygon(c(x2,zv),c(y2,0), col=3)
disegnogamma
gamma=pnorm(zv,3,v)
xg=seq(zv,7,0.05)
yg=dnorm(xg,3,v)
polygon(c(xg,zv),c(yg,0), col=4)
legend(’topleft’, c(’beta’,’gamma’), col=c(3,4), pch=19)
5
0. beta
gamma
4
0.
0, v) 0.3
x,
m(
dnor 0.2
1
0.
0
0.
-4 -2 0 2 4 6
x
Intervallidiconfidenzaeverificadelleipotesi
Esercizio1.
Ladegenzaospedalieraingiorniperiltrattamentodiunacertamalattiae´ stata
controllatasuuncampionedi10pazientiottenendoiseguentidati:
5.4,11.1,8.2,7.5,4.7,9.4,10.3,2.7,3.3,7.7.
Fissatounlivellodisignificativita´ del5%:
1. testaresulcampionel’ipotesichelamediasia0control’alternativachenonlo
siaetrovarel’intervallodiconfidenzaperlamediaµ.
2. testaresulcampionel’ipotesichelamediasia7control’alternativachenonlo
sia
3. calcolarelapotenzadeltestsel’ipotesialternativae´ media=8
4. calcolarelapotenzadeltestsel’ipotesialternativae´ media=18
5. cosaaccadeseilcampionee´ di50elementiconsd=2.1?
Dettodegilcampione
1. t.test(deg,alternative=’two.sided’)
2. t.test(deg,mu=7,alternative=’two.sided’)
3. power.t.test(n=10,sd=2.9,0.05,delta=1)
4. power.t.test(n=10,sd=2.9,0.05,delta=10)
5. power.t.test(50,sd=2.1,0.05,delta=1)
Esercizio2.
Inunostudioperaccertarel’efficaciadell’eserciziofisicoperdimagriresonostati
riportatiiseguentirisultati,interminidiriduzionedipesoinKg,relativiadungruppodi
9personeacuierastatoassegnatouncertoprogrammadieserciziperunmese:
5.9,3.2,-0.5,2.3,1.4,0.9,-0.5,0,3.6
assumendochelariduzionedipesoseguaunadistribuzionenormale
1. siverifichil’efficaciadelprogrammaallivellodell’1%
2. sicalcolilapotenzadeltestconalternativamedia=1
3. cosaaccadeallapotenzadeltestselasde´ 5?
Impostoilproblemacomemedia=0vsmedia>0.
Dettofisilcampione
t.test(fis, mu = 0, alternative = ’two.sided’)
sd(fis)
power.t.test(n=9,sd=2.152, 0.05, delta=1)
power.t.test(n=9,sd=5, 0.05, delta=1)
Esercizio3.
Ilphdell’acquachefuoriescedaunimpiantodidepurazionee´ specificatoallivellodi
7.0.Uncampionediacquaselezionatodaunimpiantoe´
5.2,5.7,5.8,6.5,6.1,5.9,6.9,6.0,6.7,6.7
C’e´ ragionedidubitarecircalaspecificazionedellivellodiphperl’impiantoallivello
del5%?