#####   R-code for making Figure 2.12

	dat  <- read.table("soil.dat",header=F)
	x    <- dat$V2

	postscript("fig212.ps",width=5,height=5,horizontal=F)

	fit1 <- density(x,bw =0.05)
	fit2 <- density(x,bw =0.25)
	
	plot(fit1$x,fit1$y,type="l",xlab=expression(SiO[2]),
             ylab="Density Estimates",mgp=c(2,1,0),cex=0.8,main="")
	lines(fit2$x,fit2$y,lty=2,cex=0.8)
	points(x,rep(0.05,length(x)),pch=".",cex=0.9)

	graphics.off()