Argonne National Laboratory

Adaptively Refined Dynamic Program for Linear Spline Regression

TitleAdaptively Refined Dynamic Program for Linear Spline Regression
Publication TypeJournal Article
Year of Publication2012
AuthorsGoldberg, N, Kim, Y, Leyffer, S, Veselka, T
Date Published09/2012
Other NumbersANL/MCS-P3040-0912

The linear spline regression problem is to determine a piecewise linear function for estimating a set of given points while minimizing a given measure of misfit or error. This is a classical problem in computational statistics and operations research; dynamic programming was proposed as a solution technique more than forty years ago by Bellman and Roth (1969). The algorithm requires a discretization of the solution space to define a grid of candidate breakpoints. This paper proposes an adaptive refinement scheme for the grid of candidate breakpoints in order to allow the dynamic programming method to scale for larger instances of the problem. We evaluate the quality of solutions found on small instances compared with optimal solutions determined by a novel integer programming formulation of the problem. We also consider a generalization of the linear spline regression problem to fit multiple curves that share break point horizontal coordinates, and we extend our method to solve the generalized problem. Computational experiments verify that our nonuniform grid construction schemes are useful for computing high-quality solutions for both the single-curve and two-curve linear spline regression problem.