L2WPMA
Fortran software is developed that calculates a best piecewise monotonic approximation to n univariate data contaminated by random errors. The underlying method minimizes the weighted sum of the squares of the errors by requiring k − 1 sign changes in the first divided differences of the approximation, where k is a given positive integer. Hence, the piecewise linear interpolant to the fit consists of k monotonic sections, alternately increasing and decreasing. This calculation can have about O(nk) local minima, because the positions of the turning points of the fit are integer variables of the problem. The method, however, by employing a dynamic programming technique divides the data into at most k disjoint sets of adjacent data and solves a k = 1 problem (monotonic fit or isotonic regression) for each set. So it calculates efficiently a global solution in only O(nσ + kσ2) computer operations when k ≥ 3, where σ is the number of local minima of the data, always bounded by n/2. This complexity reduces to only O(n) when k = 1 or k = 2 (unimodal case). At the end of the calculation a spline representation of the solution and the corresponding Lagrange multipliers are provided. The software package has been tested on a variety of data sets showing a performance that does provide in practice shorter computation times than the complexity indicates in theory. An application of the method on identifying turning points and monotonic trends of data from 1947--1996 on the U.K. pound over the U.S. dollar exchange rate is presented. Generally, the method may have useful applications as, for example, in estimating the turning points of a function from some noisy measurements of its values, or in image and signal processing, or in providing a preliminary or complementary smoothing phase to further analyses of the data.
(Source: http://dl.acm.org/)
This software is also peer reviewed by journal TOMS.
This software is also peer reviewed by journal TOMS.
Keywords for this software
References in zbMATH (referenced in 6 articles , 1 standard article )
Showing results 1 to 6 of 6.
Sorted by year (- Szusz, Emily K.; Willms, Allan R.: A linear time algorithm for near minimax continuous piecewise linear representations of discrete data (2010)
- Demetriou, I.C.; Vassiliou, E.E.: An algorithm for distributed lag estimation subject to piecewise monotonic coefficients (2009)
- Demetriou, Ioannis C.: Algorithm 863: L2WPMA, a fortran 77 package for weighted least-squares piecewise monotonic data approximation. (2007)
- García, Víctor M.; Vidal, V.; Verdú, G.; Garayoa, J.; Miró, R.: Parallel resolution of the two-group time dependent neutron diffusion equation with public domain ODE codes (2005)
- Vassiliou, E.; Demetriou, I.C.: An adaptive algorithm for least squares piecewise monotonic data Fitting (2005)
- Demetriou, I.C.: Signs of divided differences yield least squares data fitting with constrained monotonicity or convexity (2002)