# NOT RUN {
y <- rnbinom(100, 10, 0.5)
classBreaks(y, 10)
classBreaks(y, 10, type="quantile")
par(mfrow=c(2,2))
d <- density(y)
plot(d, type="n", main="Equal Area breaks")
polygon(d, col="cyan")
abline(v=classBreaks(y, 10))
plot(d, type="n", main="Quantile breaks")
polygon(d, col="cyan")
abline(v=classBreaks(y, 10, type="quantile"))
plot(d, type="n", main="Robust Standard Deviation breaks")
polygon(d, col="cyan")
abline(v=classBreaks(y, 10, type="std"))
plot(d, type="n", main="Geometric interval breaks")
polygon(d, col="cyan")
abline(v=classBreaks(y, 10, type="geometric"))
( y.breaks <- classBreaks(y, 10) )
cut(y, y.breaks, include.lowest = TRUE, labels = 1:10)
# }
Run the code above in your browser using DataLab