Trying to fit a submodel, What do I do wrong? lm_full <- lm(CO ~ AT + AP + AH + GTEP + TIT + TAT + TEY, data = Data)