В документации Сципи говорится:
Формула, используемая для теоретических квантилей (горизонтальная ось вероятностного графика), является оценкой Филлибена:
квантили = dist.ppf (val), для 0,5 ** (1 / n), для i = n
val = (i - 0,3175) / (n + 0,365), для i = 2, ..., n-1
1 - 0,5 ** (1 / n), для i = 1
, где i обозначает i-е упорядоченное значение, а n - общее количество значений.
Но когда я подключаю к нему некоторые числа, он дает значения, отличные от значений, вычисленных функцией probplot. Например:
np.random.seed(12)
n = np.random.randn(10) * 10
x, _ = stat.probplot(n, plot = plt)
-15.3472134 - это второе значение, поэтому я подключаю
(2 - 0,3175) / (10 + 0,365) и возвращает квантиль 0,162325 (число между 0и 1) тогда как функция отображает ее как -0,98494667 и помечает теоретические квантили оси X. Как вычисляется теоретический квантиль -0,98494667 из фактического квантиля 0,162325?