Alternatively, you can use the MONOTONE function in PROC TRANSREG to fit a linear regression model with monotonicity constraints.