r - 获得 plotmo plots 表示的概率方程

Question

我想获得由plotmo（R）表示的概率函数的方程。这是当改变一个或两个预测变量同时保持其他预测变量的平均值不变时的模型方程。我想要一种简单的方法来获得数学方程，因为必须对具有不同变量的许多模型进行制作。

如果我的模型是这样的：

glm(formula = pres_aus ~ pH_sp + Annual_prec + I(pH_sp^2) + I(Annual_prec^2), family = binomial(link = "logit"), data = puntos_calibrado)

我怎样才能做到？

Answer 1

没有提供数据示例，因此没有进行测试，但是您不能跳过符号表达式的构造并按照以下方式进行操作：

model.matrix(data.frame(one=1, dat) ) %*% coef(mdl.fit) 
# where mdl.fit is returned from glm()

在某种意义上，这是公式的 R 矩阵表示：sum(beta_i*X_1)。如果您想为特定列指定平均值，那么只需将该数据框拉开并仅使用其中的一部分进行计算。所以对于保持在平均值的第一列：

model.matrix(data.frame(one=1, mn1 =mean(dat[[1]]), dat[-1]) ) %*% 
                                                       coef(mdl.fit)