Включите стандартизованные коэффициенты в таблицу Stargazer

Простите меня, поскольку я новичок в R, и если это глупо / легко, но я искал часами, но безрезультатно. У меня есть серия моделей GLM, и я хотел бы сообщить о стандартизованных / перепараметризованных коэффициентах для каждого вместе с прямыми коэффициентами в таблице Stargazer. Я создал две отдельные модели, одну со стандартизованными коэффициентами, используя пакет arm.

require(arm)
model1 <- glm(...)
model1.2 <- standardize(model1)

Обе модели работают, находят и выдают нужные мне результаты, но я не могу придумать способ заставить Stargazer имитировать эту структуру / внешний вид: https://i.stack.imgur.com/UwZNn.png


r coefficients stargazer
person Lauren Kahn    schedule 13.02.2019    source источник

Ответы (1)


arrow_upward
0
arrow_downward

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

Приведенный ниже код должен помочь вам начать работу.

# generate fake data
x <- runif(100)
y <- rbinom(100, 1, x)

# fit the model
m1 <- glm(y~x, family = binomial())

# standardize it
m2 <- arm::standardize(m1)

# we make sure the coefficients have the same names in both models
names(m2$coefficients) <- names(m1$coefficients)

# we feed to stargazer
stargazer::stargazer(m1, m2, type = "text", 
                     column.labels = c("coef (s.e.)", "standarized coef (s.e.)"), 
                     single.row = TRUE)
#> 
#> ===========================================================
#>                              Dependent variable:           
#>                   -----------------------------------------
#>                                       y                    
#>                      coef (s.e.)    standarized coef (s.e.)
#>                          (1)                  (2)          
#> -----------------------------------------------------------
#> x                 4.916*** (0.987)     2.947*** (0.592)    
#> Constant          -2.123*** (0.506)      0.248 (0.247)     
#> -----------------------------------------------------------
#> Observations             100                  100          
#> Log Likelihood         -50.851              -50.851        
#> Akaike Inf. Crit.      105.703              105.703        
#> ===========================================================
#> Note:                           *p<0.1; **p<0.05; ***p<0.01

Создано 13 февраля 2019 г. пакетом REPEX (v0.2.1)

Обычно вы можете понять, как добиться желаемого из результатов, погрузившись в файл справки Stargazer и просмотрев эту полезную веб-страницу https://www.jakeruss.com/cheatsheets/stargazer/

person gfgm    schedule 13.02.2019
comment
Большое спасибо!!! Я копался во всем вышеперечисленном, и мне не повезло с тем, что я пытался - я ценю чрезвычайно плодотворный ответ !! - person Lauren Kahn; 13.02.2019