[Ether]

Given m-by-n matrix [A], (m>n), and m-by-1 column vector [b] for the overdetermined system of linear equations

[A][x] = [b],

What is the recommended way to use R to find the n-by-1 column vector [x] which minimizes the L1 norm of [b]-[A][x] ?

aka minx |[b]-[A][x]|₁, aka Least Absolute Deviation (LAD) ?