R kode til eksempel 8.11

 

 

# Vi leser inn dataene om alkoholprosent for et utvalg av 16 ?lsorter (jf. side 395 i D&B):

prosent=c(4.68,4.13,4.80,4.63,5.08,5.79,6.29, 6.79,4.93,4.25,5.70,4.74,5.88,6.77,6.04,4.95)

 

# Vi lager et normalfordelingsplott (jf side 396):

qqnorm(prosent)

 

# Vi finner konfidensintervallet for forventet alkoholprosent ved ? bruke formel (8.15) i l?reboka:

n=length(prosent)

m=mean(prosent)

s=sd(prosent)

t.025=qt(0.975,n-1)

l=m-t.025*s/sqrt(n)

u=m+t.025*s/sqrt(n)

print(c(l,u),3)

 

# Vi finner konfidensintervallet (og mere til) ved ? bruke kommandoen t.test:

t.test(prosent)