# rep
table(rep(Istat16$Eta.classi, times = Istat16$Totale))
summary(rep(Istat16$Età, times = Istat16$Totale))

# weighted.mean
weighted.mean(Istat16$Età, Istat16$Totale, na.rm = TRUE)

# xtabs
xtabs(Totale ~ Eta.classi, data = Istat16)
xtabs(Totale ~ Eta.classi + Rip.geog, data = Istat16)


#Hmisc
require(Hmisc)
describe(Istat16$Rip.geog, weights = Istat16$Totale)
describe(Istat16$Età, weights = Istat16$Totale)
wtd.var(Istat16$Età, weights = Istat16$Totale, na.rm = TRUE)
wtd.quantile(Istat16$Età, weights = Istat16$Totale, na.rm = TRUE)

# weights
require(weights)
stdz(Istat16$Età, weight = Istat16$Totale)
wpct(Istat16$Rip.geog, weight = Istat16$Totale, na.rm=TRUE)
wtd.hist(Istat16$Età, 
         weight = Istat16$Totale/1000,
         main = "Età al 1.1.2016",
         xlab = "",
         ylab = "migliaia")