Xが3次元

3x.py

import numpy as np
from sklearn import linear_model
import matplotlib.pyplot as plt

X = np.array([[0, 0, 1], [1, 1, 2], [2, 2, 3], [3, 3, 3.1]])
Y = np.array([0.1, 1.1, 1.8, 2.7])

regr = linear_model.LinearRegression()
regr.fit(X, Y)

print 'Coefficients:',   regr.coef_
print 'Residual sum of squares: %f' % np.mean((regr.predict(X) - Y) ** 2)
print 'Variance score %f' % regr.score(X, Y)