公開日: 04/01/2021

工程能力分析の統計的詳細

ここでは、「工程能力分析」レポートの統計量の計算方法について説明します。

変動の統計量

いくつかのオプションがありますが、正規分布に基づく工程能力分析では、基本的に同じ計算式が使用されています。オプションは、シグマ(σ)の計算方法が異なるだけです。

σ(全データから推定)

全体シグマを使用します。このオプションはPpk統計量に使用され、次式によりシグマが算出されます。

Equation shown here

注: 「長期シグマ(全データから推定)」レポートでは、他のシグマのレポートと同様に、工程能力指数のラベルに「Cp」がデフォルトで使用されています。このレポートで「PPK」のラベルを使用するには、[ファイル]>[環境設定]>[プラットフォーム]>[一変量の分布]>[Ppk工程能力ラベル設定]を選択します。

管理図シグマ

管理図の設定値により決定されるシグマを使用します。

シグマの値を、管理図起動ウィンドウにある[統計量の指定]ボタンを用いて指定すると、指定した値が工程能力指数の計算に使用されます。

[移動範囲(平均)]オプションを使用するIR管理図では、シグマの値は次のように計算されます。Equation shown here

ここで

Equation shown here は、移動範囲の平均

d2(n)は、母標準偏差が1であるn個の独立した正規分布に従う確率変数の範囲の期待値、ここでnは移動範囲の区間オプションに指定された値。

[メディアン移動範囲]オプションを使用するIR管理図では、シグマの値は次のように計算されます。

Equation shown here

ここで

MMRは、非欠測値の移動範囲のメディアン

d4(n)は、母標準偏差が1であるn個の独立した正規分布に従う確率変数の範囲の中央値、ここでnは移動範囲の区間オプションの値。

[R]オプションを使用するXBar管理図では、シグマの値は次のように計算されます。

Equation shown here

ここで

Ri = i番目のサブグループの範囲

ni = i番目のサブグループの標本サイズ

d2(ni)は、母標準偏差が1であるni個の独立した正規分布に従う確率変数の範囲の期待値

Nは、ni ≥ 2であるサブグループの数

[S]オプションを使用するXBar管理図では、シグマの値は次のように計算されます。

Equation shown here

ここで

ni = i番目のサブグループの標本サイズ

c4(ni)は、母標準偏差が1であるni個の独立した正規分布に従う確率変数の標準偏差の期待値

Nは、ni ≥ 2であるサブグループの数

siは、i番目のサブグループの標本の標準偏差

正規分布の工程能力指数

ここでは、正規分布に従うデータの工程能力指数の計算方法について具体的に説明します。

平均がμ、標準偏差がσの工程の場合、母集団に基づく工程能力指数は次のように定義されます。

Cp = Equation shown here

Cpl = Equation shown here

Cpu = Equation shown here

Cpk = Equation shown here

Cpm = Equation shown here

ここで

LSLは下限仕様限界

USLは上限仕様限界

Tは目標値

標本に基づく工程能力指数の場合、パラメータは推定値で置き換えられます。σの推定は、「工程能力分析」ウィンドウで指定した方法で行われます。変動の統計量を参照してください。

下側仕様限界と上側仕様限界のいずれかを指定しなかった場合、その仕様限界を含む工程能力指数は欠測値になります。

ヒント: 通常、工程能力指数は1.33以上であるべきと考えられています。工程能力指数が1.33である場合、正規分布に従っているならば、10万個の製品のうち不適合品となる個数の期待値は、約6個です。

工程能力指数の信頼区間

注: 工程能力指数の信頼区間は、「長期シグマ(全データから推定)」レポートにのみ表示されます。

Cpの100(1 - α)%の信頼区間は、次のように計算されます。

Equation shown here

ここで、

Equation shown here は、Cpの推定値

Equation shown here は、自由度n - 1のカイ2乗分布の(α/2)番目の分位点

nは、標本サイズ(観測値の個数)

Cpkの100(1 - α)%の信頼区間は、次のように計算されます。

Equation shown here

ここで

Equation shown here は、Cpkの推定値

Equation shown here は、標準正規分布の(1 - α/2)番目の分位点

nは、標本サイズ(観測値の個数)

CPMの100(1 - α)%の信頼区間は、次のように計算されます。

Equation shown here

ここで、

Equation shown here は、CPMの推定値

Equation shown here は、自由度γのカイ2乗分布の(α/2)番目の分位点

Equation shown here

nは、標本サイズ(観測値の個数)

Equation shown here は、標本平均

Tは、目標値

sは、全データから推定した長期シグマ

注: CPMの信頼区間は、目標値が下側仕様限界と上側仕様限界の中間にある場合のみ計算されます。

CPLとCPUの下側信頼限界と上側信頼限界は、Chou et al.(1990)の方法を使用して計算されます。

CPLの100(1 - α)%信頼限界(CPLLおよびCPLU)は、次の式を満たします。

Equation shown here ここでEquation shown here

Equation shown here ここでEquation shown here

ここで

tn-1(δ)は、自由度がn - 1、非心度パラメータがδの非心t分布に従う確率変数

Equation shown here は、Cplの推定値

CPUの100(1 - α)%信頼限界(CPULおよびCPUU)は、次の式を満たします。

Equation shown here ここでEquation shown here

Equation shown here ここでEquation shown here

ここで

tn-1(δ)は、自由度がn - 1、非心度パラメータがδの非心t分布に従う確率変数

Equation shown here は、Cpuの推定値

非正規分布の工程能力指数

この節では、非正規分布データの工程能力指数の計算方法について説明します。工程能力指数を一般化した式は、次のように定義されます。

Cp = Equation shown here

Cpk = Equation shown here

Cpm = Equation shown here

Cpl = Equation shown here

Cpu = Equation shown here

ここで

LSLは、下限仕様限界

USLは、上限仕様限界

Tは、目標値

Pαは、あてはめた分布の分位点

Cpmの計算に含まれるμσには、あてはめた分布の平均推定値と標準偏差推定値(分散推定値の平方根)が使われます。「パラメータ推定値」レポートのパラメータ間の関係と、あてはめた分布の期待値と分散の詳細については、『基本的な統計分析』の連続分布のあてはめ離散分布のあてはめを参照してください。

シグマクオリティ

シグマクオリティの値は、各比率(LSL未満、USL超え、限界外の合計数)に対し、次のように計算されます。

Equation shown here

ここで、

Pctは、「実測値%」列の値

Equation shown here は、標準正規分布の、累積確率(1 - Pct/100)における分位点

注: 「LSL未満」と「USL越え」の「実測値%」列の値を合計したものは、「限界外の合計数」の値と等しくなります。ただし、「LSL未満」と「USL越え」の「シグマクオリティ」列の値を合計しても、「限界外の合計数」のシグマクオリティにはなりません。これは、シグマクオリティが正規分布の分位点から求められ、加法的でないためです。

ベンチマークZ

ベンチマークZは、正規分析に基づく工程能力分析でのみ使用できます。ベンチマークZの値は、次のように計算されます。

Z Bench = Equation shown here

Z LSL = Equation shown here = 3 * Cpl

Z USL = Equation shown here = 3 * Cpu

ここで

LSLは、下限仕様限界

USLは、上限仕様限界

μは、標本平均

σは、標本標準偏差

Φ-11 - P(LSL) - P(USL)は、標準正規分布の(1 - P(LSL) -P(USL))番目の分位点

P(LSL) = Prob(X < LSL) = 1 - Φ(Z LSL)

P(USL) = Prob(X > USL) = 1 - Φ(Z USL)

Φは、標準正規分布の累積分布関数

より詳細な情報が必要な場合や、質問があるときは、JMPユーザーコミュニティで答えを見つけましょう (community.jmp.com).