r - 使用 GLM 模型中的多路集群鲁棒 vcov 预测 se.fit 以绘制适当的置信区间

Question

我想绘制我的概率模型的预测，并根据我使用 multiwaycov 包中的 cluster.vcov() 函数获得的双向聚类标准误差包括 95% 置信区间。

但是，那里的预测函数不允许我将正确的方差协方差矩阵传递给它们，并且似乎需要一种解决方法（请参阅Joris Mey和Michael在https://stackoverflow 上的 lm-models 的解决方法。 com/a/3793955）。

有没有人有一个想法去解决这个问题（glm-models）？或者，还有其他方法可以解决这个问题吗？