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多个变量，整个表格无法放在一页上，即使我使用single.row = TRUE，p值也会显示在系数下方。

- AnnaSonder

---

如果您提供一些玩具数据和一个 [mcve]，那么帮助您会更容易。 - user3603486

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)

- JWilliman

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来操作表格，以使其适合页面或拆分为多个页面。

- user3603486

---

谢谢dash2！由于输出将包含在我的学士论文中，因此其设计非常重要。如何将其优雅地插入Lyx中？我希望有像Stargazer这样提供LaTex编码的工具。 - AnnaSonder

只需使用 to_latex() 或 print_latex()。或者使用 options(huxtable.print=print_latex) 自动打印它。 - user3603486