Linear regression

Might be related to some "numpy parity" ticket?