library(psych)
podatki1 <- force(sat.act)
summary(podatki1$ACT)
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 3.00 25.00 29.00 28.55 32.00 36.00
podatki1$Spol <- factor(podatki1$gender,
levels = c(1, 2),
labels = c("M", "F"))
library(psych)
describeBy(podatki1$ACT, podatki1$Spol)
##
## Descriptive statistics by group
## group: M
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 247 28.79 5.06 30 29.23 4.45 3 36 33 -1.06 1.89 0.32
## ------------------------------------------------------------
## group: F
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 453 28.42 4.69 29 28.63 4.45 15 36 21 -0.39 -0.42 0.22