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

71 votes235 views5 answers
How to make loop python
Python
1 vote277 views1 answer
Can you print to multiple output files python
Python
1 vote286 views1 answer
Error popup in django not visible
Python Django
2 votes248 views1 answer
Boto signed url
Python
17 votes253 views2 answers
How to save array python
Python
3 votes324 views1 answer
Python delete saved image
Python
0 vote280 views1 answer
Sort for Linked Lists python
Python
0 vote267 views2 answers
Concatenate the squares of numbers in python
Python
0 vote250 views1 answer
Hazm
Python
0 vote241 views1 answer
Exclude serializer
Python