正規分布確率 計算機|任意区間の確率・パーセンタイル・両側/片側の一気算出
正規分布N(μ,σ²)の任意区間 P(a
正規分布確率ツール
計算式(PDFとCDF)
確率密度関数(PDF)
f(x) = (1 / (σ√(2π))) × exp( -(x-μ)² / (2σ²) )
μ(平均)で釣鐘型の中心、σ(標準偏差)で広がりを決める2パラメータの連続分布。μ±σの範囲に約68%、μ±2σに約95%、μ±3σに約99.7%が入るのが特徴です。
累積分布関数(CDF)
Φ(z) = (1/2) × (1 + erf(z/√2))
P(a < X < b) = Φ((b-μ)/σ) - Φ((a-μ)/σ)
正規分布の累積確率は誤差関数erf(x)で表現され、閉形式では書けないため数値計算します。本ツールは Abramowitz & Stegun のerf近似式(相対誤差1.5×10⁻⁷)で計算しています。
Z値への変換と標準正規分布
Z = (X - μ) / σ
Z値は「観測値が平均から何σ離れているか」を表す標準化指標です。Z=1は平均+1σ、Z=-2は平均-2σ。標準正規分布 N(0,1) にすると、平均・標準偏差の違うデータを同じスケールで比較可能。
標準化の意義
- 単位・尺度の異なるデータの直接比較(Z-score standardization)。
- 統計検定量として仮説検定に使用(|Z|>1.96で両側5%有意)。
- 異常値検出(|Z|>3で外れ値候補)。
- 機械学習の前処理(Standard Scaler)でモデル収束を高速化。
68-95-99.7ルール(経験則)
正規分布の最も重要な性質。試験得点分布・身長・測定誤差など多くの自然現象がこのルールに従います。
| 区間 | 含まれる確率 | 外側(両側) | 外側(片側) |
|---|---|---|---|
| μ ± 1σ | 68.27% | 31.73% | 15.87% |
| μ ± 1.645σ | 90.00% | 10.00% | 5.00% |
| μ ± 1.96σ | 95.00% | 5.00% | 2.50% |
| μ ± 2σ | 95.45% | 4.55% | 2.28% |
| μ ± 2.576σ | 99.00% | 1.00% | 0.50% |
| μ ± 3σ | 99.73% | 0.27% | 0.135% |
| μ ± 3.29σ | 99.90% | 0.10% | 0.05% |
| μ ± 4σ | 99.9937% | 0.006% | 0.003% |
| μ ± 6σ | 99.9999998% | 2×10⁻⁹ | 1×10⁻⁹ |
製造業のシックスシグマ品質管理では、工程能力Cpk=2.0(=6σ距離)で不良率3.4ppmを目標とします(実際は1.5σのシフトを想定した長期的評価)。
Z値-確率対応早見表
片側確率 P(Z ≤ z) の早見表(標準正規分布 N(0,1))。
| Z | Φ(Z) | Z | Φ(Z) | Z | Φ(Z) |
|---|---|---|---|---|---|
| -3.0 | 0.0013 | -1.0 | 0.1587 | 1.0 | 0.8413 |
| -2.58 | 0.0049 | -0.5 | 0.3085 | 1.28 | 0.9000 |
| -2.33 | 0.0099 | -0.25 | 0.4013 | 1.645 | 0.9500 |
| -1.96 | 0.0250 | 0.0 | 0.5000 | 1.96 | 0.9750 |
| -1.645 | 0.0500 | 0.25 | 0.5987 | 2.33 | 0.9901 |
| -1.28 | 0.1000 | 0.5 | 0.6915 | 2.58 | 0.9951 |
パーセンタイル逆算表
「上位10%は何σ以上か」等、確率からZ値を逆算した表。
| パーセンタイル | Z値 | 意味 |
|---|---|---|
| 50%(中央値) | 0.000 | 平均と同じ |
| 60% | 0.253 | 上位40%の閾値 |
| 75% | 0.674 | 四分位Q3 |
| 80% | 0.842 | 上位20%の閾値 |
| 90% | 1.282 | 上位10%の閾値 |
| 95% | 1.645 | 片側5%(有意水準) |
| 97.5% | 1.960 | 両側5%の片側閾値 |
| 99% | 2.326 | 上位1% |
| 99.5% | 2.576 | 両側1%の片側閾値 |
| 99.9% | 3.090 | 上位0.1% |
Excel・Python・Rでの計算方法
Excel/Google Sheets
- =NORM.DIST(x, μ, σ, TRUE) ― P(X ≤ x)、累積分布。
- =NORM.DIST(x, μ, σ, FALSE) ― 確率密度 f(x)。
- =NORM.INV(prob, μ, σ) ― パーセンタイル逆算(分位点)。
- =NORM.S.DIST(z, TRUE) ― 標準正規分布のCDF Φ(z)。
- =NORM.S.INV(prob) ― 標準正規分布の分位点。
Python (scipy)
from scipy.stats import norm # P(a < X < b) for N(mu, sigma) prob = norm.cdf(b, mu, sigma) - norm.cdf(a, mu, sigma) # 95%パーセンタイル q95 = norm.ppf(0.95, mu, sigma) # 標準正規のPDF density = norm.pdf(0) # 0.3989
R
# P(a < X < b) prob <- pnorm(b, mu, sigma) - pnorm(a, mu, sigma) # 分位点 q95 <- qnorm(0.95, mu, sigma) # 乱数生成 rnorm(1000, mean=0, sd=1)
実データの正規性確認
「正規分布に従うと仮定して〜」の前に、実データの正規性を確認する必要があります。
可視化による確認
- ヒストグラム:釣鐘型か目視確認。
- Q-Qプロット:実データ分位点と理論正規分位点を比較。直線に乗れば正規。
- 箱ひげ図:中央値の位置、対称性、外れ値を確認。
統計的検定
- Shapiro-Wilk検定:n≤5,000で最も検出力が高い。
- Kolmogorov-Smirnov検定:任意分布との比較にも使えるが検出力低め。
- Anderson-Darling検定:裾の逸脱に敏感。品質管理で頻用。
- Jarque-Bera検定:歪度・尖度から正規性判定。計量経済で標準的。
非正規時の対応
- ログ変換・Box-Cox変換で正規化を試みる。
- ノンパラメトリック検定(Mann-Whitney U、Kruskal-Wallis)に切替。
- ブートストラップで経験分布からp値を推定。
- ロバスト統計(M推定量、トリム平均)を使用。
業界・場面別の使い方
品質管理・製造業
工程能力指数Cp・Cpk、シックスシグマの不良率算定、SPC管理図(Shewhart Chart)のUCL/LCL設計。JIS Z 8101品質管理用語、ISO 9001品質マネジメントの標準的な統計解析手法として組み込まれています。
統計的仮説検定
Z検定・t検定・カイ二乗検定の基礎。信頼区間の設定、A/Bテストの有意性判定、サンプルサイズ設計に必須。医学・心理学・社会科学すべての実証研究の土台。
金融・リスク管理
VaR(Value at Risk)計算、ポートフォリオリスク評価、ブラック・ショールズモデル。ただし金融データはファットテールを示すことが多く、正規分布での過小評価は金融危機の一因になりました。t分布・スチューデントt・EVT(極値理論)との併用が主流。
教育・心理測定
偏差値(平均50・SD10の正規変換)、IQスコア(平均100・SD15)、TOEIC・SAT等の標準化テスト。個人得点の相対位置をパーセンタイルで表示し、母集団内の位置づけを可視化。
医療・臨床検査
血液検査・生理検査の基準範囲(平均±2SD=95%区間)、疾患スクリーニングの感度・特異度、ROC曲線分析。健康診断の基準値算定は多くが正規分布近似で設計されています。
機械学習・データサイエンス
特徴量の標準化(Standard Scaler)、異常検知(Z-score>3の外れ値検出)、Gaussian Mixture Model、ベイズ推論の事前分布。深層学習の重み初期化(Glorot/He initialization)も正規分布ベース。
よくある間違い・注意点
- 正規分布と決めつける — 実データの多くは非正規です。ヒストグラム・Q-Qプロット・正規性検定で必ず確認してから正規モデルを適用してください。
- ファットテール(重い尾)を無視する — 金融・保険・地震規模など「稀な大事象」が起きる分野では、正規分布は極端値の発生確率を大幅に過小評価します。t分布や一般化パレート分布への切替を検討。
- 平均±2SDで95%区間と表現 — 正確には95.45%です。厳密な95%区間は±1.96SDで、5%有意水準の検定と整合します。
- 片側検定と両側検定を混同 — 「上限」「下限」のみ検討する場合と、両側どちらでも検討する場合でZ値の閾値が異なります。目的を事前に明確化してください。
- Z=3で外れ値と即断 — サンプルサイズが大きければ、|Z|>3のデータは正規分布でも自然発生します(n=1000で約3個)。IQRやマハラノビス距離との併用が実務。
- 中心極限定理の乱用 — 「サンプル平均は正規分布に近づく」は真ですが、元データの分布形状・サンプル数・独立性の前提が必要です。歪んだ分布・依存データでは収束が遅い。
- 離散変数を無理に正規近似 — 二項分布・ポアソン分布は、パラメータが小さい範囲では正規近似の誤差が大きくなります。連続修正やexact法との併用を検討。
関連する規格・標準
- JIS Z 8101 ― 品質管理用語。正規分布・工程能力・管理限界の日本語定義。
- JIS Z 9020 ― シューハート管理図。正規分布に基づくSPCの実施手順。
- ISO 3534 ― 統計 - 用語及び記号。国際的な統計用語の標準。
- ISO 16269 ― データの統計的解釈。区間推定・分位点推定の標準手順。
- ISO 5479 ― 正規性からの逸脱の検定。Shapiro-Wilk、Kolmogorov-Smirnov等の標準実装。
- ASTM E2554 ― 統計的品質管理の実施ガイド。
- NIST SP 800-90A ― 乱数生成の標準。正規分布乱数の暗号学的用途。
参考文献・公的資料
より深く学ぶには、以下の信頼できる公的機関・研究機関の資料を参照してください。
- 総務省 統計局 ― 統計調査の設計・解析、日本の公的統計データ。
- 日本統計学会 ― 統計学教育・研究の中核学会。
- NIST/SEMATECH e-Handbook of Statistical Methods ― 統計工学の標準ハンドブック。正規分布章に詳細実例あり。
- SciPy - scipy.stats.norm ― Pythonの標準的な正規分布実装。CDF/PDF/PPF/rvsの完全ドキュメント。
- R Project - stats::Normal ― Rの正規分布関数(dnorm/pnorm/qnorm/rnorm)ドキュメント。
- 日本産業標準調査会(JISC) ― JIS Z 8101、Z 9020等の統計・品質管理JIS規格。
- ISO/TC 69 統計手法の応用 ― 統計手法の国際規格を策定する専門委員会。
- 日本品質管理学会(JSQC) ― 品質管理・シックスシグマ・SPCの実務ノウハウ。
- JASP 統計ソフトウェア ― フリーの統計ソフト。正規分布可視化とベイズ解析が容易。
- 日本行動計量学会 ― 心理・社会科学分野の統計手法学会。
よくある質問(FAQ)
Z=1.96の意味は?
±1.96SDで両側95%カバー。統計的仮説検定の5%有意水準の閾値で、信頼区間の幅の根拠。片側5%なら1.645、両側1%は2.576、両側0.1%は3.29が対応するZ値です。
片側検定と両側検定はどう選ぶ?
「効果があるかどうか(方向は問わない)」なら両側検定、「明確に増える/減ることのみ検討」なら片側検定。実務では基本的に両側検定を使い、片側は事前登録された明確な仮説がある場合に限定します。
正規性検定はどうやる?
Shapiro-Wilk検定(scipy.stats.shapiro、R shapiro.test)が最も検出力が高くn≤5000で推奨。他にAnderson-Darling(裾に敏感)、Kolmogorov-Smirnov、Jarque-Bera(歪度・尖度から判定)などがあります。可視化としてQ-Qプロットも併用。
ファットテールとは何?
正規分布より裾が重い分布のこと。金融リターン・地震規模・保険損失は極端値の発生確率が正規分布の予想より高く、t分布・パレート分布・一般化極値分布(GEV)でモデル化するのが実務。
68-95-99.7ルールはどこまで信用できる?
純粋な正規分布に対する厳密値です。実データが正規分布に近ければ良い目安になりますが、非正規なら大きく外れます。品質管理では長期的な工程シフト1.5σを想定するため、Cpk=2で「6σレベル」でも実質不良率は3.4ppmです。
Excelで確率を計算する関数は?
=NORM.DIST(x, μ, σ, TRUE)で累積確率、=NORM.INV(prob, μ, σ)でパーセンタイル逆算。標準正規なら=NORM.S.DIST(z, TRUE)、=NORM.S.INV(prob)を使います。Google Sheetsでも同一の関数が使えます。
Z-scoreを使った異常検知の閾値は?
慣習的に|Z|>3を外れ値候補とします。ただし大量データでは正規分布でも自然発生(n=1000で約3個)するため、IQRベース(Q1-1.5IQR〜Q3+1.5IQR)、マハラノビス距離、Isolation Forestなどとの併用が実務。
偏差値はどう計算する?
Z値を平均50・SD10に変換:偏差値 = 50 + 10 × (X-μ)/σ = 50 + 10×Z。偏差値60は上位約16%、70は上位2.3%、80は上位0.13%。母集団の分布が正規に近い前提で意味を持ちます。
中心極限定理はどんな条件で成立?
独立同分布のサンプル平均は、母分布に依存せずn→∞でN(μ, σ²/n)に収束。実務上はn≥30で正規近似OKと言われますが、母分布が極端に歪んでいる場合はn=100以上でも収束が不十分な場合があります。
t分布と正規分布の使い分けは?
母分散σ²が既知→正規分布のZ検定、未知(標本SDで代用)→t分布のt検定。t分布は自由度が小さい(n<30)ほど裾が重く、正規分布に比べて広い信頼区間となります。n≥30でほぼ同一結果。
正規分布の乱数はどう生成する?
Box-Muller変換、Marsaglia極方法、Ziggurat法などが標準アルゴリズム。Python:numpy.random.normal、R:rnorm、Excel:=NORM.INV(RAND(), μ, σ)。シミュレーション・機械学習の初期化に頻用。
対数正規分布との違いは?
Xがlog正規なら log(X) が正規分布。所得・株価・生存時間など正の値をとる歪んだ分布のモデル化に使います。正規と対数正規は密接に関連し、対数変換で正規化が図れます。