Confidence intervals from local minimums of objective function

  • Azzouz Dermoune Université of Lille 1
  • Daoud Ounaissi ESA, Angers
  • Yousri Slaoui université de Poitiers
Keywords: Gaussian model, Least absolute deviation, daily infection, Simplex algorithm, Nelder-Mead, Optim 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.

Published
2023-07-08
How to Cite
Dermoune, A., Ounaissi, D., & Slaoui, Y. (2023). Confidence intervals from local minimums of objective function. Statistics, Optimization & Information Computing, 11(4), 798-810. https://doi.org/10.19139/soic-2310-5070-1738
Section
Research Articles