如何在fixest的etable中转换系数与标准误差以报告发生率比?
高效实现fixest回归结果与发生率比同表展示
- 核心方案:利用
etable()内置的转换参数直接处理,无需修改原回归对象,彻底避免z值等统计量失效的问题 - 实操代码示例:
- 先拟合基准模型
library(fixest) # 模拟测试数据 set.seed(123) dat <- data.frame( y = rnorm(100), x = rnorm(100), fe_grp = sample(1:10, 100, replace = TRUE) ) # 拟合固定效应模型 mod <- feols(y ~ x | fe_grp, data = dat) - 分栏展示原结果与发生率比
etable(mod, mod, # 指定每列的系数转换逻辑:第一列保留原值,第二列指数化 coef.transform = list(identity, exp), # 对应标准误的转换逻辑 se.transform = list(identity, exp), # 自定义列名 colnames = c("原回归结果", "发生率比"), # 添加说明注释 note = "发生率比为系数与标准误的指数形式(e^b、e^s)" )
- 先拟合基准模型
- 参数关键说明:
coef.transform/se.transform:接受函数列表,按列指定转换规则,identity表示不做变换,exp实现指数化- 此方法下,z值、p值等统计量仍基于原回归的推断逻辑,不会因为指数化操作失效
- 进阶:同一列内并列展示原系数与发生率比
如果想在同一系数行下同时显示原值和指数值,可借助split.coef参数实现:etable(mod, split.coef = list(x = c("原系数", "发生率比")), coef.transform = list(x = c(identity, exp)), se.transform = list(x = c(identity, exp)) )
内容的提问来源于stack exchange,提问作者gaur_banarasi
相关产品推荐
相关产品推荐

