Я пытаюсь подобрать распределение Пуассона, чтобы вычислить MLE для моих данных. Я заметил, что в статистике scipy есть .fit для непрерывных функций, но для дискретных функций нет. Есть ли другой API, который имеет функцию .fit для дискретных распределений в Python?