【问题标题】:Report coefficients and exponentiated coefficients side by side with stargazer与 stargazer 并排报告系数和指数系数
【发布时间】:2016-03-04 11:08:29
【问题描述】:

如何在模型输出中添加一列,其中包括系数的取幂值旁边的系数的未取幂值?这是 add.lines 选项的用途吗?谢谢你的帮助。

#load library
library(stargazer)
#sample data
var1<-sample(c('men', 'women'), replace=TRUE, 100, prob=c(0.3, 0.7))
var2<-sample(c('left', 'right'), replace=TRUE, 100, prob=c(0.5, 0.5))
df<-data.frame(gender=var1, vote=var2)
#binomial logistic regression
model1<-glm(vote~gender, data=df, family='binomial')
#report model without
stargazer(model1)
#report model with exponentiated coefficients
stargazer(model1, apply.coef=exp)

【问题讨论】:

    标签: r latex stargazer


    【解决方案1】:

    可能有更好的方法,但是像这样的快速破解呢:

    model1Exp = model1 
    model1Exp$coefficients = exp(model1Exp$coefficients)
    
    stargazer(model1, model1Exp)
    

    您需要调整默认标签以标记第二个模型与第一个模型相同但具有指数系数的事实。

    【讨论】:

    • 我可以挖掘它。我将暂缓接受它作为答案,让某人有机会在观星者中展示一种方法。很聪明。
    • 所以,这行得通。然而,一个问题是在第二个模型中,它使用第一个模型的 p 值(即未取幂的系数与标准误差的比率)来决定是否在第二个模型中打印显着性星(系数已经指数化)。这会导致差异。我通过将参数添加到 stargazer 命令来解决此问题。 code p=list(coef(summary(model1))[,5], coef(summary(model1))[,5] code
    • 现在,唯一的事情是,我想真正弄清楚如何抑制第二个模型的标准误差,但要为第一个模型报告它们。这可能吗?
    猜你喜欢
    • 1970-01-01
    • 2016-12-26
    • 1970-01-01
    • 2017-10-28
    • 2019-07-06
    • 2013-04-12
    • 1970-01-01
    • 2015-10-11
    • 1970-01-01
    相关资源
    最近更新 更多