[…] The GAMPL procedure fits generalized additive models. I have previously shown the power of the GAMPL procedure to fit binary response data. …