library(DescTools)
library(moments) # skewness(...)

a <- c(12, 20, 17, 10, 21, 9, 7, 7, 21, 9, 19, 19, 5, 16, 11, 16, 9, 21, 4, 20, 15, 16, 3, 6, 5, 15, 9, 18)

mean(a) # Media aritmetica
Mode(a) # Moda
median(a) # Mediana
quantile(a, 0.28) # Valor en el percentil 28
tamaño=length(a) # Dimension del vector
# dim(d) # Dimensión en una matrix

var(a) # Cuasi-varianza
((length(a)-1)/length(a))*var(a) # Varianza
sd(a) # Cuasi-desviacion estandar/desviacion tipica
sqrt((length(a)-1)/length(a) * sd(a)) # Desviacion estandar/desviacion tipica
desvtipic <- sqrt((length(a)-1)/length(a) * sd(a)) # Igual pero copiandolo en variable
desvtipic/mean(a) # Coef variacion pearson
kurtosis(a) # Curtosis(leptocurtica>0,platicurtica<0,mesocurtica=0)
skewness(a) # Coef asimetria(no estoy del todo seguro)

is.numeric(a) # Check numeric
hist(a) # Histograma sobre a

Freq(a) # Tab frecuencias continua