Beschreibung

The R program to create the PNG file:

library(np)
library(mlbench)
data(BostonHousing)
s <- sort(BostonHousing[,13], index.return=T)
x <- BostonHousing[s$ix,13]
y <- BostonHousing[s$ix,14]
h   <- (4:30)/20
mse <- ase <- 0*h
for (i in 1:length(x)) {
  for (j in 1:length(h)) {
    nwi    <- npreg(y[-i]~x[-i], bws=h[j], exdat=x[i])
    ase[j] <- ase[j] + (fitted(nwi)-y[i])^2
  }
}
for (j in 1:length(h)) {
  nwi    <- npreg(y~x, bws=h[j])
  mse[j] <- mse[j] + sum((fitted(nwi)-y)^2)
}
ase <- ase / length(x)
mse <- mse / length(x)
png("CV.png")
plot(h, ase, t="l", xlab="Bandweite h", ylab="Mittlerer Quadratischer Fehler", main="Boston Housing Daten", lwd="2", ylim=c(23,30), col="blue")
lines(h, mse, col="red", lwd=2)
legend("bottomright", legend=c("Leave-one-out Kreuzvalidierung", "Resubstitutions Schätzung"), col=c("blue", "red"), lwd=2, lty=1)
abline(v=0.7)
dev.off()

 

Dieses Chart wurde mit R erstellt.

Lizenz

Ich, der Urheberrechtsinhaber dieses Werkes, veröffentliche es als gemeinfrei. Dies gilt weltweit. In manchen Staaten könnte dies rechtlich nicht möglich sein. Sofern dies der Fall ist: Ich gewähre jedem das bedingungslose Recht, dieses Werk für jedweden Zweck zu nutzen, es sei denn, Bedingungen sind gesetzlich erforderlich.