Question

请原谅我，因为我是R的新手，如果这很傻/容易，但是我一直在找几个小时，但无济于事。我有一系列的GLM模型，我想在Stargazer表中报告每个模型的标准化/重新参数化系数以及直接系数。我创建了两个单独的模型，其中一个使用arm包具有标准化系数。

require(arm)
model1 <- glm(...)
model1.2 <- standardize(model1)

两个模型都可以工作，找到并提供我想要的输出，但是我似乎无法找到让Stargazer模仿此结构/外观的方法： http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.65.699&rep=rep1&type=pdf

Answer 1

这可以通过要求观星者为两个模型产生输出并确保模型中的系数具有相同的名称来完成，以便观星者知道标准化系数和非标准化系数应位于同一行。

下面的代码应该可以帮助您入门。

# generate fake data
x <- runif(100)
y <- rbinom(100, 1, x)

# fit the model
m1 <- glm(y~x, family = binomial())

# standardize it
m2 <- arm::standardize(m1)

# we make sure the coefficients have the same names in both models
names(m2$coefficients) <- names(m1$coefficients)

# we feed to stargazer
stargazer::stargazer(m1, m2, type = "text", 
                     column.labels = c("coef (s.e.)", "standarized coef (s.e.)"), 
                     single.row = TRUE)
#> 
#> ===========================================================
#>                              Dependent variable:           
#>                   -----------------------------------------
#>                                       y                    
#>                      coef (s.e.)    standarized coef (s.e.)
#>                          (1)                  (2)          
#> -----------------------------------------------------------
#> x                 4.916*** (0.987)     2.947*** (0.592)    
#> Constant          -2.123*** (0.506)      0.248 (0.247)     
#> -----------------------------------------------------------
#> Observations             100                  100          
#> Log Likelihood         -50.851              -50.851        
#> Akaike Inf. Crit.      105.703              105.703        
#> ===========================================================
#> Note:                           *p<0.1; **p<0.05; ***p<0.01

^{由reprex package（v0.2.1）于2019-02-13创建}

通常，您可以通过查找占星家帮助文件并查看此有用的网页https://www.jakeruss.com/cheatsheets/stargazer/

来弄清楚如何从输出中实现所需的结果

在Stargazer表中包含标准化系数

1 个答案: