predict = [0.49, 0.18]
poly = PolynomialFeatures(degree=2)
X_ = poly.fit_transform(X)
predict_ = poly.fit_transform(predict)
X_ = delete(X_, 1, axis=1)
predict_ = delete(predict_, 1, axis=1)
clf = linear_model.LinearRegression()
clf.fit(X_, vector)
print(('X_ = ', X_))
print(('predict_ = ', predict_))
print(('Prediction = ', clf.predict(predict_)))
Multivariate polynomial regression with numpy