访问Zelig relogit输出中的z值和其他统计信息
我想计算罕见事件的logit回归。我决定使用访问Zelig relogit输出中的z值和其他统计信息,r,stargazer,r-zelig,R,Stargazer,R Zelig,我想计算罕见事件的logit回归。我决定使用Zelig包(relogit函数)来实现这一点 通常,我使用stargazer提取并保存回归结果。但是,这两个包()似乎存在兼容性问题 现在,我想从Zeligrelogit输出中提取以下信息: 系数,z值,p值,观测次数,对数似然,AIC 我已经设法提取了p值和系数。然而,剩下的我都失败了。但我确信这些值必须能够以某种方式访问,因为它们在summary()输出中报告(但是,我没有将summary输出存储为R对象)。不能以与常规glm摘要相同的方式处理摘
Zelig
包(relogit
函数)来实现这一点
通常,我使用stargazer
提取并保存回归结果。但是,这两个包()似乎存在兼容性问题
现在,我想从Zelig
relogit
输出中提取以下信息:
系数,z值,p值,观测次数,对数似然,AIC
我已经设法提取了p值和系数。然而,剩下的我都失败了。但我确信这些值必须能够以某种方式访问,因为它们在summary()
输出中报告(但是,我没有将summary
输出存储为R
对象)。不能以与常规glm
摘要相同的方式处理摘要()
一个可复制的示例:
##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
使用_zelig_模型中的
进行偏差,AIC
m <- from_zelig_model(z.out1)
m$aic
...
谢谢你,真管用!只有两个后续问题:1)你认为z.out1.dat.loglik吗1)不,不幸的是,我不知道如何得到这里的可能性。偏差是使用饱和模型和您的模型的loglik计算的,但这里没有给出饱和模型的偏差。2)您是指用于回归的行数吗?如果是这样的话,length(complete.cases(m$model))
应该是一种安全的方法,谢谢!你的解决方案非常有效,所以我接受你的答案。很高兴它有帮助!
m <- from_zelig_model(z.out1)
m$aic
...
z.out1$get_coef()[[1]]/z.out1$get_se()[[1]]