# multiple regression with prob 13.1 data
import numpy as np
import statsmodels.api as sm

y = [12,18,31,28,52,47,38,22,36,17]

x = [
     [174,281,189,202,149,188,215,150,167,135],
     [3,9,4,8,9,12,5,11,8,5]
     ]

def reg_m(y, x):
    x = np.array(x).T
    x = sm.add_constant(x)
    results = sm.OLS(endog=y, exog=x).fit()
    return results

print reg_m(y,x).summary()