gistfile1.py

from sklearn import linear_model

X = [
  [4,3,1,0,1],
  [5,2,1,0,1],
  [4,2,1,1,1],
  [3,1,0,1,1],
  [1,1,0,1,1]
]

y = [4, 6, 6, 3, 1]

clf = linear_model.LinearRegression()
clf.fit (X, y)

print clf.coef_
# [ 1. -1.  3.  1.  0.]