Preface to First Edition - lib

152 DENSITY ESTIMATIONR> opar rx d1 d2 f hist(x, probability = TRUE, xlab = "Waiting times (in min.)",+ border = "gray", xlim = range(rx), ylim = c(0, 0.06),+ main = "")R> lines(rx, f, lwd = 2)R> lines(rx, dnorm(rx, mean = mean(x), sd = sd(x)), lty = 2,+ lwd = 2)R> legend(50, 0.06, lty = 1:2, bty = "n",+ legend = c("Fitted two-component mixture density",+ "Fitted single normal density"))Density0.00 0.01 0.02 0.03 0.04 0.05 0.06Fitted two−component mixture densityFitted single normal density40 50 60 70 80 90 100 110Waiting times (in min.)Figure 8.7Fitted normal density and two-component normal mixture for geysereruption data.© 2010 by Taylor and Francis Group, LLC