> logistic <- glm(y~full_model, family=binomial(link="logit"))These are really, really bad numbers. The statistical method that 'should' account for interactions between variables actually makes things worse. (An AIC should be negative, and the residual deviance should be less than the degrees of freedom, not 200 times greater).
Null deviance: 20650253 on 860 degrees of freedom
Residual deviance: 118972 on 573 degrees of freedom
Unfortunately, the evolutionary process does not seem to employ the same types of linear thinking that our brains do. That makes it harder to tell a story people can understand.