Probability Distributions

import numpy as np
from scipy import stats
np.random.seed(0)
a=np.random.normal(70,10,40); b=np.random.normal(76,10,40)
print("Mean A=%.1f  B=%.1f"%(a.mean(),b.mean()))
t,p=stats.ttest_ind(a,b); print(f"t={t:.2f}, p={p:.4f} ->", "significant" if p<0.05 else "not significant")
ci=stats.t.interval(0.95,len(a)-1,loc=a.mean(),scale=stats.sem(a)); print("95% CI of A:",np.round(ci,1))