python import numpy as np import pandas as pd import statsmodels.api as sm data = pd.read_csv('data.csv') X = sm.add_constant(X) model = sm.OLS(y, X) results = model.fit() print(results.summary())