请原谅我,因为我对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
这可以通过要求观星者为两个模型产生输出并确保模型中的系数具有相同的名称以使得观察者知道标准化系数和非标准化系数应该在同一行上来完成。
下面的代码可以帮助您入门。
# 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创建于2019-02-13(v0.2.1)
您通常可以通过深入了解stargazer帮助文件并查看这个有用的网页来了解如何从输出中获得您想要的内容https://www.jakeruss.com/cheatsheets/stargazer/