Working with Generalized Additive Models (GAMs) in R: A Deep Dive into Smoothness Parameters and Choosing Between `method = "gam"` and `k` for Best Fit