我正在运行2SLS模型。我需要考虑数据的群集,并使用了cluster.robust.se()。这段代码可以完美运行,但是我想知道cluster.robust.se()函数是在第一阶段还是第二阶段中或者仅在第二阶段中聚集标准错误。
我已经阅读了该函数的文档,但是很遗憾,这对我没有帮助。
library(ivpack)
iv1 <- ivreg(formula = Y ~ X1 | Z1, data=df)
iv1_clust <- cluster.robust.se(iv1, df$id)