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