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

 

# Vi vil illustrere teorien p? sidene 478-482 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/h10/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 479 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 481 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)