Hi Ygor, excellent question! This is something I glossed over when researching this topic.

The most popular loss function for GAMs is Restricted Maximum Likelihood Estimation (REML). It seems in R's mgcv package, this is optimized using "Penalized Iteratively Re-weighted Least Squares" . I'm not sure what PyGAM does but I'd imagine its similar.

You can read the mgcv documentation here: https://www.rdocumentation.org/packages/mgcv/versions/1.8-35/topics/gam

--

--

Data Scientist | Get in touch: https://www.linkedin.com/in/adamshafi

Get the Medium app

A button that says 'Download on the App Store', and if clicked it will lead you to the iOS App store
A button that says 'Get it on, Google Play', and if clicked it will lead you to the Google Play store