r - 在 fixst::etable() 中对固定效果进行分组

Question

当报告来自 fixest 的结果时回归，我想通过在同一行中报告多个 FE 来节省空间。看着trade数据，这里是一个简单的例子:

data(trade)
est_1 <- feols(Euros ~ log(dist_km) | Origin + Destination + Product, trade)

etable产生我们所期望的:

> etable(est_1)
                                          est_1
Dependent Var.:                           Euros
                                               
log(dist_km)    -66,710,904.2*** (14,306,388.4)
Fixed-Effects:   ------------------------------
Origin                                      Yes
Destination                                 Yes
Product                                     Yes
_______________  ______________________________
S.E.: Clustered                      by: Origin
Observations                             38,325
R2                                      0.30324
Within R2                               0.05559
> 

但我希望 Origin 和 Destination 以及 Product 只是一行，名称类似于“Origin、Destination 和 Product FEs”。但是 group (和 keep 和 drop )似乎不适用于固定效果:

> etable(est_1, group = list("Origin and Destination and Product FEs" = "Origin|Dest|Product"))
                                                                 est_1
Dependent Var.:                                                  Euros
                                                                      
log(dist_km)                           -66,710,904.2*** (14,306,388.4)
Origin and Destination and Product FEs                              No
Fixed-Effects:                          ------------------------------
Origin                                                             Yes
Destination                                                        Yes
Product                                                            Yes
______________________________________  ______________________________
S.E.: Clustered                                             by: Origin
Observations                                                    38,325
R2                                                             0.30324
Within R2                                                      0.05559

即由 group 创建的新行都是没有，旧的 FE 行仍然存在。如果我试图欺骗 fixest通过在 setFixest_dict 中重复变量标签

setFixest_dict(c(Origin = "Origin and Destination and Product FEs",
                 Destination = "Origin and Destination and Product FEs",
                 Product = "Origin and Destination and Product FEs"))

etable明智地对我大喊大叫，因为我没有独特的值(value)观。
(我在这里使用的 trade 示例可能很糟糕，但请考虑对某些治疗和一系列人口固定效应(年龄、性别、地区等)进行回归工资。在这种情况下，是/否 FE称为 Demographic FEs 的行(在表格注释中适当解释)会很棒。)

Answer 1

不是很优雅，但您可以使用 etable() 的 data.frame返回。

library(fixest)
suppressMessages(library(dplyr))
data(trade)
est_1 <- feols(Euros ~ log(dist_km) | Origin + Destination + Product, trade)

# revemo all but first fixef_vars 
vars_to_remove <- est_1$fixef_vars[2:length(est_1$fixef_vars)]

# set a dict 
setFixest_dict(c(Origin = "Origin and Destination and Product FEs"))

# workd with data.frame
et <- etable(est_1)
et %>% filter(!rownames(et) %in% vars_to_remove)
#>                                                                  est_1
#> Dependent Var.:                                                  Euros
#>                                                                       
#> log(dist_km)                           -66,710,904.2*** (14,306,388.4)
#> Fixed-Effects:                          ------------------------------
#> Origin and Destination and Product FEs                             Yes
#> ______________________________________  ______________________________
#> S.E.: Clustered                                             by: Origin
#> Observations                                                    38,325
#> R2                                                             0.30324
#> Within R2                                                      0.05559

创建于 2021-03-23 由 reprex package (v1.0.0)