Is there any function that divides a sample into N quantiles? For example, for N = 2, this would be the solution: x <- rnorm(100) m <- median(x) q <- ifelse(x <= median, 1, 2) Alberto Monteiro