Using AIC regression (AIC stepwise gave the same result as well as other criterions:R-square adjusted, Mallow\'s Cp, SBC, SBIC)
library(olsrr) model <- lm(