CDF в scipy pearson3, по-видимому, не определен

Я пытаюсь построить распределение, используя scipy.stats.pearson3, но я получаю противоречивые результаты при просмотре кумулятивной функции распределения (или функции выживания).

В качестве примера:

scipy.stats.pearson3.cdf(-1, skew=0, loc=0, scale=1) 
>> 0.15865525393145707

scipy.stats.pearson3.cdf(-1, skew=-0.0001, loc=0, scale=1) 
>> 0.8413447461693605

scipy.stats.pearson3.cdf(+1, skew=-0.0001, loc=0, scale=1) 
>> 0.1586552538306324

Первый результат в порядке; например, это то, что вы должны ожидать от skew = 0 pearson3 (тот же результат, что и при нормальном распределении).

Но второй и третий результаты противоречивы: они показывают уменьшающуюся CDF, что не имеет смысла (Этого не происходит, если перекос положительный). Кажется, что распределение pearson3 с отрицательной асимметрией имеет инвертированные функции CDF и SF. Как ни странно, PDF-файл кажется правильным, что приводит к еще большей путанице.

Кто-нибудь знает, является ли это известной проблемой с scipy или я что-то упустил?


person ecortazar    schedule 05.12.2019    source источник
comment
Это ошибка. Я создал задачу для этого в репозитории SciPy на github: github.com/scipy/scipy/issues/11186   -  person Warren Weckesser    schedule 07.12.2019


Ответы (1)


Попробуйте построить график 1 - pearson3.cdf для данных с отрицательным перекосом.

CDF = F(x) для распределения с положительной асимметрией (коэффициент > 0)

CDF = 1 - F (x) для распределения с отрицательной асимметрией (перекос ‹ 0)

Ссылка: Рао Д.В. (2009). Обсуждение «Распределение типа Лог-Пирсона и его применение в анализе повторяемости паводков». I: Характеристики распределения» В. В. Гриффиса и Дж. Р. Стедингера. Journal of Hydrologic Engineering, 14(3), 301. https://doi.org/10.1061/(ASCE)1084-0699(2009)14:3(301)

person Jane Harrell    schedule 28.03.2020