stargazer的多列输出用于knitr

Question

这是一个最小工作示例 (MWE)：

library(pscl)

data("bioChemists", package = "pscl")

fm_pois <- glm(art ~ ., data = bioChemists, family = poisson)
fm_qpois <- glm(art ~ ., data = bioChemists, family = quasipoisson)
fm_nb <- glm.nb(art ~ ., data = bioChemists)
fm_zinb <- zeroinfl(art ~ . | 1, data = bioChemists, dist = "negbin")

library(stargazer)
stargazer(
  fm_pois, fm_qpois,  fm_nb, fm_zinb
  , type = "text"
)

=============================================================================
                                      Dependent variable:                    
                  -----------------------------------------------------------
                                              art                            
                   Poisson   glm: quasipoisson     negative     zero-inflated
                                link = log         binomial      count data  
                     (1)            (2)              (3)             (4)     
-----------------------------------------------------------------------------
femWomen          -0.225***      -0.225***        -0.216***       -0.216***  
                   (0.055)        (0.074)          (0.073)         (0.073)   

marMarried         0.155**        0.155*            0.150*         0.150*    
                   (0.061)        (0.083)          (0.082)         (0.082)   

kid5              -0.185***      -0.185***        -0.176***       -0.176***  
                   (0.040)        (0.054)          (0.053)         (0.053)   

phd                 0.013          0.013            0.015           0.015    
                   (0.026)        (0.036)          (0.036)         (0.036)   

ment               0.026***      0.026***          0.029***       0.029***   
                   (0.002)        (0.003)          (0.003)         (0.003)   

Constant           0.305***       0.305**           0.256*         0.256*    
                   (0.103)        (0.139)          (0.137)         (0.139)   

-----------------------------------------------------------------------------
Observations         915            915              915             915     
Log Likelihood    -1,651.056                      -1,561.958     -1,560.959  
theta                                          2.264*** (0.271)              
Akaike Inf. Crit. 3,314.113                       3,135.917                  
=============================================================================
Note:                                             *p<0.1; **p<0.05; ***p<0.01

我正在寻找像这样的多列输出：

```html

我正在寻找像这样的多列输出：

```

=============================================================================
                                      Dependent variable:                    
                  -----------------------------------------------------------
                                              art                            
                           Poisson              Negative Binomial 

                    Poisson     QuasiPoisson         NB           ZINB

                     (1)            (2)              (3)             (4)     
-----------------------------------------------------------------------------
femWomen          -0.225***      -0.225***        -0.216***       -0.216***  
                   (0.055)        (0.074)          (0.073)         (0.073)   

marMarried         0.155**        0.155*            0.150*         0.150*    
                   (0.061)        (0.083)          (0.082)         (0.082)   

kid5              -0.185***      -0.185***        -0.176***       -0.176***  
                   (0.040)        (0.054)          (0.053)         (0.053)   

phd                 0.013          0.013            0.015           0.015    
                   (0.026)        (0.036)          (0.036)         (0.036)   

ment               0.026***      0.026***          0.029***       0.029***   
                   (0.002)        (0.003)          (0.003)         (0.003)   

Constant           0.305***       0.305**           0.256*         0.256*    
                   (0.103)        (0.139)          (0.137)         (0.139)   

-----------------------------------------------------------------------------
Observations         915            915              915             915     
Log Likelihood    -1,651.056                      -1,561.958     -1,560.959  
theta                                          2.264*** (0.271)              
Akaike Inf. Crit. 3,314.113                       3,135.917                  
=============================================================================
Note:                                             *p<0.1; **p<0.05; ***p<0.01

1. 第一行前两列应该是 Poisson ，下两列是Negative Binomial。
2. 第二行的列名称应该是Poisson，Quasi Poisson，Negative Binomial和Zero Inflated Negative Binomial。

我找到了这个链接，但不知道如何得到这个表格。

- MYaseen208

---

我认为使用stargazer或texreg都无法完成这个任务。你需要修改源代码以生成所需的输出。除非你经常这样做，否则手动编辑输出可能更容易些。 - Nick Kennedy

Answer 1

像 Nick Kennedy 一样，我认为 stargazer 不能直接生成您想要的输出。

因此，这里是一个解决方法：将 stargazer 表保存在一个对象中，并手动添加所需的行。 我在这里硬编码了它；用更多的努力应该可以自动居中相应列上方的文本。 请注意，我稍微修改了您的 stargazer 调用，以隐藏（错误的）模型名称。

library(pscl)
library(stargazer)

data("bioChemists", package = "pscl")

fm_pois <- glm(art ~ ., data = bioChemists, family = poisson)
fm_qpois <- glm(art ~ ., data = bioChemists, family = quasipoisson)
fm_nb <- glm.nb(art ~ ., data = bioChemists)
fm_zinb <- zeroinfl(art ~ . | 1, data = bioChemists, dist = "negbin")

byLine <- 
  do.call("c", 
          strsplit(
            capture.output(
              stargazer(fm_pois, fm_qpois,  fm_nb, fm_zinb, 
                        type = "text", model.names = FALSE)
              ),
            "\n"))

result <- append(
  byLine, 
  c(
  "                        Poisson               Negative Binomial",
  "",
  "                   Poisson  QuasiPoisson      NB          ZINB"
  ), 
  after = c(4, 5, 6))

cat(paste(result, collapse = "\n"))

# ==================================================================
#                                 Dependent variable:               
#                   ------------------------------------------------
#                                         art                       
#                         Poisson               Negative Binomial
# 
#                    Poisson  QuasiPoisson      NB          ZINB
#                      (1)        (2)          (3)           (4)    
# ------------------------------------------------------------------
# femWomen          -0.225***  -0.225***    -0.216***     -0.216*** 
#                    (0.055)    (0.074)      (0.073)       (0.073)  
#                                                                   
# marMarried         0.155**    0.155*        0.150*        0.150*  
#                    (0.061)    (0.083)      (0.082)       (0.082)  
#                                                                   
# kid5              -0.185***  -0.185***    -0.176***     -0.176*** 
#                    (0.040)    (0.054)      (0.053)       (0.053)  
#                                                                   
# phd                 0.013      0.013        0.015         0.015   
#                    (0.026)    (0.036)      (0.036)       (0.036)  
#                                                                   
# ment               0.026***  0.026***      0.029***      0.029*** 
#                    (0.002)    (0.003)      (0.003)       (0.003)  
#                                                                   
# Constant           0.305***   0.305**       0.256*        0.256*  
#                    (0.103)    (0.139)      (0.137)       (0.139)  
#                                                                   
# ------------------------------------------------------------------
# Observations         915        915          915           915    
# Log Likelihood    -1,651.056              -1,561.958    -1,560.959
# theta                                  2.264*** (0.271)           
# Akaike Inf. Crit. 3,314.113               3,135.917               
# ==================================================================
# Note:                                  *p<0.1; **p<0.05; ***p<0.01

- CL.