python import statsmodels.api as sm import numpy as np np.random.seed(0) X = np.random.normal(0, 1, size=(100, 2)) y = 2 * X[:, 0] + 3 * X[:, 1] + 0.5 * np.random.normal(0, 1, size=100) X = sm.add_constant(X) model = sm.OLS(y, X) results = model.fit() print(results.summary())