print("Model coefficients: ") for i in range(X.shape[1]): print(X.columns[i], "=", model.coef_[i].round(4))