Distribution of Regression Residuals: Is this a normal distribution?
You can perform a statistical test to confirm your data is normally distributed
Try:
from scipy import stats

np.random.seed(42)
x = np.random.normal(2, 1, size=1000)
k2, p = stats.normaltest(x)
alpha = 0.001
print(&quot;p = {:g}&quot;.format(p))

if p &lt; alpha:  # null hypothesis: x comes from a normal distribution
    print(&quot;The null hypothesis can be rejected&quot;)
else:
    print(&quot;The null hypothesis cannot be rejected&quot;)


This function tests the null hypothesis that a sample comes from a normal distribution. It is based on D’Agostino and Pearson’s test that combines skew and kurtosis to produce an omnibus test of normality.