bm_cubist function

Fit Cubist models (M5)