3 $\sigma$ kritérium

import numpy as np
import matplotlib.pyplot as plt

#naměřená data
data=np.array([1.656,1.589,1.719,8.188,1.648,1.636,1.711,1.642,1.597,1.675,1.690,-8.147,1.669,1.652,1.514,1.499,1.660,1.678,1.669,1.618,1.617,1.559,1.142,1.014,1.403,1.656,1.546,1.619,1.658,1.647])

ndata=np.size(data)
print(f"počet dat {ndata}")

počet dat 30

mean=np.mean(data) #aritmetický průměr
stdev=np.std(data,ddof=1) #standardní odchylka (nepředpojatý odhad)
print(f"aritmetický průměr : {mean:8.5f}")   
print(f"standardní odchylka: {stdev:8.5f} (nepředpojatý odhad)")   

aritmetický průměr :  1.48080
standardní odchylka:  2.18675 (nepředpojatý odhad)

n_hruba_chyba=0
for i in range(ndata):
    if abs(data[i]-mean)>3*stdev:
        print(f"hrubá chyba: {i+1}. hodnota {data[i]}") #pokud se hodnota liší od aritmetického průměru o více než 3 sigma vypíšeme ji
        n_hruba_chyba+=1
print(f"počet hrubých chyb: {n_hruba_chyba}")   

hrubá chyba: 4. hodnota 8.188
hrubá chyba: 12. hodnota -8.147
počet hrubých chyb: 2