1with(tephra, {
2logit <- log(Al2O3/(100-Al2O3))
3par(mfrow=c(1,2))
4qqnorm(logit)
5qqline(logit)
6cat("ISE statistic:", nise(logit),"\n")
7sm <- sm.density(logit)
8y  <- sm$eval.points
9sd <- sqrt(hnorm(logit)^2 + var(logit))
10lines(y, dnorm(y, mean(logit), sd), lty = 3)
11par(mfrow=c(1,1))
12})
13
14