test datasets
data.test <- data
data.test$z <- 1 - data.test$z
fit <- stan4bart(