Warning: file_get_contents(/data/phpspider/zhask/data//catemap/4/r/76.json): failed to open stream: No such file or directory in /data/phpspider/zhask/libs/function.php on line 167

Warning: Invalid argument supplied for foreach() in /data/phpspider/zhask/libs/tag.function.php on line 1116

Notice: Undefined index: in /data/phpspider/zhask/libs/function.php on line 180

Warning: array_chunk() expects parameter 1 to be array, null given in /data/phpspider/zhask/libs/function.php on line 181
访问Zelig relogit输出中的z值和其他统计信息_R_Stargazer_R Zelig - Fatal编程技术网

访问Zelig relogit输出中的z值和其他统计信息

访问Zelig relogit输出中的z值和其他统计信息,r,stargazer,r-zelig,R,Stargazer,R Zelig,我想计算罕见事件的logit回归。我决定使用Zelig包(relogit函数)来实现这一点 通常,我使用stargazer提取并保存回归结果。但是,这两个包()似乎存在兼容性问题 现在,我想从Zeligrelogit输出中提取以下信息: 系数,z值,p值,观测次数,对数似然,AIC 我已经设法提取了p值和系数。然而,剩下的我都失败了。但我确信这些值必须能够以某种方式访问,因为它们在summary()输出中报告(但是,我没有将summary输出存储为R对象)。不能以与常规glm摘要相同的方式处理摘

我想计算罕见事件的logit回归。我决定使用
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]]