多元逻辑回归模型的变量选择，包括R中的二次项和交互项

Question

当我在做一个项目时，我陷入了错误提示

glm.fit：发生数字0或1的拟合概率

当我使用完整模型尝试glm()时。在将二次项和交互项纳入模型之前，它可以正常工作。我该如何解决？

我的数据temp具有11个变量和1个响应变量（12列），具有201行（观察）。

以下是我的代码：

test2 <- glm(Y ~ (Age_Q16 + COPD + HighSch + Widodiv + Married + DISABLED + CEREBROVASCULAR_DISEASE + VALVE_DISEASE + NITRATES + Q9 + Q14Down)^2
             + Age_Q16sq + Q9sq + Q14Downsq,temp,family=binomial("logit"),maxit=100)

summary(test2)