I’ve been attempting to use the cv.glmnet function from the glmnet package in R when making a cv LASSO regression. I want to standardize my explanatory variables as they are all on differing scales. I first used the argument standardize = TRUE and then cross checked the answer by scaling my data beforehand. However, it appears that I’m getting different lambda and MSE results from both methods. Please see below.