如何使用强大的群集标准错误运行逻辑和/或序数逻辑回归？
我认为Zelig会这样做 - 但robust和cluster功能已不再可用！
示例（非常简单）
ordinal <- zelig(Y ~ X1+X2, model = "ologit", robust=TRUE, **cluster**="id", data = df))
当然还有其他方法可以做到这一点。例如，在robcov包中使用rms函数。但是如果你运行500个不同的模型（比如我现在做的那样），这太费时了。那么有没有其他方法可以获得类似于Zelig以前的强大的集群标准错误？
非常感谢提前！