# The Chondrites data: Table 10.1
#
DF = data.frame(
    laboratory=as.factor( rep( c( 'A', 'B', 'C', 'D', 'E', 'F' ), 3 ) ),
    month=as.factor( c( rep( '1', 6 ), rep( '1', 6 ), rep( '1', 6 ) ) ),
    value=c(
        221.1, 208.8, 211.1, 208.3, 221.1, 224.2,
        224.2, 206.9, 198.4, 214.1, 208.8, 206.9,
        217.8, 205.9, 213.0, 209.1, 211.1, 198.4
    )
)

# Apply density estimation with a Gaussing kernal and several different bandwidths:
#
plot( density( DF$value, bw=5.16, adjust=2, kernel="gaussian" ) ) # one mode
plot( density( DF$value, bw=5.16, adjust=0.75, kernel="gaussian" ) ) # two modes
plot( density( DF$value, bw=5.16, adjust=0.5, kernel="gaussian" ) ) # three modes