example_data <- simulate_spareg_data(n = 200, p = 2000, ntest = 100)
spar_res <- spar(example_data$x, example_data$y, xval = example_data$xtest,
yval = example_data$ytest, nummods=c(5, 10, 15, 20, 25, 30),
rp = rp_gaussian(control = list(sd = 1/sqrt(ncol(example_data$x)))))
Run the code above in your browser using DataLab