確率分布のPythonへの実装
メニューを表示するにはスワイプしてください
二項分布
二項分布は、各試行の成功確率が p である n 回の独立した試行において、ちょうど k 回成功する確率をモデル化。
123456789101112131415161718from scipy.stats import binom import matplotlib.pyplot as plt # number of trials n = 100 # probability of success p = 0.02 # number of successes k = 3 binom_prob = binom.pmf(k, n, p) # Vizualization x_vals = range(0, 15) y_vals = binom.pmf(x_vals, n, p) plt.bar(x_vals, y_vals, color='skyblue') plt.title(f'Binomial probability: {binom_prob:.4f}') plt.show()
n = 100- 100本の棒を検査;p = 0.02- 棒が不良品である確率は2%;k = 3- ちょうど3本が不良品である確率;binom.pmf()は確率質量関数を計算。
一様分布
一様分布は、$a$ から $b$ までのすべての値が等しく発生する連続変数をモデル化。
1234567891011121314151617from scipy.stats import uniform import matplotlib.pyplot as plt import numpy as np a = 49.5 b = 50.5 low, high = 49.8, 50.2 uniform_prob = uniform.cdf(high, a, b - a) - uniform.cdf(low, a, b - a) # Vizualization x = np.linspace(a, b, 100) pdf = uniform.pdf(x, a, b - a) plt.plot(x, pdf, color='black') plt.fill_between(x, pdf, where=(x >= low) & (x <= high), color='lightgreen', alpha=0.5) plt.title(f'Uniform probability: {uniform_prob:.1f}') plt.show()
a, b- 棒の長さの全範囲;low, high- 注目する区間;- CDF値の差分が区間内の確率。
正規分布
正規分布は、平均値 $\mu$ の周囲に値が集中し、標準偏差 $\sigma$ によって広がりが測定される分布。
1234567891011121314151617181920import numpy as np import matplotlib.pyplot as plt from scipy.stats import norm mu = 200 sigma = 5 lower, upper = 195, 205 norm_prob = norm.cdf(upper, mu, sigma) - norm.cdf(lower, mu, sigma) z1 = (lower - mu) / sigma z2 = (upper - mu) / sigma # Vizualization x = np.linspace(mu - 4*sigma, mu + 4*sigma, 200) pdf = norm.pdf(x, mu, sigma) plt.plot(x, pdf, color='black') plt.fill_between(x, pdf, where=(x >= lower) & (x <= upper), color='plum', alpha=0.5) plt.title(f'Normal probability: {norm_prob:.4f}\nZ-scores: {z1}, {z2}') plt.show()
mu- 平均ロッド重量sigma- 標準偏差- 確率 - CDF の差分
- Zスコアは、境界が平均からどれだけ離れているかを示す
実世界での応用
- 二項分布 - 不良な棒が特定の本数である確率
- 一様分布 - 棒の長さが許容範囲内かどうか
- 正規分布 - 棒の重さが期待されるばらつき内かどうか
これらを組み合わせることで、品質管理は欠陥の特定、精度の確保、製品の一貫性の維持を実現。
すべて明確でしたか?
フィードバックありがとうございます!
セクション 5. 章 11
AIに質問する
AIに質問する
何でも質問するか、提案された質問の1つを試してチャットを始めてください
セクション 5. 章 11