Get z-scores of coefficients using glmnet

Question

I am using the glmnet package to get LASSO estimates, as follows:

model <- cv.glmnet(data, return, standardize = TRUE)

I am able to extract the coefficients using coef(model), however, I couldn't figure out a way of getting the standard error and Z-score for each variable.