moda <- function(x) {
t <- table(x)
return(as.numeric(names(t)[t == max(t)]))
}
# Se generan 500 números enteros aleatorios con distribución normal (media=100, dt=2) y se calcula su moda
set.seed(10)
x <- trunc(rnorm(500,100,2))
print(moda(x))
## [1] 99