R kode som illustrerer konfidensintervall for et tilfeldig utvalg fra en eksponentialfordeling (avsnitt 8.1)

 

 

# Vi tar for oss eksempel 8.5 p? side 382-383 i l?reboka.

# Her ser en p? hvor lang tid (i minutter) det tar f?r en isolerende v?ske bryter sammen.

# Vi modellerer tidene som et tilfeldig utvalg fra eksponentialfordelingen med parameter l.

# Vi vil bestemme et 95% konfidensintervall for l og for forventningsverdien 1/l.

 

# Vi leser inn de 10 tidene:

tid=c(41.53,18.73,2.99,30.34,12.33,117.52,73.02,223.63,4.00,26.78)

 

# Vi betemmer nedre og ?vre grense for konfidensintervallet for l. (jf. eksempel 8.5):

n=length(tid)

x.025=qchisq(0.025,2*n)

x.975=qchisq(0.975,2*n)

l=x.025/(2*sum(tid))

u= x.975/(2*sum(tid))

print(c(l,u),4)

 

 

# Vi bestemmer nedre og ?vre grense for konfidensintervallet for forventningsverdien 1/l:

print(c(1/u,1/l),4)