Scipy check normal distribution

from scipy import stats
pts = 1000
np.random.seed(28041990)
a = np.random.normal(0, 1, size=pts)
b = np.random.normal(2, 1, size=pts)
x = np.concatenate((a, b))

k2, p_value = stats.normaltest(x)

alpha = 1e-3
print("p_value = {:g}".format(p))
# p_value = 3.27207e-11

if p_value < alpha:  # null hypothesis: x comes from a normal distribution
  print("The null hypothesis can be rejected")
else:
  print("The null hypothesis cannot be rejected")
# The null hypothesis can be rejected

Source: docs.scipy.org

Similar codes

45 votes100 views10 answers
Youtube
Python
0 vote105 views1 answer
Python boxplot show mean
Python
5 votes110 views2 answers
How to plot pie chart in python
Python
0 vote114 views1 answer
AttributeError: module 'django.db.models' has no attribute 'ArrayField'
Python
1 vote105 views2 answers
Python dictionary default value if key doesn't exist
Python
-1 vote94 views1 answer
Find the area of a circle in python
Python
1 vote141 views1 answer
Passing user instance django form submission
Python Django
6 votes102 views3 answers
Flask cors
Python
-2 votes105 views2 answers
String pythhon
Python
1 vote112 views1 answer
Sleep 5 seconds py
Python