Confidence intervals from local minimums of objective function
Abstract
The weighted median plays a central role in the least absolute deviations (LAD). We propose a nonlinear regression using (LAD). Our objective function $f(a, l, s)$ is non-convex with respect to the parameters a, l, s, and is such that for each fixed l, s the minimizer of $a\to f (a, l, s)$ is the weighted median $med(x(l, s), w(l, s))$ of a sequence $x(l, s)$ endowed with the weights $w(l, s)$ (all depend on $l$, $s$). We analyse and compare theoretically the minimizers of the function $(a, l, s)\to f (a, l, s)$ and the surface $(l, s) \to f (med(x(l, s), w(l, s)), l, s)$. As a numerical application we propose to fit the daily infections of COVID 19 in China using Gaussian model. We derive confident interval for the daily infections from each local minimum.References
Barmparis, G. D. and Tsironis, G. P. (2020). Estimating the infection horizon of COVID-19 in eight countries with a data-driven approach, Chaos, Solitons and Fractals, 135, 109842.
Bloomfield, P. and Steiger, W. L. (1980). Least Absolute Deviations Curve-Fitting. SIAM J. Sci. Stat. Comput., 2, 290-301.
Dielman, T., E. (2005). Least absolute value regression: recent contributions. J. Stat. Comput. Simul., 75, 263-286.
Eisenhart, C. (1961). Boscovitch and the combination of observations. In Roger Joseph Boscovitch, ed. L.L. Whyte. New York: Fordham University Press.
Gao, F. and Lixing H. (2010). Implementing the Nelder-Mead simplex algorithm with adaptive parameters. Comput. Optim. Appl., 51, 259-277.
Karst, O. J. (1958). Linear curve fitting using least deviations. J. Amer. Stat. Assoc., 53, 118-132.
Koenker, R. and Bassett, G. (1985). On Boscovich's estimator. Ann. Statist., 13, 1625-1628.
Lagarias, J. C., Reeds, J. A, Wright, M.W and Wright, P. E. (1998). Convergence properties of the Nelder--Mead simplex method in low dimensions. SIAM J. Optim., 9, 112-147.
Li, Y.; Arce, G. R. (2004). A maximum likelihood approach to least absolute deviation regression. EURASIP J. Appl. Signal Process. 12, 1762-1769.
Narula, S.C. and Wellington, J.F. (1977). Algorithm AS 108. Multiple linear regression with minimum sum of absolute errors. J. R. Stat. Soc. Ser. C. Appl. Stat., 26, 106-111.
Nelder J. A and Mead R. (1965). A simplex method for function minimization. Comput. J., 7, 308--313.
Novoselac, V. (2020). The Component Weighted Median Absolute Deviations Problem. Eur. J. Pure Appl. Math., 13, 964-976.
Rhodes, E. C. (1930). Reducing observations by the method of minimum deviations. Phil Mag. Seventh Series, 9, 974-992.
Robert, C. (1995). Simulation of truncated normal variables. Stat. Comput., 5, 121-125.
Sadovski, A. N. (1974). Algorithm As 74. L1-norm fit of straight line. J. R. Stat. Soc. Ser. C. Appl. Stat., 23, 244-248.
- Authors retain copyright and grant the journal right of first publication with the work simultaneously licensed under a Creative Commons Attribution License that allows others to share the work with an acknowledgement of the work's authorship and initial publication in this journal.
- Authors are able to enter into separate, additional contractual arrangements for the non-exclusive distribution of the journal's published version of the work (e.g., post it to an institutional repository or publish it in a book), with an acknowledgement of its initial publication in this journal.
- Authors are permitted and encouraged to post their work online (e.g., in institutional repositories or on their website) prior to and during the submission process, as it can lead to productive exchanges, as well as earlier and greater citation of published work (See The Effect of Open Access).