Многомерная полиномиальная регрессия с numpy

У меня есть много выборок (y_i, (a_i, b_i, c_i)), где y Предполагается, что изменяется как многочлен от a,b,cдо определенной степени. Например, для заданного набора данных и степени 2 я мог бы создать модель

y = a^2 + 2ab - 3cb + c^2 +.5ac

. Это можно сделать с помощью метода наименьших квадратов, и это небольшое расширение процедуры полифита numpy. Есть ли стандартная реализация где-нибудь в экосистеме Python?