An homotopy method for p regression provably beyond self-concordance and in input-sparsity time

Abstract

We consider the problem of linear regression where the 2n norm loss (i.e., the usual least squares loss) is replaced by the pn norm. We show how to solve such problems up to machine precision in O*(n|1/2 - 1/p|) (dense) matrix-vector products and O*(1) matrix inversions, or alternatively in O*(n|1/2 - 1/p|) calls to a (sparse) linear system solver. This improves the state of the art for any p∈ \1,2,+∞\. Furthermore we also propose a randomized algorithm solving such problems in input sparsity time, i.e., O*(Z + poly(d)) where Z is the size of the input and d is the number of variables. Such a result was only known for p=2. Finally we prove that these results lie outside the scope of the Nesterov-Nemirovski's theory of interior point methods by showing that any symmetric self-concordant barrier on the pn unit ball has self-concordance parameter (n).

0

Turn this paper into a lesson

ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…