Abstract
In this paper, a cubic integral smoothing spline with roughness penalty for restoring a function by integrals is described. A mathematical method for building such a spline is described in detail. The method is based on cubic integral spline with a penalty function, which minimizes the sum of squares of the difference between the observed integrals of the unknown function and the integrals of the spline being constructed, plus an additional penalty for the nonlinearity (roughness) of the spline. This method has a matrix form, and this paper shows in detail how to fill in each matrix. The parameter \( \alpha \) governs the desired smoothness of the restored function. Spline knots can be chosen independently of observations, and a weight can be defined for each observation for more control over the resulting spline shape. An implementation in the R language as function int_spline is given. The function int_spline is easy to use, with all arguments completely described and corresponding examples given. An example of the application of the method in rare event analysis and forecasting is given.
Supplemental Material
Available for Download
Software for Restoration of Function by Integrals with Cubic Integral Smoothing Spline in R
- [1] . 2001. A Practical Guide to Splines (Revised ed.). Springer Verlag, Berlin.Google Scholar
- [2] . 1985. Numerical solutions of the adsorption integral equation utilizing the spline functions. Thin Solid Films 123, 3 (1985), 245–272. Google ScholarCross Ref
- [3] . 2016. Application of integrodifferential splines to solving an interpolation problem. Computational Mathematics and Mathematical Physics. 54, 3 (2016), 19–34. Google ScholarCross Ref
- [4] . 2018. ibs: Integral of B-Spline Functions. https://CRAN.R-project.org/package=ibs.
R package version 1.4. Google Scholar - [5] . 1976. Numerical algorithm for recursively-defined convolution integrals involving distribution functions. Management Science 22, 10 (1976), 1138–1146.Google ScholarDigital Library
- [6] . 1978. Smoothing noisy data with spline functions - Estimating the correct degree of smoothing by the method of generalized cross-validation. Numer. Math. 31, 4 (1978), 377–403. Google ScholarDigital Library
- [7] . 2019. Spline quasi-interpolating projectors for the solution of nonlinear integral equations. J. Comput. Appl. Math. 354 (2019), 360–372. Google ScholarDigital Library
- [8] . 2016. Method of creation of a spline with the integral equal to integral of function of two variables on area of its definition. Science Almanac (In Russian) 33, 1–3(15) (2016), 31–35.Google Scholar
- [9] . 2011. Edge Functions for Spectral Element Methods.
Lecture Notes in Computational Science and Engineering , Vol. 76 LNCSE. 199–207. Google ScholarCross Ref - [10] . 1994. Nonparametric Regression and Generalized Linear Models. A Roughness Penalty Approach. Chapman & Hall, New York.Google ScholarCross Ref
- [11] . 1992. Analysis of discrete ill-posed problems by means of the L-curve. SIAM Rev. 34, 4 (1992), 561–580. Google ScholarDigital Library
- [12] . 2001. The L-curve and its use in the numerical treatment of inverse problems. (2001), 119–142.Google Scholar
- [13] . 2012. The theory and application of penalized methods or reproducing kernel Hilbert spaces made easy. Statistics Surveys 6, 1 (2012), 113–141.Google Scholar
- [14] . 2014. High order geometric methods with exact conservation properties. J. Comput. Phys. 257, PB (2014), 1444–1471. Google ScholarDigital Library
- [15] . 2012. Linear and nonlinear inverse problems with practical applications. Society for Industrial and Applied Mathematics, Philadelphia. Google ScholarCross Ref
- [16] . 1971. Some results on Tchebychefian spline functions. Journal Of Mathematical Analysis and Applications 33, 1 (1971), 82–95. Google ScholarCross Ref
- [17] . 2014. Integro-Differential Information Processing Method and Its Application in Numerical Analysis. IPI RAS (In Russian), Moscow.Google Scholar
- [18] . 2016. System structure of the economy: Qualitative time-space analysis. Fronteiras 5, 2 (2016), 61–81. Google ScholarCross Ref
- [19] . 2015. Capacity method determination consumption rate function. Ekonomika i Menedzhment Sistem Upravleniya (In Russian) 15, 1.1 (2015), 140–150.Google Scholar
- [20] . 2017. The causes analysis and error estimation of the anomalous pictures in the capacity method for the analysis of rare events. Ekonomika i Upravlenie: Problemy, Resheniya (In Russian) 8, 6 (2017), 8–12.Google Scholar
- [21] . 2019. Error of the capacity method of rare events analysis, remoteness from the end user. The News of KBSC of RAS. (In Russian) 3, 89 (2019), 48–77. Google ScholarCross Ref
- [22] . 2020. Accuracy research of the capacity method for analyzing rare events from uncertainty within the process of event formation (Part 1). The News of KBSC of RAS. (In Russian) 97, 5 (2020), 49–67. Google ScholarCross Ref
- [23] . 2020. The function restoration method by integrals for analysis and forecasting of rare events in the economy. Economics and Mathematical Methods. (In Russian) 56, 3 (2020), 113–124. Google ScholarCross Ref
- [24] . 1984. Methods for Solving Incorrectly Posed Problems. Springer-Verlag, New York. Google ScholarCross Ref
- [25] . 2003. Circle approximation using integral B-splines. Computer-Aided Design 35, 6 (2003), 601–607. Google ScholarCross Ref
- [26] . 1991. A fast efficient technique for the estimation of frequency. Biometrika 78, 3 (1991), 489–497.Google ScholarCross Ref
- [27] . 2001. The Estimation and Tracking of Frequency. Cambridge University Press, Cambridge.Google ScholarCross Ref
- [28] . 1991. Piecewise integral splines of low degree. Computers and Geosciences 17, 9 (1991), 1255–1263. Google ScholarDigital Library
- [29] . 1967. Smoothing by spline functions. Numer. Math10 (1967), 177–183.Google ScholarDigital Library
- [30] . 2001. Spline collocation for singular integro-differential equation. Izv. Vyssh. Uchebn. Zaved. Mat. (In Russian). 54, 11 (2001), 35–45.Google Scholar
- [31] . 2018. splines2: Regression Spline Functions and Classes. https://CRAN.R-project.org/package=splines2.
R package version 0.2.8. Google Scholar - [32] . 2004. Reconstruction of signal from samples of its integral in spline subspaces. Lecture Notes in Computer Science (including subseries Lecture Notes in Artificial Intelligence and Lecture Notes in Bioinformatics)3037 (2004), 574–577.Google Scholar
- [33] . 2020. A new B-spline collocation method for singular integro-differential equations of higher orders. J. Comput. Appl. Math. 380 (2020), 112949. Google ScholarCross Ref
Index Terms
- Algorithm 1023: Restoration of Function by Integrals with Cubic Integral Smoothing Spline in R
Recommendations
Integral form of the spline function in approximation problems
The article examines a new (integral) approach to the construction of a spline-approximation function. The proposed approach simplifies the process of spline construction. The integral form of the spline yields an analytical representation of the spline ...
Shape-preserving approximation of multiscale univariate data by cubic L1 spline fits
Spline fits are calculated by minimizing a data fitting functional over a manifold of splines. Smoothing splines are calculated by minimizing a linear combination of a data fitting functional and an interpolating spline functional. For multiscale data, ...
Restoration of the Product Consumption Rate with Integral Cubic Smoothing Spline, Study of the Best Smoothing Parameter Choice
AbstractThis article shows how the consumption rate function of a certain customer can be restored from a sequence of discrete purchases. For this, purchases are considered as integrals of the unknown consumption function with some uncertainties. To ...
Comments