inst/scripts/te_hcvsj.q

In sm: Smoothing Methods for Nonparametric Regression and Density Estimation

with(tephra, {
logit <- log(Al2O3/(100-Al2O3))
par(mfrow=c(1,2))
h.cv <- hcv(logit, display = "line", ngrid = 32)
n  <- length(logit)
sd <- sqrt(var(logit))
h  <- seq(0.003, 0.054, length=32)
lines(h, nmise(sd, n, h) - 5.5, lty = 3)
sm.density(logit, h.cv)
sm.density(logit, lty = 3, add = TRUE)
par(mfrow=c(1,1))
})