Доступ к z-значению и другой статистике в выходных данных Zelig relogit
Я хочу вычислить логит-регрессию для редких событий. Я решил использовать Zelig
пакет (relogit
функция), чтобы сделать это.
Обычно я использую stargazer
извлечь и сохранить результаты регрессии. Однако, похоже, что есть проблемы совместимости с этими двумя пакетами ( Использование Stargazer с Zelig).
Теперь я хочу извлечь следующую информацию из Zelig
relogit
выход:
Коэффициенты, значения 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]]