R kode som illustrerer test og konfidensintervall for sammenligning av to populasjonsgjennomsnitt ved store utvalg (avsnitt 10.1)

 

# Vi vil illustrere teorien p� sidene 490-495 i l�reboka ved � se p� vekten til et tilfeldig utvalg av 2116 nyf�dte gutter og 1950 nyf�dte jenter.

 

# Vi vil (i) teste om det er forskjell p� forventet f�dselsvekt for gutter og jenter og (ii) bestemme et (tiln�rmet) 95% konfidensintervall for forskjellen i forventet f�dselsvekt.

 

# Vi leser inn datafilen som inneholder f�dselsvektene (og mere til):

fvekt=read.table("http://www.uio.no/studier/emner/matnat/math/STK1110/h12/fvekt.txt",header=T)

 

# Vi trekker ut f�dselsvektene for gutter og jenter:

fvekt.gutt=fvekt$vekt[fvekt$kjonn==1]

fvekt.jente=fvekt$vekt[fvekt$kjonn==2]

 

# Vi lager boxplott av f�dselsvektene:

boxplot(fvekt.gutt,fvekt.jente,names=c("Gutt","Jente"))

 

# Vi betemmer testobservatoren  [jf side 491 i l�reboka]

m=length(fvekt.gutt)

n=length(fvekt.jente)

s1=sd(fvekt.gutt)

s2=sd(fvekt.jente)

z=(mean(fvekt.gutt)-mean(fvekt.jente))/sqrt(s1^2/m+s2^2/n)

 

# Vi betemmer 95% konfidensintervall for forskjellen i forventet f�dselsvekt mellom gutter og jenter [jf side 493 i l�reboka]

l=mean(fvekt.gutt)-mean(fvekt.jente)-1.96*sqrt(s1^2/m+s2^2/n)

u=mean(fvekt.gutt)-mean(fvekt.jente)+1.96*sqrt(s1^2/m+s2^2/n)

print(c(l,u),2)