I have just run the following linear regression using sklearn:
model4 = LinearRegression().fit(fama_french_factors_including_mkt_premium,portfolio_returns)
