Как в Python лучше всего подогнать квадратичный многочлен к p-мерным данным и вычислить его градиент и гессиан? ⇐ Python

[Anonymous]

Я пытался использовать библиотеку scikit-learn для решения этой проблемы. Примерно:

Код: Выделить всё

from sklearn.preprocessing import PolynomialFeatures
from sklearn.linear_model import LinearRegression

# Make or load an n x p data matrix X and n x 1 array y of the corresponding
# function values.

poly = PolynomialFeatures(degree=2)
Xp = poly.fit_transform(X)
model = LinearRegression()
model.fit(Xp, y)

# Approximate the derivatives of the gradient and Hessian using the relevant
# finite differences equations and model.predict.

Как показано выше, sklearn делает выбор в пользу разделения полиномиальной регрессии на PolynomialFeatures и LinearReгрессия, а не объединения их в одну функцию. Такое разделение имеет концептуальные преимущества, но также имеет и серьезный недостаток: оно фактически не позволяет модели предлагать атрибутыgradient_ и hessian_, а модель была бы значительно полезнее, если бы она это делала. .
Как указано выше, мой текущий обходной путь заключается в использовании уравнений конечных разностей и model.predict для аппроксимации элементов градиента и гессиана (как описано здесь). . Но мне не нравится такой подход. Есть ли более элегантный и точный метод в Python? Мне подойдет тот, который использует другую библиотеку.

Подробнее здесь: https://stackoverflow.com/questions/790 ... ensional-d