R kode til eksempel 8.15

 

 

# Vi leser inn dataene om alkoholprosent for et utvalg av 16 ?lsorter (jf. eksempel 8.11 p? 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 finner konfidensintervallet for variansen til alkoholprosentene ved ? bruke resultatet p? side 403 i l?reboka:

n=length(prosent)

s2=var(prosent)

chi.975=qchisq(0.025,n-1)

chi.025=qchisq(0.975,n-1)

l=(n-1)*s2/chi.025

u=(n-1)*s2/chi.975

print(c(l,u),3)

 

# Vi finner konfidensintervallet for standardavviket:

print(sqrt(c(l,u)),3)