Two-dimensional smoothing of scattered normal or non-normal (GLM) responses using tensor product P-splines.
Useful links