Los datos

datos <- c(17, 18, 20, 21, 20, 16, 15, 15, 19, 20, 22, 21, 20, 16, 15, 17, 20, 19, 15, 16, 20, 20, 19, 18, 17, 15, 16, 21, 16, 20, 22, 21, 20, 19, 20, 21, 20, 21, 20, 19, 18, 17, 17, 19, 20, 22, 15, 17, 20, 16)

# ordenados

sort(datos)
##  [1] 15 15 15 15 15 15 16 16 16 16 16 16 17 17 17 17 17 17 18 18 18 19 19
## [24] 19 19 19 19 20 20 20 20 20 20 20 20 20 20 20 20 20 20 21 21 21 21 21
## [47] 21 22 22 22
n <- length(datos) # cuantos son

n
## [1] 50

Frecuencias

frecuencia <- table(datos)
frecuencia
## datos
## 15 16 17 18 19 20 21 22 
##  6  6  6  3  6 14  6  3

Tabla de frecuencias

tabla <- transform(frecuencia)


tabla <- cbind(tabla, "Freq. Rel" = tabla$Freq / n, "Freq. %" = tabla$Freq / n * 100)

tabla
##   datos Freq Freq. Rel Freq. %
## 1    15    6      0.12      12
## 2    16    6      0.12      12
## 3    17    6      0.12      12
## 4    18    3      0.06       6
## 5    19    6      0.12      12
## 6    20   14      0.28      28
## 7    21    6      0.12      12
## 8    22    3      0.06       6

media mediana y moda

media <- mean(datos)

media
## [1] 18.56
mediana <- median(datos)

mediana
## [1] 19
moda <- max(frecuencia)

moda
## [1] 14