问题描述
Cluster.robust.se似乎改变了iv回归的数据类型,并通过观星仪减少了观察次数,有什么提示可以使它们重新投入使用?
示例:
library(ivpack)
library(stargazer)
df<- tibble(
indiv_id = sample(10,100,rep = TRUE),group_id = sample(10,cluster= sample(2,rep=TRUE)
)
iv<-cluster.robust.se(ivreg(df$indiv_id~df$group_id),df$cluster)
stargazer(iv,type="text")
解决方法
这是怎么回事,cluster.robust.se()
将模型的类从ivreg
更改为coeftest
。检查一下:
> class(ivreg(df$indiv_id~df$group_id))
[1] "ivreg"
> class(iv)
[1] "coeftest"
这是两个不同的对象,因此观星者会得出不同的结果,而从coeftest
输出中排除观察计数。
一种可能的解决方案是将两个对象都包含在同一表中,如下所示:
iv<-cluster.robust.se(ivreg(df$indiv_id~df$group_id),df$cluster)
myivreg <- ivreg(df$indiv_id~df$group_id)
stargazer(myivreg,iv,type="text")