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