podatki <- data.frame("ID" = c(1, 2, 3, 4, 5),
"visina" = c(170, 175, 177, 180, 182),
"starost" = c(20, 22, 22, 23, 24),
"spol" = c("Z", "Z", "M", "M", "M"))
mean(podatki$visina)
## [1] 176.8
summary(podatki[ , c(-1, -4) ])
## visina starost
## Min. :170.0 Min. :20.0
## 1st Qu.:175.0 1st Qu.:22.0
## Median :177.0 Median :22.0
## Mean :176.8 Mean :22.2
## 3rd Qu.:180.0 3rd Qu.:23.0
## Max. :182.0 Max. :24.0
podatki$teza <- c(65, 60, 70, 72, 80)
podatki$BMI <- round(podatki$teza/(podatki$visina/100)^2, 2)
library(pastecs)
round(stat.desc(podatki[-1, -4]), 3)
## ID visina starost teza BMI
## nbr.val 4.000 4.000 4.000 4.000 4.000
## nbr.null 0.000 0.000 0.000 0.000 0.000
## nbr.na 0.000 0.000 0.000 0.000 0.000
## min 2.000 175.000 22.000 60.000 19.590
## max 5.000 182.000 24.000 80.000 24.150
## range 3.000 7.000 2.000 20.000 4.560
## sum 14.000 714.000 91.000 282.000 88.300
## median 3.500 178.500 22.500 71.000 22.280
## mean 3.500 178.500 22.750 70.500 22.075
## SE.mean 0.645 1.555 0.479 4.113 0.939
## CI.mean.0.95 2.054 4.947 1.523 13.089 2.987
## var 1.667 9.667 0.917 67.667 3.524
## std.dev 1.291 3.109 0.957 8.226 1.877
## coef.var 0.369 0.017 0.042 0.117 0.085
library(psych)
describe(podatki$visina)
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 5 176.8 4.66 177 176.8 4.45 170 182 12 -0.29 -1.73 2.08