Modeling Heterogeneity With Bayesian Additive Regression Trees