inst/tinytest/stata/clean.R

In marginaleffects: Predictions, Comparisons, Slopes, Marginal Means, and Hypothesis Tests

library(haven)
library(ivreg)
library(betareg)
library(MASS)
library(data.table)
library(readxl)

results <- list()

get_stata <- function(name, n) {
    sprintf(here::here("inst/tinytest/stata/results/%s.txt"), name) |>
        read.table(sep = "\t", skip = 4) |>
        setNames(c("term", "dydxstata", "std.errorstata")) |>
        head(n) |>
        transform(dydxstata = as.numeric(dydxstata),
                  std.errorstata = as.numeric(std.errorstata))
}


#  stats::glm
tmp <- read.table(here::here("inst/tinytest/stata/results/stats_glm_01.txt"), sep = "\t", skip = 4) |>
    setNames(c("term", "dydxstata", "std.errorstata")) |>
    head(2)
results[["stats_glm_01"]] <- tmp

#  stats::lm
tmp <- read.table(here::here("inst/tinytest/stata/results/stats_lm_01.txt"), sep = "\t", skip = 4) |>
    setNames(c("term", "dydxstata", "std.errorstata")) |>
    head(2)
results[["stats_lm_01"]] <- tmp

#  survival::coxph
tmp <- read.table(here::here("inst/tinytest/stata/results/survival_coxph_01.txt"), sep = "\t", skip = 4) |>
    setNames(c("term", "dydxstata", "std.errorstata")) |>
    head(1)
results[["survival_coxph_01"]] <- tmp

#  truncreg::truncreg
tmp <- read.table(here::here("inst/tinytest/stata/results/truncreg_truncreg_01.txt"), sep = "\t", skip = 4) |>
    setNames(c("term", "dydxstata", "std.errorstata")) |>
    head(2)
results[["truncreg_truncreg_01"]] <- tmp

#  quantreg::rq
tmp <- read.table(here::here("inst/tinytest/stata/results/quantreg_rq_01.txt"), sep = "\t", skip = 4) |>
    setNames(c("term", "dydxstata", "std.errorstata")) |>
    head(4)
results[["quantreg_rq_01"]] <- tmp

#  pscl::zeroinfl
tmp <- read.table(here::here("inst/tinytest/stata/results/pscl_zeroinfl_01.txt"), sep = "\t", skip = 4) |>
    setNames(c("term", "dydxstata", "std.errorstata")) |>
    head(3)
results[["pscl_zeroinfl_01"]] <- tmp

#  MASS::polr 
tmp <- read.table(here::here("inst/tinytest/stata/results/MASS_polr_01.txt"), sep = "\t", skip = 4) |>
    setNames(c("group", "dydxstata_x1", "std.errorstata_x1", "dydxstata_x2", "std.errorstata_x2")) |>
    head(3) |>
    transform(group = 2:4) |>
    data.table()
cols <- names(tmp)
tmp[, (cols) := lapply(.SD, as.numeric)]
tmp <- melt(tmp, id.vars = "group")
tmp[, c("statistic", "term") := tstrsplit(variable, "_")]
tmp[, variable := NULL]
tmp <- dcast(tmp, group + term ~ statistic, value.var = "value")
results[["MASS_polr_01"]] <- tmp

#  nnet::multinom 
tmp <- read.table(here::here("inst/tinytest/stata/results/nnet_multinom_01.txt"), sep = "\t", skip = 4) |>
    setNames(c("group", "dydxstata_x1", "std.errorstata_x1", "dydxstata_x2", "std.errorstata_x2")) |>
    head(3) |>
    transform(group = 2:4) |>
    data.table()
cols <- names(tmp)
tmp[, (cols) := lapply(.SD, as.numeric)]
tmp <- melt(tmp, id.vars = "group")
tmp[, c("statistic", "term") := tstrsplit(variable, "_")]
tmp[, variable := NULL]
tmp <- dcast(tmp, group + term ~ statistic, value.var = "value")
results[["nnet_multinom_01"]] <- tmp

#  ivreg::ivreg
tmp <- read.table(here::here("inst/tinytest/stata/results/ivreg_ivreg_01.txt"), sep = "\t", skip = 4) |>
    setNames(c("term", "dydxstata", "std.errorstata")) |>
    head(2)
results[["ivreg_ivreg_01"]] <- tmp

#  betareg::betareg 
tmp <- read.table(here::here("inst/tinytest/stata/results/betareg_betareg_01.txt"), sep = "\t", skip = 4) |>
    setNames(c("term", "dydxstata", "std.errorstata")) |>
    head(1)
results[["betareg_betareg_01"]] <- tmp

#  MASS::glm.nb
tmp <- read.table(here::here("inst/tinytest/stata/results/MASS_glm_nb.txt"), sep = "\t", skip = 4) |>
    setNames(c("term", "dydxstata", "std.errorstata")) |>
    head(3)
tmp$term <- gsub("\\.cyl", "", tmp$term)
tmp$term[2:3] <- paste0("factor(cyl)", tmp$term)[2:3]
results[["mass_glm_nb"]] <- tmp

#  AER::tobit
tmp <- read.table(here::here("inst/tinytest/stata/results/AER_tobit.txt"), sep = "\t", skip = 4) |>
    setNames(c("term", "dydxstata", "std.errorstata")) |>
    head(5)
results[["aer_tobit"]] <- tmp

#  AER::tobit (right-censured)
tmp <- read.table(here::here("inst/tinytest/stata/results/AER_tobit_right4.txt"), sep = "\t", skip = 4) |>
    setNames(c("term", "dydxstata", "std.errorstata")) |>
    head(5)
results[["aer_tobit_right"]] <- tmp

#  estimatr::lm_robust
tmp <- read.table(here::here("inst/tinytest/stata/results/estimatr_lm_robust.txt"), sep = "\t", skip = 4) |>
    setNames(c("term", "dydxstata", "std.errorstata")) |>
    head(3)
results[["estimatr_lm_robust"]] <- tmp

#  estimatr::iv_robust
tmp <- read.table(here::here("inst/tinytest/stata/results/estimatr_iv_robust.txt"), sep = "\t", skip = 4) |>
    setNames(c("term", "dydxstata", "std.errorstata")) |>
    head(2)
tmp$std.errorstata <- as.numeric(gsub("\\(|\\)", "", tmp$std.errorstata))
results[["estimatr_iv_robust"]] <- tmp


# fixest::feols 
tmp <- read.table(here::here("inst/tinytest/stata/results/fixest_feols_01.txt"), sep = "\t", skip = 4) |>
    head(2) |>
    setNames(c("term", "dydxstata", "std.errorstata"))
tmp$dydxstata <- as.numeric(tmp$dydxstata)
results[["fixest_feols"]] <- tmp


# fixest::fepois 
tmp <- read.table(here::here("inst/tinytest/stata/results/fixest_fepois_01.txt"), sep = "\t", skip = 4) |>
    head(2) |>
    setNames(c("term", "dydxstata", "std.errorstata"))
tmp$dydxstata <- as.numeric(tmp$dydxstata)
results[["fixest_fepois"]] <- tmp

# plm
tmp <- read.table(here::here("inst/tinytest/stata/results/plm_pooling_01.txt"), sep = "\t", skip = 4) |>
    head(2) |>
    setNames(c("term", "dydxstata", "std.errorstata"))
tmp$dydxstata <- as.numeric(tmp$dydxstata)
results[["plm_pooling"]] <- tmp

tmp <- read.table(here::here("inst/tinytest/stata/results/plm_sa_01.txt"), sep = "\t", skip = 4) |>
    head(2) |>
    setNames(c("term", "dydxstata", "std.errorstata"))
tmp$dydxstata <- as.numeric(tmp$dydxstata)
results[["plm_sa"]] <- tmp



# lme4:
tmp <- read.table(here::here("inst/tinytest/stata/results/lme4_01.txt"), sep = "\t", skip = 4) |>
    head(2) |>
    setNames(c("term", "dydxstata", "std.errorstata"))
tmp$dydxstata <- as.numeric(tmp$dydxstata)
results[["lme4_lmer"]] <- tmp

# lme4:
tmp <- read.table(here::here("inst/tinytest/stata/results/lme4_02.txt"), sep = "\t", skip = 4) |>
    head(2) |>
    setNames(c("term", "dydxstata", "std.errorstata"))
tmp$dydxstata <- as.numeric(tmp$dydxstata)
results[["lme4_glmer"]] <- tmp

# elasticity
results[["stats_lm_elasticity_eyex"]] <- get_stata("stats_lm_elasticity_eyex", 2)
results[["stats_lm_elasticity_dyex"]] <- get_stata("stats_lm_elasticity_dyex", 2)
results[["stats_lm_elasticity_eydx"]] <- get_stata("stats_lm_elasticity_eydx", 2)

results[["stats_glm_elasticity_eyex"]] <- get_stata("stats_glm_elasticity_eyex", 2)
results[["stats_glm_elasticity_dyex"]] <- get_stata("stats_glm_elasticity_dyex", 2)
results[["stats_glm_elasticity_eydx"]] <- get_stata("stats_glm_elasticity_eydx", 2)

# save
saveRDS(results, file = here::here("inst/tinytest/stata/stata.rds"))