Abstract:Representation of generalized additive models (GAM's) using penalized regression splines allows GAM's to be employed in a straightforward manner using penalized regression methods. Not only is inference facilitated by this approach, but it is also possible to integrate model selection in the form of smoothing parameter selection into model fitting in a computationally efficient manner using well founded criteria such as generalized cross-validation. The current fitting and smoothing parameter selection methods for such models are usually effective, but do not provide the level of numerical stability to which users of linear regression packages, for example, are accustomed. In particular the existing methods cannot deal adequately with numerical rank deficiency of the GAM fitting problem, and it is not straightforward to produce methods that can do so, given that the degree of rank deficiency can be smoothing parameter dependent. In addition, models with the potential flexibility of GAM's can also present practical fitting difficulties as a result of indeterminacy in the model likelihood: Data with many zeros fitted by a model with a log link are a good example. In this article it is proposed that GAM's with a ridge penalty provide a practical solution in such circumstances, and a multiple smoothing parameter selection method suitable for use in the presence of such a penalty is developed. The method is based on the pivoted QR decomposition and the singular value decomposition, so that with or without a ridge penalty it has good error propagation properties and is capable of detecting and coping elegantly with numerical rank deficiency. The method also allows mixtures of user specified and estimated smoothing parameters and the setting of lower bounds on smoothing parameters. In terms of computational efficiency, the method compares well with existing methods. A simulation study compares the method to existing methods, including treating GAM's as mixed models.

An introduction to modeling longitudinal data with generalized additive models: Applications to single-case designs.

Modeling Multivariate Count Time Series Data with a Vector Poisson Log-Normal Additive Model: Applications to Testing Treatment Effects in Single-Case Designs

Linear mixed effects models for non‐Gaussian continuous repeated measurement data

Generalized additive models with flexible response functions

Statistical Inference for Generalized Additive Models: Simultaneous Confidence Corridors and Variable Selection

Estimating Structural Mean Models with Multiple Instrumental Variables Using the Generalised Method of Moments

Methods for Modeling Autocorrelation and Handling Missing Data in Mediation Analysis in Single Case Experimental Designs (SCEDs)

Missing not at random models for latent growth curve analyses.

Multilevel modeling in single-case studies with zero-inflated and overdispersed count data

On Bayesian Generalized Additive Models

Selecting a linear mixed model for longitudinal data: Repeated measures analysis of variance, covariance pattern model, and growth curve approaches.

Generalized additive models: Building evidence of air pollution, climate change and human health

Estimating curvilinear time-varying treatment effects: Combining g-estimation of structural nested mean models with time-varying effect models for longitudinal causal inference.

Bayesian Nonparametric Regression Analysis of Data with Random Effects Covariates from Longitudinal Measurements

Generalized Estimating Equations in Longitudinal Studies: A Non-Parametric Alternative for Two-Way Repeated Measures Mixed ANOVA

Longitudinal Modeling of Age-Dependent Latent Traits with Generalized Additive Latent and Mixed Models

Generalised additive mixed models for dynamic analysis in linguistics: a practical introduction

Structural Nested Mean Models Under Parallel Trends Assumptions

Statistical modelling of COVID-19 data: Putting generalized additive models to work

Stable and Efficient Multiple Smoothing Parameter Estimation for Generalized Additive Models

A generalized additive model approach to time-to-event analysis