library(KernSmooth) set.seed<-1 x<-c(rnorm(20,sd=0.8,mean=0),rnorm(15,sd=1,mean=4),rnorm(15,sd=0.3,mean=-3)) n<-20+15+15 xx<-(-500:600)/100 y<-(20*dnorm(xx,sd=0.8,mean=0)+15*dnorm(xx,sd=1,mean=4)+15*dnorm(xx,sd=0.3,mean=-3))/n plot(xx,y,type="l") points(x,x*0,col=2,pch=4) par(mfrow=c(2,2)) plot(bkde(x,bandwidth=0.02),type="l") points(x,x*0,col=2,pch=4) plot(bkde(x,bandwidth=0.1),type="l") points(x,x*0,col=2,pch=4) plot(bkde(x,bandwidth=0.5),type="l") points(x,x*0,col=2,pch=4) plot(bkde(x,bandwidth=2),type="l") points(x,x*0,col=2,pch=4)