Revisiting fitting monotone polynomials to data

Kevin Murray*, Samuel Müller, Berwin A. Turlach

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

14 Citations (Scopus)

Abstract

We revisit Hawkins' (Comput Stat 9(3):233-247, 1994) algorithm for fitting monotonic polynomials and discuss some practical issues that we encountered using this algorithm, for example when fitting high degree polynomials or situations with a sparse design matrix but multiple observations per x-value. As an alternative, we describe a new approach to fitting monotone polynomials to data, based on different characterisations of monotone polynomials and using a Levenberg-Marquardt type algorithm. We consider different parameterisations, examine effective starting values for the non-linear algorithms, and discuss some limitations. We illustrate our methodology with examples of simulated and real world data. All algorithms discussed in this paper are available in the R Development Core Team (A language and environment for statistical computing, R Foundation for Statistical Computing, Vienna, 2011) package MonoPoly.

Original languageEnglish
Pages (from-to)1989-2005
Number of pages17
JournalComputational Statistics
Volume28
Issue number5
DOIs
Publication statusPublished - Oct 2013
Externally publishedYes

Keywords

  • Monotone polynomial
  • Monotone regression

Fingerprint

Dive into the research topics of 'Revisiting fitting monotone polynomials to data'. Together they form a unique fingerprint.

Cite this