Non-convex isotonic regression via the Myersonian approach

Isotonic regression refers to a class of regression models with order constraints. It is widely used in maximum likelihood estimation of ordered parameter, testing of distributions with ordered means, multistage production systems, and machine learning. A vast majority of the literature considers the isotonic regression problems with convex or piece-wise convex objective functions (or those that can be converted to such functions). We connect a class of isotonic regression problems with the so-called ironing problem in mechanism design, by establishing a discrete version of the Myerson's ironing method. We use such a connection to solve an isotonic regression problem with non-convex objective functions. We also prove the optimality of pool adjacent violator (PAV) algorithm in such a case.