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

11 votes143 views5 answers
Python Ordered Dictionary
Python
0 vote147 views4 answers
Python bool
Python
2 votes205 views1 answer
Dataframe from two series
Python
1 vote105 views1 answer
Beautifulsoup remove element
Python
0 vote154 views1 answer
Using shebang python
Python
0 vote151 views1 answer
Unreadable Notebook: jupyter
Python
3 votes173 views2 answers
Np load csv
Python
0 vote162 views1 answer
How to separate string in python by blank line
Python
1 vote136 views2 answers
Remove 1st column pandas
Python
4 votes145 views2 answers
Pandas most frequent value
Python