Как дифференцировать функцию, суммирующую две матрицы? ⇐ Python

[Anonymous]

Я пытаюсь реализовать собственную функцию отрицательных биномиальных потерь для моей модели LGBM. Вот мой код:

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

def negative_binomial_objective(y_true, y_pred):
p = 0.5

def loss(x, t):
loss = (
gamma(x) + gamma(t + 1) - gamma(x + t) - x * np.log(p) - t * np.log(1 - p)
)
return loss

def compute_grad(x):
return derivative(f=lambda x: loss(x, y_true), x=x)

grad = compute_grad(x=y_pred)
hess = derivative(f=compute_grad, x=y_pred)  # second derivative

return grad, hess

Проблема в том, что производная функция scipy расширяет x в соответствии с параметром order (по умолчанию установлено значение 8). Затем я получаю сообщение об ошибке, что мои операнды не могут быть переданы вместе (например, в пределах гаммы (x + t)), поскольку они имеют формы (N,8) и (N,).
Как мне найти производную этой функции?

Подробнее здесь: https://stackoverflow.com/questions/798 ... o-matrices