bm_mars function

Fit Multivariate Adaptive Regression Splines models