Straightforward intermediate rank tensor product smoothing in mixed models

Simon N. Wood, Fabian Scheipl, Julian J. Faraway

Research output: Contribution to journalArticlepeer-review

74 Citations (SciVal)
290 Downloads (Pure)

Abstract

Tensor product smooths provide the natural way of representing smooth interaction terms in regression models because they are invariant to the units in which the covariates are measured, hence avoiding the need for arbitrary decisions about relative scaling of variables. They would also be the natural way to represent smooth interactions in mixed regression models, but for the fact that the tensor product constructions proposed to date are difficult or impossible to estimate using most standard mixed modelling software. This paper proposes a new approach to the construction of tensor product smooths, which allows the smooth to to be written as the sum of some fixed effects and some sets of i.i.d. Gaussian random effects: no previously published construction achieves this. Because of the simplicity of this random effects structure, our construction is useable with almost any flexible mixed modelling software, allowing smooth interaction terms to be readily incorporated into any Generalized Linear Mixed Model. To achieve the computationally convenient separation of smoothing penalties, the construction differs from previous tensor product approaches in the penalties used to control smoothness, but the penalties have the advantage over several alternative approaches of being explicitly interpretable in terms of function shape. Like all tensor product smoothing methods, our approach builds up smooth functions of several variables from {\em marginal} smooths of lower dimension, but unlike much of the previous literature we treat the general case in which the marginal smooths can be {\em any} quadratically penalized basis expansion, and there can be any number of them. We also point out that the imposition of identifiability constraints on smoothers requires more care in the mixed model setting than it would in a simple additive model setting, and show how to deal with the issue. An interesting side effect of our construction is that an ANOVA-decomposition of the smooth can be read off from the estimates, although this is not our primary focus. We were motivated to undertake this work by applied problems in the analysis of abundance survey data, and two examples of this are presented.
Original languageEnglish
Pages (from-to)341-360
Number of pages20
JournalStatistics and Computing
Volume23
Early online date11 Feb 2012
DOIs
Publication statusPublished - 1 May 2012

Fingerprint

Dive into the research topics of 'Straightforward intermediate rank tensor product smoothing in mixed models'. Together they form a unique fingerprint.

Cite this