title: Median-based linear regression in Python
body: <p>I would like to perform one-dimensional linear regression by minimizing the median absolute error.</p>
  
  <p>While initially assuming that it should be a fairly standard use case, a quick search surprisingly revealed that all regression and interpolation functions use the mean squared error.</p>
  
  <p>Therefore my question: is there a function that can perform median error based linear regression for one dimension?</p>