Linear regression in Python: Use of numpy, scipy, and statsmodels

As we can see, the statsmodels library allows us to generate highly detailed output on a level similar to R, with additional statistics such as skew, kurtosis, R-Squared and AIC. While these readings can be generated through scipy or sklearn, doing so is a more intensive process and in many cases these statistics must be calculated individually.