https://www.rhosignal.com/posts/fitting-linear-models-with-Polars/