cubist.default function

Fit a Cubist model