Проблема вычитает две гистограммы и вычисляет CDF в Python
У меня проблема с вычитанием двух гистограмм и вычислением их CD
F и его Kolmogorov Smirnov Test
. У меня есть данные двух файлов, один из которых самый большой, и у меня было масштабирование. Затем складываю эти файлы в гистограммы, пытаюсь вычесть и вычислитьCDF
, но я не могу.
input_file1 = "F:/JavierRengifo/Datos/muon_counter_data_Test_Th15_Tamb_LuzNormal.txt" input_file2 = "F:/JavierRengifo/Datos/muon_counter_data_Sr01uCi_dist05cm_Th15_Det1_NM.txt" input_file3 = "F:/JavierRengifo/Datos/muon_counter_data_Ce025uCi_dist05cm_Th15_Det1_NM.txt"
df1 = pd.read_table(input_file1,delimiter = '', skiprows = 5, header = 0)
df2 = pd.read_table(input_file2,delimiter = '', skiprows = 5, header = 0)
df3 = pd.read_table(input_file3,delimiter = '', skiprows = 5, header = 0)
arr1 = df1 ["SiPM[мВ]"]
arr2 = df2 ["SiPM[мВ]"]
arr3 = df3 ["SiPM[мВ]"]
BINS = 100
(количество, ячейки) = np.histogram(arr1, bin=BINS)
коэффициент = 1/6
arr1n = plt.hist(интервалы [:-1], интервалы, веса = фактор * количество, кумулятивное значение =False, плотность = True, color = 'blue', alpha = 0.75, histtype = 'step')
arr2n = plt.hist(arr2, bins = BINS, cumulative=False, density = True, color = 'red', alpha = 0.75, histtype = 'step')
arr3n = plt.hist(arr3, bins = BINS, cumulative=False, density = True, color = 'orange', alpha = 0.75, histtype = 'step')
dfAB = arr2n[0]-arr1n[0]
dfAC = arr3n[0]-arr1n[0]
lala1 = plt.hist(diffABv0, бины = 100, кумулятивное =True, density=True, color='red', alpha=0,75, histtype = 'step')
lala2 = plt.hist(diffACv0, бины = 100, кумулятивное = True, density = True, color = 'purple', alpha = 0.75, histtype = 'step')
print ('1 намерение Ks: '+str(ks_2samp(dfAB,dfAC)))