R kode som illustrerer konfidensintervall for et populasjonsgjennomsnitt ved store utvalg (avsnitt 8.2)

 

 

# Vi vil illustrere teorien p? sidene 385-386 i l?reboka ved ? se p? f?dselsvekten til 50 gutter (valgt tilfeldig fra en st?rre populasjon).

# Vi vil bestemme et (tiln?rmet) 95% konfidensintervall for forventet vekt m for nyf?dte gutter.

 

# Vi leser inn de 50 f?dselsvektene:

fvekt.gutt=c(3590,3300,3020,3450,3480,3900,3670,3650,4190,3650, 3905,2720,3560,3600,3850,3610,3900,3660,3130,4040, 3450,3060,3760,3735,3440,3450,4020,3430,3400,1900, 3320,4250,3490,3600,3400,4200,3470,3900,3520,2660, 3290,2750,3480,2700,3750,3230,4100,4020,4300,3240)

 

# Vi betemmer nedre og ?vre grense for konfidensintervallet [jf formel (8.8) i l?reboka]

n=length(fvekt.gutt)

s=sd(fvekt.gutt)

l=mean(fvekt.gutt)-1.96*s/sqrt(n)

u=mean(fvekt.gutt)+1.96*s/sqrt(n)

print(c(l,u),4)