如何在R中使用stargazer在同一行中輸出多個變量

Question

我想從同一行中的幾個回歸中輸出交互項，並將其稱為“交互”。 到目前為止，我所擁有的是交互術語顯示在稱為“交互”的兩個不同行中（請參閱下面的代碼）。

這個問題已經在這里問過了，但我的分數還不夠高，無法對其進行投票或評論： https : //stackoverflow.com/questions/28859569/several-coefficients-in-one-line 。

library("stargazer")
stargazer(attitude)
stargazer(attitude, summary=FALSE)
# 2 OLS models with Interactions
linear.1 <- lm(rating ~ complaints + privileges + complaints*privileges
           , data=attitude)
linear.2 <- lm(rating ~ complaints + learning + complaints*learning, data=attitude)
stargazer(linear.1, linear.2, title="Regression Results", type="text", 
      covariate.labels=c("Complaints", "Privileges", "Interaction", "Learning", "Interaction"))

感謝您的幫助。

Answer 1

我認為這不是本機支持的，因為這不是一個好主意。 您要求混淆表格中數字的含義，這對您的讀者沒有幫助。

現在聲明的警告，您可以通過修改lm對象的內容來做到這一點：

# copy objects just for demonstration
m1 <- linear.1
m2 <- linear.2

# see names of coefficients
names(m1$coefficients)
# [1] "(Intercept)"           "complaints"            "privileges"            "complaints:privileges"
names(m2$coefficients)
# [1] "(Intercept)"         "complaints"          "learning"            "complaints:learning"

# replace names
names(m1$coefficients)[names(m1$coefficients) == "complaints:privileges"] <- "interaction"
names(m2$coefficients)[names(m2$coefficients) == "complaints:learning"] <- "interaction"

結果：

> stargazer(m1, m2, title="Regression Results", type="text")

Regression Results
==========================================================
                                  Dependent variable:     
                              ----------------------------
                                         rating           
                                   (1)            (2)     
----------------------------------------------------------
complaints                       1.114**         0.307    
                                 (0.401)        (0.503)   

privileges                        0.434                   
                                 (0.570)                  

learning                                        -0.171    
                                                (0.570)   

interaction                       -0.007         0.006    
                                 (0.008)        (0.009)   

Constant                          -7.737        31.203    
                                 (27.409)      (31.734)   

----------------------------------------------------------
Observations                        30            30      
R2                                0.692          0.713    
Adjusted R2                       0.657          0.680    
Residual Std. Error (df = 26)     7.134          6.884    
F Statistic (df = 3; 26)        19.478***      21.559***  
==========================================================
Note:                          *p<0.1; **p<0.05; ***p<0.01

Answer 2

以下回應：

reg ~ felm(....)
rownames(reg$coefficients)[rownames(reg$coefficients)=='oldname']<-'newname'
rownames(reg$beta)[rownames(reg$beta)=='oldname']<-'newname'

似乎適用於大多數情況。

雖然我有時會遇到問題。 當 IV 與 felm 一起使用時，這是必需的。 雖然區分適合和不適合 IV 的變量很好，但與其他模型相比，這些表格會很麻煩！ 這種語法是有幫助的。

Answer 3

如果有人想知道，我需要它用於felm包的不同目的。 為此需要以下代碼：

reg ~ felm(....)
rownames(reg$coefficients)[rownames(reg$coefficients)=='oldname']<-'newname'
rownames(reg$beta)[rownames(reg$beta)=='oldname']<-'newname'