GLM是处理（数字IV +分类IV =二进制DV）数据集的最佳方法吗？如何绘制？

Question

我有一个包含1个数字IV + 5个分类IV的数据集，DV也是二进制的，因此我使用glmer探索了谓词效果：

$model <- glmer(Y~(Xn+Xc1+Xc2+Xc3+Xc3)^2+(1|id),data,family=binomial)

但这是最好的方法吗？以及如何绘制主要效果和相互作用？