我目前在MRS中使用rxLogit模型作为MRO中标准GLM模型的替代方案（约300,000行，但有3个因子，分别为200,400和5000）。这可以将我的计算时间从26小时减少到7小时，使用128 GB RAM的2x6核心Xeon。
估算后我需要/想要计算群集的健壮标准误差。在标准R中，我可以通过多种方式完成此操作，例如，通过multiwayvcoc包中的cluster.vcoc（）。
我尝试将带有as.glm（）的rxLogit对象转换为glm对象并随后执行我的分析，但后来我收到以下错误：
Error in qr.lm(object) : lm object does not have a proper 'qr' component.
Rank zero or should not have used lm(.., qr=FALSE).
调用时也会出现此错误，例如摘要（as.glm（rxLogitObject））。
有没有人遇到类似/相同的问题并且可以提供解决方案？
非常感谢，
尼古拉斯