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

13 votes212 views3 answers
Permutation python
Python
-1 vote250 views1 answer
Django cheat sheet pdf
Python Django
5 votes254 views1 answer
How to override save method in django
Python Django
9 votes229 views2 answers
Matplotlib axis rotate xticks
Python
0 vote238 views1 answer
Python string: index error
Python
0 vote222 views1 answer
Matplotlib colorbar reverse direction
Python
1 vote214 views1 answer
Find the value counts for the column 'your_column'
Python
1 vote224 views1 answer
Python moving average of list
Python
0 vote227 views1 answer
Queue using linked list in python
Python
0 vote264 views1 answer
Cant import flask mail
Python Flask