从cpglmm解释输出和预测

Question

我正在使用cpglmm包（link）中的cplm函数，帮助同事在R中拟合化合物-泊松广义线性混合模型。该模型涉及三向交互，我想计算一些可解释的数量。到目前为止，我已经尝试计算出几率，但是我不确定这是正确的方法吗？

# Fit model with three-way interaction in fixed effects #

m <- cpglmm(ncs ~ diversity_index*diversity_speciality*n_authors + selfcit + 
n_refs + (1|region), data = diversity)

# Calculate Odds-ratio #
se <- sqrt(diag(vcov(m)))

tab <- cbind(Est = m$fixef, 
             S.E. = se, 
             LL = m$fixef - 1.96 * se, 
             UL = m$fixef + 1.96 * se)

print(exp(tab), digits=3)

我还想计算一些预测值，例如预测概率等，但是我无法让predict()为cpglmm工作。我可以使用任何功能吗？