Why not just use basic gbm?
Note that it is a bad idea to encode your multinomial outcome as an integer unless the function specifically requires it. You did ask for a "regression tree" at first but wanted a multinomial outcome, which implies classification rather than regression. Factors are the preferred method of storing a categorical outcome in R.
# this works:
data.gbm$Y <- factor(paste0("cls", data.gbm$Y))
fit.step <- gbm(
Y ~ .,
data = data.gbm,
distribution = "multinomial",
interaction.depth = 5,
shrinkage = 0.005,
bag.fraction = 1,
n.minobsinnode = 2
)
Note the argument names are consistent with gbm too.