Compute a compact NA-cleaned numeric summary: N, min, Q1, median, mean, Q3, max.
distrib <- function(x, k = 6) {
# NA removing
x <- x[is.na(x) == FALSE]
Nx <- length(x)
sum <- round(summary(x), k)
res <- c(Nx, sum)
names(res) <- c("N", names(sum))
return(res)
} # end of function distrib