R kode til eksempel 9.9

 

 

# Vi leser inn dataene:

maksvekt=c(25.8,36.6,26.3,21.8,27.2)

 

# Vi vil teste H₀ : m =25 mot H_a : m > 25

# Vi beregner testobservatoren ved ? bruke resultatet p? side 436 i l?reboka:

n=length(maksvekt)

s=sd(maksvekt)

t.stat=(mean(maksvekt)-25)/(s/sqrt(n))

 

# Vi utf?rer testingen ved kommandoen t.test :

t.test(maksvekt,alternative="greater",mu=25)