statmode <- function(x) { names(which.max(table(x))) } x <- c("a", "b", "b", "c", "c", "c", "d", "d", "e") y <- c(1, 1, 2, 2, 2, 3, 3, 3, 3, 4, 4, 4, 5, 5) z <- c(1.1, 1.1, 2.2, 2.3, 2.4, 3.1, 3.2, 3.4, 3.5, 4.8, 4.8, 4.8, 5.6, 5.6) (mode <- statmode(x)) (mode <- statmode(y)) (mode <- statmode(z))