source("hdr.txt")
x=c(9,14,15,15,16,18,18,19,19,20,21,22,22,24,24,26,27,29,30,32)
n <- length(x)
nu <- n-1
xbar <- mean(x)
S <- sum((x-xbar)^2)
s <- sd(x)
tvalues <- ht(0.95,19)
cat("Interval for mu                                   ",
    (s/sqrt(n))*tvalues+xbar,"\n")
cat("Values of chi^2 corresponding to HDR for log chi^2",
    hchisq(0.95,nu),"\n")
cat("Interval for phi                                  ",
    sort(S/hchisq(0.95,nu)),"\n")
cat("Interval for phi assuming mean known              ",
    sort(S/hchisq(0.95,n)),"\n")