Доступ к z-значению и другой статистике в выходных данных Zelig relogit

Я хочу вычислить логит-регрессию для редких событий. Я решил использовать Zelig пакет (relogit функция), чтобы сделать это.

Обычно я использую stargazer извлечь и сохранить результаты регрессии. Однако, похоже, что есть проблемы совместимости с этими двумя пакетами ( Использование Stargazer с Zelig).

Теперь я хочу извлечь следующую информацию из Zeligrelogit выход:

Коэффициенты, значения z, значения p, количество наблюдений, логарифмическая вероятность, AIC

Мне удалось извлечь р-значения и коэффициенты. Однако мне не удалось в остальном. Но я уверен, что эти значения должны быть как-то доступны, потому что они указаны в summary() вывод (однако мне не удалось сохранить summary вывод как R объект). Резюме не может быть обработано так же, как обычный glm резюме ( https://stats.stackexchange.com/questions/176821/relogit-model-from-zelig-package-in-r-how-to-get-the-estimated-coefficients)

Воспроизводимый пример:

##Initiate package, model and data

require(Zelig)

data(mid)

z.out1 <- zelig(conflict ~ major + contig + power + maxdem + mindem + years,
                data = mid, model = "relogit")

##Call summary on output (reports in console most of the needed information)

summary(z.out1)

##Storing the summary fails and only produces a useless object

summary(z.out1) -> z.out1.sum

##Some of the output I can access as follows

z.out1$get_coef() -> z.out1.coeff
z.out1$get_pvalue() -> z.out1.p
z.out1$get_se() -> z.out1.se

Однако я не нашел аналогичных команд для других элементов, таких как z values, AIC и т.д. Однако, как они показаны в summary() позвоните, они должны быть доступны как-то.

summary результат вызова:

Model: 

Call:
z5$zelig(formula = conflict ~ major + contig + power + maxdem + 
    mindem + years, data = mid)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-3.0742  -0.4444  -0.2772   0.3295   3.1556  

Coefficients:
             Estimate Std. Error z value Pr(>|z|)
(Intercept) -2.535496   0.179685 -14.111  < 2e-16
major        2.432525   0.157561  15.439  < 2e-16
contig       4.121869   0.157650  26.146  < 2e-16
power        1.053351   0.217243   4.849 1.24e-06
maxdem       0.048164   0.010065   4.785 1.71e-06
mindem      -0.064825   0.012802  -5.064 4.11e-07
years       -0.063197   0.005705 -11.078  < 2e-16

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 3979.5  on 3125  degrees of freedom
Residual deviance: 1868.5  on 3119  degrees of freedom
AIC: 1882.5

Number of Fisher Scoring iterations: 6

Next step: Use 'setx' method

1 ответ

Решение

Использование from_zelig_model за девиацию, АПК.

m <- from_zelig_model(z.out1)
m$aic
...

Z-значения являются коэффициент / SD.

z.out1$get_coef()[[1]]/z.out1$get_se()[[1]]
Другие вопросы по тегам