我正在尝试使用4个不同的预测变量来运行逻辑模型：深度（连续），DPM（2个级别），邻近度（3个级别）和区域（2个级别）。我知道这很多交互会让事情变得复杂，但是想知道是否有可能生成一个包含所有成对交互（和深度）的汇总表。非常感谢您的帮助。代码是：
x-coord.at(row) - x-coord.at(col)
y-coord.at(row) - y-coord.at(col)
我非常感谢任何帮助。
一些数据来修复想法：
mod <- glm(cbind(Towards, Total-Towards) ~ DPM*Proximity*Depth*Area,data=Data, family=binomial('logit'))