与 stargazer 并排报告系数和取幂系数
Report coefficients and exponentiated coefficients side by side with stargazer
如何向模型输出添加一列，其中包括系数的指数值旁边的系数的指数值？这是 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)
可能有更好的方法，但像这样的快速破解怎么样：
model1Exp = model1
model1Exp$coefficients = exp(model1Exp$coefficients)
stargazer(model1, model1Exp)
您需要调整默认标签以标记第二个模型与第一个模型相同，但具有指数系数这一事实。
如何向模型输出添加一列，其中包括系数的指数值旁边的系数的指数值？这是 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)
可能有更好的方法，但像这样的快速破解怎么样：
model1Exp = model1
model1Exp$coefficients = exp(model1Exp$coefficients)
stargazer(model1, model1Exp)
您需要调整默认标签以标记第二个模型与第一个模型相同，但具有指数系数这一事实。