library(mlbench)
 data(BostonHousing)
 s <- sort(BostonHousing[,13], index.return=T)
 x <- BostonHousing[s$ix,13]
 png("NWWhi.png")
 d <- -0.02*runif(length(x))-0.01
 plot(x, d, cex=0.5, ylim=c(-0.03, 0.12), xlab="Anteil der Unterschichtbevölkerung pro Bezirk (in %)", ylab="Gewichtsfunktionen W_hi(x)",  main="Boston Housing Daten")
 xk <- 30
 h  <- 2.1
 k <- dnorm((x-xk)/h)
 k <- k/sum(k)
 lines(x,k, col="red", lwd=2)
 xk <- 20
 k <- dnorm((x-xk)/h)
 k <- k/sum(k)
 lines(x,k, col="blue", lwd=2)
 xk <- 10
 k <- dnorm((x-xk)/h)
 k <- k/sum(k)
 lines(x,k, col="green", lwd=2)
 h  <- 0.7
 xk <- 30
 k <- dnorm((x-xk)/h)
 k <- k/sum(k)
 lines(x,k, col="red", lwd=2, lty=3)
 xk <- 20
 k <- dnorm((x-xk)/h)
 k <- k/sum(k)
 lines(x,k, col="blue", lwd=2, lty=3)
 xk <- 10
 k <- dnorm((x-xk)/h)
 k <- k/sum(k)
 lines(x,k, col="green", lwd=2, lty=3)
 legend("topleft", legend=c("W_hi(10), h=2,1", "W_hi(20), h=2,1", "W_hi(30), h=2,1", "W_hi(10), h=0,7", "W_hi(20), h=0,7", "W_hi(30), h=0,7"),  col=c("green", "blue", "red", "green", "blue", "red"), lwd=2, lty=c(1,1,1,3,3,3))
 dev.off()