Блог

Всякий раз, когда я пытаюсь вычислить логическую вероятность распределения Tweedie, я получаю только nans. Следующий пример имеет мощность tweedie p = 2 (которая является гамма-распределением) и вывод гаммы в качестве сравнения:

 import statsmodels as sm

aux_gamma = sm.genmod.families.Gamma()
print(aux_gamma.loglike(1000, 10))

aux_tweedie = sm.genmod.families.Tweedie(var_power=2)
print(aux_tweedie.loglike(1000, 10))
  

Из документации:

Функция логической вероятности не реализована из-за сложности вычисления бесконечной серии суммирования. Степень отклонения может быть оценена с помощью функции estimate_tweedie_power, которая является частью statsmodels.genmod.generalized_linear_model.Класс GLM.

Единственная функция правдоподобия, доступная в этом случае, — это расширенная функция квази-правдоподобия, предложенная Нелдером и Прегибоном (1987). Его можно использовать, установив eql=True внутри .Tweedie() функции.