交互效果图与speedglm输出

Question

我正在处理一个足够大的数据集，以便使用比stats glm更先进的算法来计算二项式回归模型。 glm需要一周的时间来计算。

Large fixed effects binomial regression in R建议rxGlm和speedglm。我尝试了两种方法，并设法实现了速度的显着提高（低至2小时以下）。但是，我面临以下问题：

在这两种情况下，我都无法使用effects函数plot(effect())创建类似于以下内容的交互图，因为effects无法处理这些类型缺少对象和/或所需数据。

我已经在rxGlm（How to plot interaction effects from extremely large data sets (esp. from rxGlm output)）上发布了一个关于赏金的问题，但无济于事。所以我想使用rxGlm无法做到这一点。

现在我的问题是：我可以使用speedglm 绘制此类互动效果字词，如果是，如何？