“…Some previous work has estimated monotone exposure-response functions for exposure to air pollution or weather and a health outcome assessed on the same day (Powell et al, 2012;Wilson et al, 2014), and there is a rich statistical literature on shape constrained regression for a general regression function. Approaches for shape constrained regression in a general regression setting include piecewise linear functions (Hildreth, 1954;Brunk, 1955), kernel smoothers Mammen (1991), a large number of spine based approaches (Ramsay, 1988;Neelon and Dunson, 2004;Meyer, 2008;Wang and Li, 2008;Meyer et al, 2011;Meyer, 2012;Powell et al, 2012) and Bernstein polynomial methods (Chang et al, 2005(Chang et al, , 2007Curtis and Ghosh, 2011;Wilson et al, 2014;Ding and Zhang, 2016;Wilson et al, 2020). Chipman et al (2021) proposed a monotone regression model based on the popular nonparametric Bayesian additive regression trees (BART) framework (Chipman et al, 2010) that constrains a general regression function to be monotone with respect to some or all predictors.…”