確率統計の勉強と同時にnumpyとmatplotlibの勉強も進めようと思い一般正規分布の描画をしてみました。
python
1import numpy as np 2import matplotlib.pyplot as plt 3 4np.random.seed(1000) 5 6μ = 5 # 平均値 7σ = 3 # 標準偏差 8N = 100000 # データの個数 9 10r = np.random.normal(μ, σ, N) # 平均μ,標準偏差σの正規分布をN個生成 11r2 = r[np.where((r <= µ+1.96*σ)&(r >= µ-1.96*σ))] # rの95%の範囲 12# print(len(r2) / len(r) * 100) # 94.979% 13 14plt.hist(r, bins=100, color='#ff0000', alpha=0.5) 15plt.hist(r2, bins=100, color='#000000', alpha=0.5) 16plt.show()
赤色の一般正規分布の95%に当たる部分を黒のヒストグラムで重ねようと思ったのですが、横軸についてはちゃんとできているのですが、縦軸についてが上手く行っていません。
コードを見たところ問題無いように思うのですが、どうして95%の部分の上部が黒で描画されていないのでしょうか。
5-1.963 <= 横軸 <= 5+1.963
の範囲は全部黒になると思うのですが、上手く行っていません。
どなたかわかる方教えていただけませんでしょうか。
よろしくお願いします。
回答1件
あなたの回答
tips
プレビュー
バッドをするには、ログインかつ
こちらの条件を満たす必要があります。
退会済みユーザー
2017/07/18 15:55
2017/07/18 16:35 編集
退会済みユーザー
2017/07/19 02:30
2017/07/19 03:12