STATA: Genera descritivos rápidos con bases de datos gigantes!

data_stata_loop.R

library(tidyverse)
library(haven)
library(labelled)

#crear data frame
set.seed(22)

dta <- data_frame("edad" = floor(rnorm(6,25,10)), 
                  "peso" = floor(rnorm(6,60,10)), 
                  "sexo" = gl(2,1,6, labels = c("f","m")),
                  "altura" = gl(2,2,6, labels = c("bajo","alto")),
                  "outcome" = rep(c("primaria","secundaria","universidad"),each=2)) 
#regla
!grepl("^[A-Za-z_]{1}[A-Za-z0-9_]{0,31}$", names(dta))
#fuente:
#https://github.com/tidyverse/haven/blob/master/R/haven.R

#con error
dta %>% 
  mutate_if(is.character,as.factor) %>% 
  write_dta("data/test.dta")

stata_loop.do

*importar base de datos
use "data/test.dta"

*ver los tipos de variable
describe

*crear un "vector/lista" de variables según el tipo
*ejecutar un loop con funciones descriptivas
*la línea if+condition+continue permite evitar ciertas variables dentro de la lista creada

*loop para variables categóricas (guardadas con tipo long)

ds, has(type long)
foreach var in `r(varlist)' {
    if "`var'" == "outcome" continue
    codebook `var', detail
    tab `var' outcome, chi col
}

*loop para variables numéricas (guardadas con tipo double)

ds, has(type double)
foreach var in `r(varlist)' {
    codebook `var', detail
    tabstat `var', s(n min max mean sd q sk k) by(outcome)
    oneway `var' outcome, bonferroni
    *kwallis2 `var', by(community1)
}