将系数与标准偏差相乘

Question

在R中，stargazer包提供了将函数应用于系数、标准误差等的可能性：

dat <- read.dta("http://www.ats.ucla.edu/stat/stata/dae/nb_data.dta")
dat <- within(dat, {
    prog <- factor(prog, levels = 1:3, labels = c("General", "Academic", "Vocational"))
    id <- factor(id)
})
m1 <- glm.nb(daysabs ~ math + prog, data = dat)
transform_coef <- function(x) (exp(x) - 1)
stargazer(m1, apply.coef=transform_coef)

我如何应用一个函数，其中乘以的因子取决于变量，例如该变量的标准差？

Answer 1

这可能不是您所希望的，但您可以变换系数，并为 stargazer 自定义 list 系数。例如，如果您想报告系数乘以每个变量的标准差，您的示例的以下扩展可能有效：

library(foreign)
library(stargazer)
library(MASS)

dat <- read.dta("http://www.ats.ucla.edu/stat/stata/dae/nb_data.dta")
dat <- within(dat, {
  prog <- factor(prog, levels = 1:3, labels = c("General", "Academic", "Vocational"))
  id <- factor(id)
})
m1 <- glm.nb(daysabs ~ math + prog, data = dat)

# Store coefficients (and other coefficient stats)
s1 <- summary(m1)$coefficients

# Calculate standard deviations (using zero for the constant)
math.sd  <- sd(dat$math)
acad.sd  <- sd(as.numeric(dat$prog == "Academic"))
voc.sd   <- sd(as.numeric(dat$prog == "Vocational"))
int.sd   <- 0

# Append standard deviations to stored coefficients
StdDev   <- c(int.sd, math.sd, acad.sd, voc.sd)
s1       <- cbind(s1, StdDev)

# Store custom list
new.coef <- s1[ , "Estimate"] * s1[ , "StdDev"]

# Output
stargazer(m1, coef = list(new.coef))

您可能需要考虑一些关于在 stargazer 中输出系数的原始问题之外的问题。在乘以标准偏差时是否应该报告截距？你的标准误和推论会和这个转换一样吗？

将系数与标准偏差相乘

Multiply coefficients with standard deviation

r

stargazer