我使用glm.nb函数创建了一个负二项式广义线性模型：
glm_mod <- glm.nb(obs ~ offset(log(time_interval)) + pred1*pred2, data = df)
我想基于pred1，pred2及其相互作用对速率（以单位时间间隔进行的观测）建模。这部分似乎工作正常。
接下来，我想使用pred1:pred2包绘制两个预测变量effects之间的交互作用。这是我遇到麻烦的原因，因为目前看来该程序包无法处理公式中定义的偏移量，并且glm.nb函数不接受和偏移变量作为参数，就像标准glm函数那样
effects和glm.nb之间是否存在这种不兼容的解决方法？
还有其他我可以使用的相关软件包吗？