## Cross validation of hierarchical linear regression model
data(GaussianExample)
## 5-fold cross validation
cv_xrnet <- tune_xrnet(
x = x_linear,
y = y_linear,
external = ext_linear,
family = "gaussian",
control = xrnet_control(tolerance = 1e-6)
)
## Get coefficient estimates at optimal penalty combination
coef_opt <- coef(cv_xrnet)
Run the code above in your browser using DataLab