x <- c(2,4,2,5,8,9,4,5,2,1,4,4,6,5,2,8,5)
k <- 1 + 3.322 * log10(length(x))
table(cut(x, k))plot(table(cut(x, k)), type = "l")barplot(table(x))barplot(table(cut(x, k)))