tune_xrnet function

k-fold cross-validation for hierarchical regularized regression