On a class of L-splines of order 4: fast algorithms for interpolation and smoothing

摘要

In this paper a special class of one-dimensional L-splines of order 4 is studied, which naturally appear in the computation of interpolation and smoothing with multivariate polysplines. Fast algorithms are provided for interpolation and smoothing with this class of L-splines, as well as a generalization of the Reinsch algorithm to this setting. The explicit description of all mathematical expressions permits a simple and direct numerical implementation. Applications are provided to financial data of the index S&P500, for the fast calculation of statistically interesting quantities, as cross validation (scores), generalized cross validation (scores) for finding the best smoothing parameter alpha, and the residual sum of squares.