More efficient approximation of smoothing splines via space-filling basis selection

We consider the problem of approximating smoothing spline estimators in a nonparametric regression model. When applied to a sample of size , the smoothing spline estimator can be expressed as a linear combination of basis functions, requiring computational time when the number of predictors . Such a sizable computational cost hinders the broad applicability of smoothing splines. In practice, the full sample smoothing spline estimator can be approximated by an estimator based on randomly-selected basis functions, resulting in a computational cost of . It is known that these two estimators converge at the identical rate when is of the order , where depends on the true function , and depends on the type of spline. Such is called the essential number of basis functions. In this article, we develop a more efficient basis selection method. By selecting the ones corresponding to roughly equal-spaced observations, the proposed method chooses a set of basis functions with a large diversity. The asymptotic analysis shows our proposed smoothing spline estimator can decrease to roughly , when . Applications on synthetic and real-world datasets show the proposed method leads to a smaller prediction error compared with other basis selection methods.
View on arXiv