Блог

Я хочу ввести средние предельные эффекты в stargazer таблицу в R.

Проблема в том, что stargazer рассматривает первое значение предельных эффектов как константу в таблице. Средние предельные эффекты, которые я рассчитываю, не включают константу. Поэтому я хочу исключить его из таблицы. Использование omit аргумента в stargazer просто полностью исключает этот термин. Поскольку он рассматривает первое значение предельных эффектов как константу, он фактически снижает коэффициент.

Как я могу заставить звездочета опустить константу, заставляя его распознавать, что первое значение на самом деле является коэффициентом, а не константой?

Вот пример проблемы:

 #devtools::install_github("vincentarelbundock/marginaleffects")
library(marginaleffects)
library(stargazer)

Y <- sample(c(0,1), 1000, replace = T)
X <- sample(seq(1,50,1), 1000, replace = T)
Z <- sample(seq(1,50,1), 1000, replace = T)

DF <- data.frame(Y,X, Z)
Test <- glm(Y ~ X   Z, data = DF, family=binomial())
Test_DF <- data.frame(summary(marginaleffects(Test)))

Test_DF

#********************* Table **********************
#Fake Regression - to Make Table
Fake_RE <- glm(Y ~ X   Z, data = DF, family=binomial())

stargazer(list(Fake_RE), type = "text",
          coef = list(Test_DF$estimate),
          se = list(Test_DF$std.error))

 

Средние Предельные Эффекты:

 > Test_DF
      type term     estimate   std.error statistic   p.value     conf.low    conf.high
1 response    X  0.001540192 0.001109148  1.388627 0.1649462 -0.000633697 0.0037140816
2 response    Z -0.001215714 0.001087828 -1.117561 0.2637546 -0.003347817 0.000916389

 

Последующая Таблица:

 > stargazer(list(Fake_RE), type = "text",
            coef = list(Test_DF$estimate),
            se = list(Test_DF$std.error))

=============================================
                      Dependent variable:    
                  ---------------------------
                               Y             
---------------------------------------------
X                           -0.001           
                            (0.001)          
                                             
Z                                            
                                             
                                             
Constant                     0.002           
                            (0.001)          
                                             
---------------------------------------------
Observations                 1,000           
Log Likelihood             -690.097          
Akaike Inf. Crit.          1,386.194         
=============================================
Note:             *p<0.1; **p<0.05; ***p<0.01

 

stargazer Пакет официально не поддерживает объекты этого типа, и поскольку этот пакет не обновлялся с 2018 года, я не думаю, что разумно ожидать поддержки в ближайшее время.

Если вы готовы рассмотреть альтернативный пакет, вы, возможно, захотите попробовать modelsummary пакет (отказ от ответственности: я-автор):

 #devtools::install_github("vincentarelbundock/marginaleffects")
library(marginaleffects)
library(modelsummary)

Y <- sample(c(0,1), 1000, replace = T)
X <- sample(seq(1,50,1), 1000, replace = T)
Z <- sample(seq(1,50,1), 1000, replace = T)

DF <- data.frame(Y,X, Z)
Test <- glm(Y ~ X   Z, data = DF, family=binomial())
mfx <- marginaleffects(Test)

modelsummary(mfx)