R语言 模型中的固定效应计数总结

93ze6v8z  于 2023-01-28  发布在  其他
关注(0)|答案(1)|浏览(227)

我有三个固定效应回归的modelsummary,如下所示:

remotes::install_github("lrberge/fixest")
remotes::install_github("vincentarelbundock/modelsummary")

library(fixest)
library(modelsummary)

mod1 <- feols(mpg ~ hp | cyl, data = mtcars)
mod2 <- feols(mpg ~ wt | cyl, data = mtcars)
mod3 <- feols(mpg ~ drat | cyl, data = mtcars)

modelsummary(list(mod1, mod2, mod3), output = "markdown")

| | 模型1|模型2|模型3|
| - ------|- ------|- ------|- ------|
| 马力|-0.024|||
| | (0.015)|||
| 重量||-3.206||
| | | (1.188)||
| 德拉特|||一点七九三|
| | | | (1.564)|
| Num.Obs.|三十二|三十二|三十二|
| R2|0.754|0.837|零点七四五|
| R2调整|零点七二七|0.820|零点七一八|
| R2范围内|0.080|0.392|0.048|
| R2调整范围内|0.047|0.371|0.014|
| AIC公司|一百六十七点九|一百五十四点六|一百六十九点零万|
| BIC|一百七十三点八|一百六十点五|一百七十四点九|
| 均方根误差|二点九四|二点三十九分|二点九九|
| Std.Errors|通过:cyl|通过:cyl|通过:cyl|
| FE:气缸|十|十|十|

    • 除了让表格只显示某些固定效应是否存在之外,是否有可能显示估计的固定效应的数量?**

原始模型包含以下信息:

> mod1
OLS estimation, Dep. Var.: mpg
Observations: 32 
Fixed-effects: cyl: 3
Standard-errors: Clustered (cyl) 
    Estimate Std. Error  t value Pr(>|t|) 
hp -0.024039   0.015344 -1.56664  0.25771 
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
RMSE: 2.94304     Adj. R2: 0.727485
                Within R2: 0.07998
wgx48brx

wgx48brx1#

是的,你需要定义一个glance_custom.fixest()方法。请参阅文档的这一节以获得详细的说明和许多示例:
https://vincentarelbundock.github.io/modelsummary/articles/modelsummary.html#customizing-existing-models-part-i
下面是一个fixest的例子:

library(fixest)
library(tibble)
library(modelsummary)

models <- list(
    feols(mpg ~ hp | cyl, data = mtcars),
    feols(mpg ~ hp | am, data = mtcars),
    feols(mpg ~ hp | cyl + am, data = mtcars)
)

glance_custom.fixest <- function(x, ...) {
    tibble::tibble(`# FE` = paste(x$fixef_sizes, collapse = " + "))
}

modelsummary(models, gof_map = c("nobs", "# FE"))

| | (1) |(二)|(三)|
| - ------|- ------|- ------|- ------|
| 马力| -0.024 | -0.059 | -0.044 |
| | (0.015)|(千)|(0.016)|
| Num.Obs. |三十二|三十二|三十二|
| FE数量|三个|第二章|三加二|

相关问题