在括号中显示p值而不是标准误差

Question

当使用stargazer包时，我想要更改在系数下方括号中出现的值。默认情况下，该包将输出标准误差。如何在括号中包含实际的p值？

- user3799604

---

有些选项需要您定义正在使用的数据列和表中模型的数量。没有可用的示例将不会很有帮助。 - IRTFM

Answer 1

如Stargazer Omit test statistics所述，自5.0版本以来，stargazer已包含report参数，允许用户选择要报告的统计数据。因此，要显示p值而不是SEs，您需要执行以下操作：

require(stargazer)
linear.1 <- lm(rating ~ complaints + privileges + learning 
                        + raises + critical, data=attitude)
## put in the p-values rather than the se's
stargazer(linear.1, type="text", report=('vc*p'))

这将输出：

> stargazer(linear.1,  type="text", report=('vc*p'))

========================================
                 Dependent variable:    
             ---------------------------
                       rating           
----------------------------------------
complaints            0.692***          
                     p = 0.0002         

privileges             -0.104           
                      p = 0.450         

learning                0.249           
                      p = 0.132         

raises                 -0.033           
                      p = 0.870         

critical                0.015           
                      p = 0.918         

Constant               11.010           
                      p = 0.357         

----------------------------------------
Observations             30             
R2                      0.715           
Adjusted R2             0.656           
F Statistic           12.060***         
========================================
Note:        *p<0.1; **p<0.05; ***p<0.01

这种方法比使用se参数更安全,不会影响显著性星号。

另请参阅：

• StarGazer中的不精确p值

- landroni

---

这个问题尚未得到回答：我能用括号显示这些p值吗？“p =”在期刊上比较不常见。 - Jakob

这很可能是@Marek的一个功能请求。您可以手动将SE替换为p值，如Stargazer中不精确的p值所示，但您必须非常小心处理t.auto和p.auto参数... - landroni

Answer 2

没有简单的方法（除非其中一个style选项解决了这个问题）。但是您可以用p值替换标准误差，并硬编码p值，以便正确数量的星号出现如下。

library(stargazer)
#> 
#> Please cite as:
#>  Hlavac, Marek (2018). stargazer: Well-Formatted Regression and Summary Statistics Tables.
#>  R package version 5.2.2. https://CRAN.R-project.org/package=stargazer

models <- list()
models[[1]] <- lm(mpg ~ cyl + disp, data = mtcars)
models[[2]] <- lm(mpg ~ cyl + disp + wt, data = mtcars)

get_ts <- function(fm) {
  summary(fm)$coefficients[,3]
}

get_pvals <- function(fm) {
  summary(fm)$coefficients[,4]
}

ts <- lapply(models, get_ts)
pvals <- lapply(models, get_pvals)

stargazer(models, type = "text", report=('vc*s'), se = pvals, p = pvals)
#> 
#> =================================================================
#>                                  Dependent variable:             
#>                     ---------------------------------------------
#>                                          mpg                     
#>                              (1)                    (2)          
#> -----------------------------------------------------------------
#> cyl                        -1.587**              -1.785***       
#>                            (0.034)                (0.007)        
#>                                                                  
#> disp                       -0.021*                 0.007         
#>                            (0.054)                (0.533)        
#>                                                                  
#> wt                                               -3.636***       
#>                                                   (0.002)        
#>                                                                  
#> Constant                  34.661***              41.108***       
#>                            (0.000)                (0.000)        
#>                                                                  
#> -----------------------------------------------------------------
#> Observations                  32                     32          
#> R2                          0.760                  0.833         
#> Adjusted R2                 0.743                  0.815         
#> Residual Std. Error    3.055 (df = 29)        2.595 (df = 28)    
#> F Statistic         45.808*** (df = 2; 29) 46.424*** (df = 3; 28)
#> =================================================================
#> Note:                                 *p<0.1; **p<0.05; ***p<0.01

^{本文档由reprex package (v2.0.0)于2021-05-11创建}

- Ian Gow