mvjump

multivariate version

mvjump(ic = 0, rate = 0.4, nc = 2, ...)