# dati
library(carData)

# MEDIANA
median(SLID$age)

# mediana per variabile con dati mancanti
median(SLID$education, na.rm = TRUE)

# MEDIA ARITMETICA
mean(SLID$age)

# media per variabile con dati mancanti
mean(SLID$education, na.rm = TRUE)

# trimmed mean
x <- c(2, 3, 3, 5, 5, 6, 6, 7, 8, 60)
mean(x)

mean(x, trim = 0.1)

# MEDIA QUADRATICA
sqrt(mean(SLID$age^2))

# MEDIA GEOMETRICA

n = length(SLID$age)
prod(SLID$age)^(1/n)

# oppure
prod(SLID$age)^(1/length(SLID$age))

# media geometrica con valori mancanti

# numero di casi validi
n = sum(!is.na(SLID$education))

# oppure
n = nval(SLID$education)    # pacchetto LabRS

# calcolo della media
prod(SLID$education, na.rm = TRUE)^(1/n)

# MEDIA ARMONICA
1/mean(1/SLID$age)