python - 标记曲线之间的最大偏差/距离

Question

我只想使用 matplotlib 标记两条曲线之间发生最大偏差的位置。请帮助我。

垂直距离适用于 Kolmogorov–Smirnov test

import numpy as np 
%matplotlib inline
import matplotlib.pyplot as plt
import scipy.stats as stats
#------------------------------------
data=np.random.uniform(low=1,high=10,size=300)
standardized_data=np.sort(data-np.mean(data))/np.std(data)
probs=np.arange(1.0, 301)/300
plt.plot(standardized_data,probs) #curve1
plt.plot(stats.norm.ppf(probs),probs) #curve2
plt.show()

输出:

像这样

或者像这样

Answer 1

您需要将其中一条曲线插值到另一条曲线的 x 值上。然后找出它们的差值的最大值。

import numpy as np 
import matplotlib.pyplot as plt
import scipy.stats as stats
#------------------------------------
data=np.random.uniform(low=1,high=10,size=300)
x1= np.sort(data-np.mean(data))/np.std(data)
y = np.arange(1.0, 301)/300

x2 = stats.norm.ppf(y)
yc = np.interp(x1, x2, y)

ind_max = np.argmax((yc-y)**2)


plt.plot(x1, y) #curve1
plt.plot(x2, y) #curve2

plt.axvline(x1[ind_max], color="red", linestyle="dashed", alpha=0.4)
plt.plot([x1[ind_max], x1[ind_max]], [y[ind_max], yc[ind_max]], color="red")
plt.show()