3  Descriptivos

3.1 Descriptivos básicos

De forma más sencilla, para empezar se puede hacer …

https://rpkgs.datanovia.com/rstatix/reference/get_summary_stats.html

Para todos los casos, cargar los siguientes paquetes …

pacman::p_load(rstatix,kableExtra,expss)
data <- read_spss("data/ejemplogbw8.sav")

Del mismo modo, se pueden aádir tras el listado de variables

3.1.1 Sin grupo y sin peso

desc.result <- get_summary_stats(data, P16)
kable(desc.result)
variable n min max median q1 q3 iqr mad mean sd se ci
P16 264 304.21 780.88 330.06 315.03 390.17 75.14 26.746 377.217 104.12 6.408 12.618

3.1.2 Con grupo y sin peso

desc.result <- get_summary_stats(data %>% group_by(SEXO), P16)
kable(desc.result)
SEXO variable n min max median q1 q3 iqr mad mean sd se ci
1 P16 156 304.21 780.88 336.07 318.03 397.682 79.653 35.657 387.673 116.128 9.298 18.366
2 P16 108 306.01 630.61 330.06 312.02 351.095 39.075 26.746 362.114 81.986 7.889 15.639

3.1.3 Sin grupo y con peso

dataw <- weight_by(data, weight=P15)
desc.result <- get_summary_stats(dataw, P16)
kable(desc.result)
variable n min max median q1 q3 iqr mad mean sd se ci
P16 884 304.21 780.88 342.08 318.03 420.22 102.19 40.104 389.526 110.016 3.7 7.262

3.1.4 Con grupo y con peso

dataw <- weight_by(data, weight=P15)
desc.result <- get_summary_stats(dataw %>% group_by(SEXO), P16)
kable(desc.result)
SEXO variable n min max median q1 q3 iqr mad mean sd se ci
1 P16 548 304.21 780.88 342.08 324.04 420.22 96.18 40.104 396.592 119.689 5.113 10.043
2 P16 336 306.01 630.61 330.06 318.03 420.22 102.19 26.746 378.001 91.087 4.969 9.775