library(densEstBayes) ; data(OldFaithful2011)
hist(OldFaithful2011,breaks = 25,col = "gold",
xlab = "time interval between geyser eruptions (minutes)")
dest <- densEstBayes(OldFaithful2011,method="SMFVB")
plot(dest,xlab = "time interval between geyser eruptions (minutes)")
rug(OldFaithful2011,col = "dodgerblue",quiet = TRUE)
Run the code above in your browser using DataLab