R Stargazer 报告系数、置信区间和精确 p 值

Question

我正在运行一个国家固定效应模型，该模型还包括一个虚拟变量。 我正在使用 Stargazer 软件包，但似乎无法弄清楚如何报告置信区间和确切的 p 值。如果我运行我的模型：

stargazer(dummy_CPP, title = "xx",align = TRUE, 
          no.space = TRUE, ci = TRUE, 
          report = ('vc*p'), single.row = TRUE)

我的输出中只有 p 值。此外，由于我有 30 多个变量，整个表格无法放入 1 页，因为 p 值在系数下方报告，即使我使用 single.row = TRUE。

Answer 1

来自stargazer 文档；

report = 仅包含“v”、“c”、“s”、“t”、“p”、“*”元素的字符串，用于确定变量名称（“v” )、系数 ("c")、标准误差/置信区间 ("s")、检验统计量 ("t") 和 p 值 ("p") 应在回归表中报告。

因此报告 CI 和 p 值包括字母“s”和“p”。不过，我不知道如何在一行中获得所有信息（命令single.row 仅适用于标准错误/CI，但不适用于我认为的 p 值）。

stargazer(dummy_CPP, title = "xx",align = TRUE, 
          no.space = TRUE, ci = TRUE, 
          report = ('vcsp'), single.row = TRUE)

Answer 2

不使用 stargazer 而是使用我的 huxtable 包：

library(plm)
library(dplyr)
library(huxtable)
dfr <- tibble( id = rep(1:10, 10), x = rnorm(100), y = x + id/10 + rnorm(100), time = rep(1:10, each = 10))

dummy_CPP <- plm(y ~ x, dfr, index = c("id", "time"), effect = "twoways")

huxreg(dummy_CPP, error_pos = 'right', error_format = "({std.error}) [{conf.low} - {conf.high}]", ci_level = 0.95, statistics = "nobs")
──────────────────────────────────────────────────
                          (1)                     
         ─────────────────────────────────────────
  x          1.100 ***   (0.093) [0.917 - 1.282]  
         ─────────────────────────────────────────
  nobs     100                                    
──────────────────────────────────────────────────
  *** p < 0.001; ** p < 0.01; * p < 0.05.         

您可以在error_format 字符串中包含带有{p.value} 的p 值，并且您可以使用标准R 子集或通过更改font_size 来操作表格，以使其适合页面或拆分为多个页面.