library(car)
Prestige[1, "income"] <- NA
numSummary(Prestige[,c("income", "education")],
statistics=c("mean", "sd", "quantiles", "cv", "skewness", "kurtosis"))
numSummary(Prestige[,c("income", "education")], groups=Prestige$type)
remove(Prestige)
Run the code above in your browser using DataLab