投稿者: QCプラネッツ

「群内変動と群間変動の分散の推定がわからない」、「QC検定®1級で頻出な群内変動と群間変動の問題が解けない」と困っていませんか？

こういう疑問に答えます。

本記事のテーマ

• ①R管理図ではσ=R/d2でいい理由
• ②群内変動と群間変動の分散を推定する公式
• ③群内変動と群間変動の分散を推定する演習問題1
• ④群内変動と群間変動の分散を推定する演習問題2

記事の信頼性

記事を書いている私は、QC検定®１級合格しましたが、管理図の係数表、群内変動・群間変動の解き方に疑問が残りました。そこで、管理図の理論を研究しました。その成果をブログで解説します。

①R管理図ではσ=R/d2でいい理由

R管理図に準拠した確率密度関数からσが推定できる

R管理図の確率密度関数や係数\(d_2\),\(d_3\)の導出については、関連記事で詳細に解説しています。

【必読】R管理図の変数d2,d3の導出が(半分)わかる R管理図の係数d2,d3はどうやって求めるか説明できますか？本記事では、範囲Rの確率密度関数を順序統計量の同時分布を使って導出し、途中までですが、d2,d3の導出方法を解説します。管理図をマスターしたい方は必見です。

Ｒ管理図のデータは確率密度関数に従うと仮定するため、
データによらず、σ=R/\(d_2\)と推定してよいのです。
単なる便利な式ではないことに注意しましょう。

なお、s管理図の場合は、E[s]=\(c_4\)σですね。
σ=s/\(c_4\)で推定してよいとなります。

R管理図の管理範囲内だけ使える

全体変動＝群内変動＋群間変動
で計算しますが、
R管理図で見ている群内変動だけです。

\(σ_T^2\)= \(σ_w^2\)+ \(σ_b^2\)
\(σ_w\)=R/\(d_2\)

群間変動、全変動はデータから別途から求めます。QC検定®１級なら、値が与えられています。

②群内変動と群間変動の分散を推定する公式

全体Xの場合

全体平均\(\bar{X}\)の場合

この区別が混合すると、わけがわからなくなります。
実際に演習問題で練習しましょう。見て⇒慣れて⇒ものにしましょう！

③群内変動と群間変動の分散を推定する演習問題1

演習問題1

解法

推定すべき変動の関係式をまず書きます。
\(σ^2\)=\(σ_W^2\)+\(σ_M^2\)+\(σ_L^2\)

(1)
●R管理図から、工作機内変動\(σ_W^2\)が公式で求められます。
\(\widehat{σ_w^2}\)=\((\frac{\bar{R}}{d_2 (n=4)})^2\)
=\((\frac{3.0+2.4}{2×2.059})^2\)=\(1.311^2\)=1.720

(2)
●R管理図から、工作機内変動と工作機間変動の和が公式で求められます。
\(\widehat{σ_w^2}\)+\(\widehat{σ_M^2}\)=\((\frac{R}{d_2 (n=8)})^2\)
(右辺)= \((\frac{R}{d_2})^2\)= \((\frac{4.3}{2.847})^2\)=\(1.510^2\)=2.281

工作機間変動\(σ_M^2\)は、
\(σ_M^2\)=2.281-1.720=0.5612

(3)
ロット間変動は、残りになります。
\(σ_L ^2\)=\(σ^2\)-\(σ_M^2\)-\(σ_L^2\)
\(σ_L ^2\)=\(2.2^2\)-1.720-0.5616=2.556

●まとめると、
(1)工作機内変動\(σ_W^2\)=1.720
(2)工作機間変動\(σ_M^2\)=0.5616
(3)ロット間変動\(σ_L^2\)=2.556

慣れるまで、何回も見てポイントを身に着けていきましょう。

④群内変動と群間変動の分散を推定する演習問題2

演習問題１と解き方は同じですが、もう１問見て、ポイントをおさえましょう。

演習問題2

変動の計算で、全体の場合と、全体平均の場合があります。この例題として演習問題を解いてみましょう。違いを意識して解きましょう。

全体Xの場合

全体平均\(\bar{X}\)の場合

●(1) 測定誤差分散の推定値\(\widehat{σ_M^2}\)は、
\(\widehat{σ_M^2}\)=\((\frac{R_{ij}}{2×10×d2(n=3)})^2\)
=\((\frac{15}{20×1.693})^2\)=0.1963

●(2-A) サンプル間変動の推定値\(\widehat{σ_S^2}\)は、
\(\widehat{σ_S^2}\)=\((\frac{R_i}{10×d2(n=2)})^2\)-\(\widehat{σ_M^2}\)
=\((\frac{8}{10×1.128})^2\)-0.1963=0.3067

●(2-B) サンプル間平均変動の推定値\(\widehat{σ_\bar{S}^2}\)は、
\(\widehat{σ_\bar{S}^2}\)=\((\frac{R_i}{10×d2(n=2)})^2\)-\(\frac{1}{3} \widehat{σ_M^2}\)
=\((\frac{8}{10×1.128})^2\)-0.1963/3=0.4376

●(3-A) ロット間変動の推定値\(\widehat{σ_L^2}\)は、
\(\widehat{σ_L^2}\)=\((\frac{R}{d2(n=10)})^2\)-\((\widehat{σ_S^2}+\widehat{σ_M^2})\)
=\((\frac{5}{3.078})^2\)-(0.3067+0.1963)=2.136

●(3-B) ロット間平均変動の推定値\(\widehat{σ_\bar{L}^2}\)は、
\(\widehat{σ_\bar{L}^2}\)=\((\frac{R}{d2(n=10)})^2\)-\(\frac{1}{2}(\widehat{σ_\bar{S}^2}+\frac{1}{3}\widehat{σ_M^2})\)
=\((\frac{5}{3.078})^2\)-\(\frac{1}{2}(0.4376+\frac{1}{3}\)0.1963)=2.387

●(4-A) 全変動\(σ_T^2\)は、
\(σ_T^2\)=\(σ_M^2\)+\(σ_S^2\)+\(σ_L^2\)
=0.1963+0.3067+2.136
=2.639

●(4-B) 全成分の平均変動\(σ_\bar{T}^2\)は、
\(σ_\bar{T}^2\)=\(σ_M^2\)+\(σ_\bar{S}^2\)+\(σ_\bar{L}^2\)
=0.1963+0.4376+2.387
=3.021

公式を活用して、それぞれの変動を導出できました。

QC検定®1級対策(管理図)

まとめ

計量値管理図の群内変動と群間変動の分散が推定できる方法を解説しました。

• ①R管理図ではσ=R/d2でいい理由
• ②群内変動と群間変動の分散を推定する公式
• ③群内変動と群間変動の分散を推定する演習問題1
• ④群内変動と群間変動の分散を推定する演習問題2

管理図

「R管理図の変数d2,d3の導出がわからない」と困っていませんか？

こういう疑問に答えます。

本記事のテーマ

d2,d3の式

●\(d_2\)=\(\displaystyle \int_{-\infty}^{\infty} [1-(1-φ(x))^n-(φ(x))^n]dx\)
●\(d_3\)=\(\sqrt{2\displaystyle \int_{-\infty}^{\infty} \displaystyle \int_{-\infty}^{y}f(x,y)dxdy-d_2^2 }\)
\(f(x,y)=1-φ(y)^n-(1-φ(x))^n+(φ(y)-φ(x))^n\)

読んでも理解ができない超難関な式です。でも、これをR管理図の係数表としてよく見かけます。でも、どうやってこの式になったのか？と気になるのは当然！

注意！

R管理図よりs管理図の方を使った方がよいかもです。
s管理図については理論式の導出ができるからです。

本記事のテーマ

• ①範囲Rの特性
• ②順序統計量の同時分布を確率密度関数とする
• ③d2の導出(わかる範囲で)
• ④d3の導出(わかる範囲で)
• ⑤係数\(d_2\),\(d_3\)の参考文献

記事の信頼性

記事を書いている私は、管理図の係数表、群内変動・群間変動の解き方に疑問が残りました。そこで、管理図の理論を研究しました。その成果をブログで解説します。

●Youtube動画にも解説しています。ご確認ください。

①範囲Rの特性

範囲Rはいい点も、困る点もあります。

>範囲Rのいい点

1. 「(最大)と(最小)の差」と計算しやすい
2. 理解しやすい、使いやすい

なので、計量値を扱う管理図のほとんどが\(\bar{X}\)-R管理図です。

>範囲Rの困る点

1. (最大)―(最小) ≥ 0と範囲Rの分布は0以上と限定
2. xが正のみな分布を表現する確率密度関数が超複雑になる

その結果、d2,d3が
●\(d_2\)=\(\displaystyle \int_{-\infty}^{\infty} [1-(1-φ(x))^n-(φ(x))^n]dx\)
●\(d_3\)=\(\sqrt{2\displaystyle \int_{-\infty}^{\infty} \displaystyle \int_{-\infty}^{y}f(x,y)dxdy-d_2^2 }\)
\(f(x,y)=1-φ(y)^n-(1-φ(x))^n+(φ(y)-φ(x))^n\)

とわけのわからない式になってしまいます。

②順序統計量の同時分布を確率密度関数とする

範囲Rに合う確率密度関数を探す

範囲Rは
●0以上の値
●大きい値と小さい値の2つの差
という特性があります。

これを表現できる確率密度関数が、
順序統計量の同時分布です。

1. 順序統計量って何？　難しい？⇒難しいです。
2. 同時分布って何？　2変数で表現する難解な式です。

急にレベルが高くなりました。順序統計量、順序統計量の同時分布については関連記事で解説します。まずは、「こういう関数があるんだ！」でかまいません。まずは使ってみることです。

順序統計量の同時分布

●「順序統計量が難しい」
●「同時分布はもっと難しい」
と、2段階で難しい話ですが、式だけ追いましょう。

順序統計量の同時分布を表現する確率密度関数は次のようになります。
\(f_{X(k) X(l)} (u,v)\)
=\(\frac{n!}{(k-1)!(l-k-1)!(n-l)!}\)\(F_{X(u)}^{k-1} f_{X(u)}\)\([F_{X(v)}-F_{X(u)}]^{l-k-1}\)\(f_{X(v)}\)\([1-F_{X(v)}]^{n-l}\)

範囲Rはx(n)-x(1)の差

順序統計量とは、個々の変数\(x_{i}\)について、
\(x_{1}\) ≤ \(x_{2}\) ≤ … ≤ \(x_{n}\)
の関係が成り立ちます。

範囲Rは最大と最小の差ですから、

R= \(x_{n}\) – \(x_{1}\) ≥ 0
が成り立ちます。
さらに\(x_{n}\) = \(x_{2}\), \(x_{1}\) = \(x_{1}\)にも注目しましょう。

範囲Rの確率密度関数

順序統計量の同時分布を表現する確率密度関数について、k=1,l=nを代入します。

\(f_{X(1) X(n)} (u,v)\)
=\(\frac{n!}{(1-1)!(n-1-1)!(n-n)!}\)\(F_{X(u)}^{1-1} f_{X(u)}\)\([F_{X(v)}-F_{X(u)}]^{n-1-1}\)\(f_{X(v)}\)\([1-F_{X(v)}]^{n-n}\)

=\(\frac{n!}{(n-2)!}\)\(f_{X(u)}\)\([F_{X(v)}-F_{X(u)}]^{n-2}\)\(f_{X(v)}\)
=\(n(n-1)\)\(f_{X(u)}\)\([F_{X(v)}-F_{X(u)}]^{n-2}\)\(f_{X(v)}\)

③d2の導出(わかる範囲で)

E[R]の立式

辻褄合わせですが、範囲Rを
R=\(x_{2}\)-\(x_{1}\) = (u-v)σ
に変えて積分します。(ちょっと無理があるけど)

期待値と分散の公式

よって、Rの期待値E[R]は次の式となります。
E[R]= \(\displaystyle \int R f(R) dR\)
=\(\displaystyle \int R \)\(n(n-1)\)\(f_{X(u)}\)\([F_{X(v)}-F_{X(u)}]^{n-2}\)\(f_{X(v)}dR\)
=σ\(\displaystyle \int \int (u-v) n(n-1)\)\(f_{X(u)}\)\([F_{X(v)}-F_{X(u)}]^{n-2}\)\(f_{X(v)}dudv\)

●次に順序統計量について関係式を使います。
\(\displaystyle \frac{d F_{X(u)}}{dx}\)=\( f_{X(u)}\)
\(\displaystyle d F_{X(u)}\)=\( f_{X(u)} dx\)
と変形し、これを使います。

Rの期待値E[R]は
E[R]=σ\(n(n-1)\displaystyle \int d F_{X(u)} \int (u-v)\)\([F_{X(v)}-F_{X(u)}]^{n-2}\)\( d F_{X(u)}\)

積分区間と、u⇒x1,v⇒x2と表記を変えます。
E[R]=σ\(n(n-1)\displaystyle \int_{-\infty}^{\infty} d F_1 \int_{-\infty}^{x_1} (x_1-x_2)\)\([F_1-F_2]^{n-2} dF_2\)

この式が、「新編統計数値表　河出書房　1952」P207と同じ式です。ここまでの導出は理解したのですが、ここからがまだわかっていません。

E[R]の導出

「新編統計数値表　河出書房　1952」P207によって、係数d2を導出します。

E[R]=σ\(n(n-1)\displaystyle \int_{-\infty}^{\infty} d F_1 \int_{-\infty}^{x_1} (x_1-x_2)\)\([F_1-F_2]^{n-2} dF_2\)
を\([F_1-F_2]^{n-2}\)を展開して、部分積分をしたうえで、まとめると次の式になるようです。

E[R]=σ\(n! \sum_{r=0}^{n-2} \frac{(-1)^r}{(r+1)!(n-r-1)!}\)\( \displaystyle \int_{-\infty}^{\infty}(1-F_1^{n-r-1})F_1^{r+1} dx_1\)
=σ\( \displaystyle \int_{-\infty}^{\infty}(1-F_1^n-(1-F_1)^n)dx_1\)
=\(d_2\)σ

とすると、係数\(d_2\)の式になるようです。実際\(F_1\)を正規分布の確率密度関数\(φ(x)\)に置き換えるとOKです。

ただし、完全に導出できたかは、今も研究中です。わかり次第、ブログを更新します。

④d3の導出(わかる範囲で)

E[\(R^2\)]の立式

期待値と分散の公式

E[R]=σ\(n(n-1)\displaystyle \int_{-\infty}^{\infty} d F_1 \int_{-\infty}^{x_1} (x_1-x_2)\)\([F_1-F_2]^{n-2} dF_2\)
でしたね。

●E[\(R^2\)]は
E[\(R^2\)]=\(σ^2 n(n-1)\displaystyle \int_{-\infty}^{\infty} d F_1 \int_{-\infty}^{x_1} (x_1-x_2)^2\)\([F_1-F_2]^{n-2} dF_2\)
となります。

この式も、ここからの導出は研究中ですが、「新編統計数値表　河出書房　1952」P207によると、次の結果になるそうです。

E[\(R^2\)]=\(σ^2\)2\(\displaystyle \int_{-\infty}^{\infty} \displaystyle \int_{-\infty}^{x_1} (1-F_1^n-(1-F_2)^n-(F_1-F_2)^n) dx_1 dx_2\)
Fをφに書き換えると
E[\(R^2\)]=\(σ^2\)2\(\displaystyle \int_{-\infty}^{\infty} \displaystyle \int_{-\infty}^{x_1} (1-φ(x_1)^n-(1-φ(x_2))^n-(φ(x_1)-φ(x_2))^n) dx_1 dx_2\)
と置きますね。

分散Vの立式

分散Vは、
V[R]= V[R]= E[\(R^2\)]-\(E[R]^2\)
=\(d_3 σ^2\)

平方根を取ると、
●D[R]=\(\sqrt{d_3}\)σ
●\(d_2\)=\(\displaystyle \int_{-\infty}^{\infty} [1-(1-φ(x))^n-(φ(x))^n]dx\)
●\(d_3\)=\(\sqrt{2\displaystyle \int_{-\infty}^{\infty} \displaystyle \int_{-\infty}^{y}f(x,y)dxdy-d_2^2 }\)
\(f(x,y)=1-φ(y)^n-(1-φ(x))^n+(φ(y)-φ(x))^n\)

と書けます。

⑤係数\(d_2\),\(d_3\)の参考文献

完全に導出できていませんので、参考文献を紹介します。
係数係数\(d_2\),\(d_3\)が分かったら教えてください。

引き続き研究して参ります。

参考文献

●「新編統計数値表　河出書房　1952」P207
最も詳細に書いていますが、計算の途中経過がいまいちよくわかりません。

●管理図法―品質管理教程 (1962年)
1986年改訂版もありますが、数式や理論は1962年の初版の方が詳しく解説している印象があります。

参考サイト

いくつか紹介しますが、導出過程まで解説したものはありません。

【水増し係数と割引き係数：不偏標準偏差と管理図係】数
c4,d2,d3についての解説があるが、導出は無い。

【エクセルＱＣ館　管理図の係数】
管理図の係数の式は紹介されているが、導出は無い。

【Deriving Control Chart Constants】
管理図の係数の式は紹介されているが、導出は無い。

【管理図係数計算】
管理図の係数の式は紹介されているが、導出は無い。

係数\(d_2\),\(d_3\)の導出を一番詳しく書いているのは、
「新編統計数値表　河出書房　1952」P207ですが、
導出の途中経過までは書いていないため、自分で調べる必要があるのが現状です。

また、シューハートの論文などを読みましたが、計算過程が分からず…でした。

わかった内容をすべてお伝えし、さらに導出過程の解明に努めていきます。

まとめ

R管理図の管理限界線の係数\(d_2\),\(d_3\)の導出を解説しました。

• ①範囲Rの特性
• ②順序統計量の同時分布を確率密度関数とする
• ③d2の導出(わかる範囲で)
• ④d3の導出(わかる範囲で)
• ⑤係数\(d_2\),\(d_3\)の参考文献

管理図

「s管理図の変数c4の導出がわからない」、「χ2乗分布の確率密度関数が変形できない」など困っていませんか？

こういう疑問に答えます。

本記事のテーマ

• ①s管理図の変数c4について
• ②不偏標準偏差vと母標準偏差σはχ2乗分布でつなぐ
• ③s管理図の変数c4の導出(χ2乗分布に慣れよう！)
• ④s管理図の管理限界の期待値E[v]と標準偏差D[v]がわかる

本記事でわかること

記事の信頼性

記事を書いている私は、管理図の係数表、群内変動・群間変動の解き方に疑問が残りました。そこで、管理図の理論を研究しました。その成果をブログで解説します。

困ったときは、最新のJIS規格を見ましょう。ちょっと高価ですが。

●Youtube動画でも解説しています。ご確認ください。

①s管理図の変数c4について

関連記事【重要】管理図(計量値)の変数の導出がわかるのように、管理図の管理限界は、
●E(a)±kD(a)
で表現され、aに変数、kに倍数を入れて管理します。

s管理図の変数aは不変標準偏差vです。不偏標準偏差vの確率密度関数から、期待値E(v),標準偏差D(v)を導出します。

その時に必要な変数が、c4です。

s管理図の変数c4

●\(c_4\)=\(\sqrt{\frac{2}{n-1}}\frac{Γ(\frac{n}{2})}{Γ(\frac{n-1}{2})}\) (v:不偏標準偏差の場合)
●\(c_4\)=\(\sqrt{\frac{2}{n}}\frac{Γ(\frac{n}{2})}{Γ(\frac{n-1}{2})}\) (s:標準偏差の場合)

訳が分からない式ですが、導出できます。

②不偏標準偏差vと母標準偏差σはχ2乗分布でつなぐ

関連記事【簡単】χ2乗分布がすぐ使いこなせる【初心者向け】のように、s,vとσの関係はχ2乗分布でつなぎます。χ2乗分布に慣れていく重要な式です。

\(χ^2\)=\(\frac{S(平方和)}{σ^2}\)
=\(\frac{v^2(不偏標準偏差)(n-1)}{σ^2}\)
から、vとσがつながります。

なお、標準偏差sの場合は、
\(χ^2\)=\(\frac{S(平方和)}{σ^2}\)
=\(\frac{s^2(標準偏差)(n)}{σ^2}\)
から、sとσがつながります。

χ2乗分布の確率分布関数を出しましょう。これが複雑ですが、そういう式と思ってください。簡単に説明すると、正規分布(\(e^{-x^2}\))の2乗和した関数がχ2乗分布ですね。

正規分布、χ2乗分布、t分布、F分布の関係も復習しましょう。関連記事で確認ください。

【簡単】χ2乗分布とt分布とF分布がすぐわかる【初心者向け】 「χ2乗分布とt分布とF分布がわからない」など、分布の特性を知らずに検定や推定、分散分析をしていませんか？本記事では、正規分布、χ2乗分布、t分布とF分布について、わかりやすく理解すべきポイントを解説します。χ2乗分布とt分布とF分布が理解したい方は必見です。

χ2乗分布の確率分布関数

χ2乗分布は2変数用意します。自由度nと変数xですね。

●\(χ^2(x)\)=\(\frac{1}{2^{n/2}Γ(\frac{n}{2})} x^{\frac{n}{2}-1}e^{-\frac{x}{2}}\)
Γ(z)=\(\displaystyle \int_{0}^{\infty} t^{z-1} e^{^t}dt\)

そういう関数だという理解でOKですし、変数c4の導出には、χ2乗分布の確率分布関数の一部を使って式変形するだけなので、慌てずに読み進めてください。

本記事は、①の「乗っかる」だけで行けます！

③s管理図の変数c4の導出(χ2乗分布に慣れよう！)

s管理図の管理限界の期待値E(s)は、
E(s)=c4σ
と書けます。このc4を導出してみましょう。

導出のポイント

1. 不偏標準偏差vとσをχ2乗分布でつなぐ

x=\(\frac{S(平方和)}{σ^2}\)
=\(\frac{v^2(不偏標準偏差)(n-1)}{σ^2}\)
とすると、
xは自由度(n-1)のχ2乗分布に従います。

よって、
v(不偏標準偏差)=\(\sqrt{\frac{x}{n-1}}σ\)
と書けます。
(ちなみに、s(標準偏差)= \(\sqrt{\frac{x}{n}}σ\)
とも書けます。)

また、準備として、自由度(n-1)のχ2乗分布の確率分布関数の式を用意します。
\(χ_{n-1}^2(x)\)=\(\frac{1}{2^{\frac{n-1}{2}}Γ(\frac{n-1}{2})} x^{\frac{n-1}{2}-1}e^{-\frac{x}{2}}\)

期待値E[v]の立式

期待値E[v]は定義では、
●E[v]= \(\displaystyle \int_{0}^{\infty} vf(v) dv \)
＊f(v)にχ2乗分布を代入
＊積分区間:χ2乗分布は2乗だけに負はないため、[0,∞]で積分です。

期待値と分散の公式

よって、E[v]は
●E[v]= \(\displaystyle \int_{0}^{\infty} vf(v) dv\)
=E[\(\sqrt{\frac{x}{n-1}}σ\)]
= \(\displaystyle \int_{0}^{\infty} \sqrt{\frac{x}{n-1}}σ \frac{1}{2^{\frac{n-1}{2}}Γ(\frac{n-1}{2})} x^{\frac{n-1}{2}-1}e^{-\frac{x}{2}} dx\)

ちなみに、標準偏差sについて立式すると　n-1⇒nに代わります。

●E[s]= \(\displaystyle \int_{0}^{\infty} sf(s) ds\)
=E[\(\sqrt{\frac{x}{n}}σ\)]
= \(\displaystyle \int_{0}^{\infty} \sqrt{\frac{x}{n}}σ \frac{1}{2^{\frac{n-1}{2}}Γ(\frac{n-1}{2})} x^{\frac{n-1}{2}-1}e^{-\frac{x}{2}} dx\)
標準偏差sと不偏標準偏差vで区別します。ややこしいですが。

期待値E[v]の導出

●E[v] =\(\displaystyle \int_{0}^{\infty} \sqrt{\frac{x}{n-1}}σ \frac{1}{2^{\frac{n-1}{2}}Γ(\frac{n-1}{2})} x^{\frac{n-1}{2}-1}e^{-\frac{x}{2}} dx\) (式1)
の定数項を積分の外に出しましょう。

(式1)= \(\sqrt{\frac{1}{n-1}}σ\displaystyle \int_{0}^{\infty} \sqrt{x} \frac{1}{2^{\frac{n-1}{2}}Γ(\frac{n-1}{2})} x^{\frac{n-1}{2}-1}e^{-\frac{x}{2}} dx\) (式2)

次にxをまとめます。

(式2)= \(\sqrt{\frac{1}{n-1}}σ\displaystyle \int_{0}^{\infty} \frac{1}{2^{\frac{n-1}{2}}Γ(\frac{n-1}{2})} x^{\frac{n}{2}-1}e^{-\frac{x}{2}} dx\) (式3)

ここでよく積分の中を見ると、複雑な式ですが、元のχ2乗分布の式に似ていることがわかります。つまり、すべて計算せずに、χ2乗分布の関数をそのまま使えばよいのです。

(式3)= \(\sqrt{\frac{1}{n-1}}σ \frac{2^{\frac{n}{2}}Γ(\frac{n}{2})}{2^{\frac{n}{2}}Γ(\frac{n}{2})}\displaystyle \int_{0}^{\infty} \frac{1}{2^{\frac{n-1}{2}}Γ(\frac{n-1}{2})} x^{\frac{n}{2}-1}e^{-\frac{x}{2}} dx\) (式4)

(式4)のように、分母分子に同じ \(2^{\frac{n}{2}}Γ(\frac{n}{2})\)を入れます。
分母にある、\(2^{\frac{n}{2}}Γ(\frac{n}{2})\)と\(2^{\frac{n-1}{2}}Γ(\frac{n-1}{2})\)を入れ替えます。

(式4)= \(\sqrt{\frac{1}{n-1}}σ \frac{2^{\frac{n}{2}}Γ(\frac{n}{2})}{2^{\frac{n-1}{2}}Γ(\frac{n-1}{2})}\displaystyle \int_{0}^{\infty} \frac{1}{2^{\frac{n}{2}}Γ(\frac{n}{2})} x^{\frac{n}{2}-1}e^{-\frac{x}{2}} dx\) (式5)

(式5)の
●\(\frac{2^{\frac{n}{2}}}{2^{\frac{n-1}{2}}}\)=\(\sqrt{2}\)
●\(\displaystyle \int_{0}^{\infty} \frac{1}{2^{\frac{n}{2}}Γ(\frac{n}{2})} x^{\frac{n}{2}-1}e^{-\frac{x}{2}} dx\)
=\(\displaystyle \int_{0}^{\infty} χ_n^2(x) dx\)
=1
となります。よって、

(式5)= \(\sqrt{\frac{2}{n-1}} \frac{Γ(\frac{n}{2})}{Γ(\frac{n-1}{2})}\)σ (式6)

まとめると、

ちなみに、標準偏差sについて解くと、

χ2乗分布の式が複雑ですが、ゆっくり導出を読めば、難しくないことがわかりますね。数式の見た目でビビらないことも重要です。

④s管理図の管理限界の期待値E[v]と標準偏差D[v]がわかる

導出のポイント

より期待値E[v]がわかりました。D[v]も導出しましょう。

期待値と分散の公式をもう一度確認しましょう。複雑な式でも公式通り代入しているだけです。

期待値と分散の公式

D[v]の式を作ります。

期待値D[s]の導出

D[v]は、
D[s]= E[\(v^2\)]-\(E[v]^2\)
です。

期待値E[\(v^2\)]の導出

さて、E[\(v^2\)]はどうしましょうか？

vはそもそも不偏標準偏差の期待値は母分散σとみてよいでしょう。
よって,
E[\(v^2\)]=\(σ^2\)

ちなみに、標準偏差sの場合のE[\(s^2\)]は、
もともと
\(χ^2\)=\(\frac{S(平方和}{σ^2}\)
から
\(χ^2\)=\(\frac{s(標準偏差)^2 n}{σ^2}\)=\(\frac{v(不偏標準偏差)^2 (n-1)}{σ^2}\)
の関係が成り立つので、
\(s^2\)=\(\frac{n-1}{n} v^2\)
となります。

なお、E[\(v^2\)]=\(σ^2\)より、
E[\(s^2\)]=\(\frac{n-1}{n} \) E[\(v^2\)]=\(\frac{n-1}{n} σ^2\)
となります。

まとめると

不偏標準偏差vの場合は、
D[v]= E[\(v^2\)]-\(E[v]^2\)
=\(σ^2-c_4^2 σ^2\)
=\((1-c_4 ^2 )σ^2\)
(\(c_4\)=\(\sqrt{\frac{2}{n-1}} \frac{Γ(\frac{n}{2})}{Γ(\frac{n-1}{2})}\))

標準偏差sの場合は、
D[v]= E[\(v^2\)]-\(E[v]^2\)
=\(\frac{n-1}{n}σ^2-c_4^2 σ^2\)
=\((\frac{n-1}{n}-c_4 ^2 )σ^2\)
(\(c_4\)=\(\sqrt{\frac{2}{n}} \frac{Γ(\frac{n}{2})}{Γ(\frac{n-1}{2})}\)

c4,期待値E、標準偏差Dのまとめ

s管理図における管理限界の導出を解説しました。

標準偏差sと不偏標準偏差vによって、式が若干異なります。
なお、JISZ9020では不偏標準偏差vの場合の値が管理限界の係数表に載っています。

教科書やwebサイトによっては、標準偏差sと不偏標準偏差vの両方が書いていますので、両方解説しました。

c4

E

D

まとめ

s管理図の管理限界線の係数の導出を解説しました。

• ①s管理図の変数c4について
• ②不偏標準偏差vと母標準偏差σはχ2乗分布でつなぐ
• ③s管理図の変数c4の導出(χ2乗分布に慣れよう！)
• ④s管理図の管理限界の期待値E[v]と標準偏差D[v]がわかる

「シューハート管理図の係数表にあるA2,D3,d2など変数って何？」、「係数表の変数はどうやって求めるのかわからない」など困っていませんか？

こういう疑問に答えます。

本記事のテーマ

• ①変数の基本は期待値と標準偏差から導出
• ②\(\bar{X}\)管理図の場合
• ③s管理図の場合
• ④R管理図の場合

記事の信頼性

記事を書いている私は、管理図の係数表、群内変動・群間変動の解き方に疑問が残りました。そこで、管理図の理論を研究しました。その成果をブログで解説します。

困ったときは、最新のJIS規格を見ましょう。ちょっと高価ですが。

●Youtube動画でも解説しています。ご覧ください。

①変数の基本は期待値と標準偏差から導出

変数の一覧

JISZ9020-2の表2「管理限界線を計算するための係数」から変数一覧を出します。

変数Aは簡単に導出できそうですが、
変数d2,d3は意味不明な式ですね。

もし、式や値が間違っていたらどうしますか？ JISに文句言っても、あなたの顧客は満足しませんよね。自分が使うものは、対象の良し悪しや意味を理解して使うべきです。

変数基本形

変数の種類が多いですが、すべて同じ形で導出します。

簡単な式ですね。それともう１つ。

変位Xの確率密度関数は正規分布なので簡単ですが、
標準偏差sと範囲Rについては別途関連記事で確率密度関数からの導出を解説します。

では、個別に解説していきます。

②\(\bar{X}\)管理図の場合

A,\(A_2\),\(A_3\)の3つについてです。

変位Xは正規分布で考える

正規分布を描きます。これから解説する、標準偏差s、範囲Rの確率密度関数と比較してください。

正規分布における、期待値Eと標準偏差Dは、
●E(X)=\(\bar{X}\)
●D(X)=σ
ですね。これは正規分布を勉強すればわかりますので導出は割愛します。

変数A,\(A_2\),\(A_3\)の導出

変数A の導出

E(X)±kD(X)
=\(\bar{X}\)±\(\frac{k}{\sqrt{n}}σ\)
= \(\bar{X}\)±Aσ

●A=\(\frac{k}{\sqrt{n}}\)

変数,\(A_2\)の導出

E(X)±kD(X)
= \(\bar{X}\)±\(\frac{k}{\sqrt{n}}σ\)
= \(\bar{X}\)±\(\frac{k}{\sqrt{n}} \frac{R}{d_2} \)
= \(\bar{X}\)±\(\frac{k}{d_2 \sqrt{n}} R\)
= \(\bar{X}\)±\(A_2\)σ

●\(A_2\)=\(\frac{k}{d_2 \sqrt{n}} R\)

変数,\(A_3\)の導出

E(X)±kD(X)
= \(\bar{X}\)±\(\frac{k}{\sqrt{n}}σ\)
=\(\bar{X}\)±\(\frac{k}{\sqrt{n}}\frac{s}{c_4}\)
=\(\bar{X}\)±\(\frac{k}{c_4 \sqrt{n}}s\)
= \(\bar{X}\)±\(A_3 s\)

●\(A_3\)=\(\frac{k}{c_4 \sqrt{n}} \)

③s管理図の場合

変数(B_3),(B_4),(B_5),(B_6)を導出します。

標準偏差sはχ2乗分布で考える

χ2乗分布を描きます。これから解説する、標準偏差s、範囲Rの確率密度関数と比較してください。

χ2乗分布における、期待値Eと標準偏差Dの導出は関連記事で詳細に解説します。
本記事は結果だけ使います。
●E(s)=\(c_4 σ\)
●D(s)=\(\sqrt{1-c_4^2}σ\)
になります。

関連記事にあるように、s管理図の管理限界を導出する際、標準偏差sと不偏標準偏差vを使った場合では若干式が異なります。本記事では、JISZ9020の管理限界係数表に準拠して、不偏標準偏差vを使った場合を解説します。

具体的には、

c4

E

D

あまり、この区別が教科書や他のwebサイトでは書いていないため、違いがわかるように解説します。

【必読】s管理図の変数c4と管理限界の導出がわかる s管理図の管理限界を求めるc4と管理限界値の導出を解説します。χ2乗分布、平方和、標準偏差の関係式を使って、意外と簡単に係数c4が導出できます。さらに、標準偏差と不偏標準偏差によって、若干式が異なる点も詳しく解説します。管理図をマスターしたい方は必見です。

変数\(B_3\),\(B_4\),\(B_5\),\(B_6\)の導出

変数\(B_3\),\(B_4\)の導出

●σが既知の場合、変数\(B_3\),\(B_4\)を使います。

E(s)±kD(s)
=\(c_4 σ\)±k \(\sqrt{1-c_4^2}σ\)
=(1±\(\frac{k}{c_4}\sqrt{1-c_4^2}\))\(c_4\)σ
(ここで、s=\(c_4\)σ)

よって
●\(B_3\)=\(1-\frac{k}{c_4}\sqrt{1-c_4^2}\)
●\(B_4\)=\(1+\frac{k}{c_4}\sqrt{1-c_4^2}\)

なお、\(B_3\)は0以上としたいので、まとめます。
●\(B_3\)=\(max(0,1-\frac{k}{c_4}\sqrt{1-c_4^2})\)
●\(B_4\)=\(1+\frac{k}{c_4}\sqrt{1-c_4^2}\)

変数\(B_5\),\(B_6\)の導出

●σが未知の場合、変数\(B_5\),\(B_6\)を使いますが、変数\(B_3\),\(B_4\)の導出方法から変わりません。

E(s)±kD(s)
=\(c_4 σ\)±k\(\sqrt{1-c_4^2}σ\)
=\((c_4 \)±k\(\sqrt{1-c_4^2})σ\)
(ここで、σは未知なので、sとσの関係式はありません。)

よって
●\(B_5\)=\(c_4 \)-k\(\sqrt{1-c_4^2}\)
●\(B_6\)=\(c_4 \)+k\(\sqrt{1-c_4^2}\)

なお、\(B_5\)は0以上としたいので、まとめます。
●\(B_5\)=\(max(0, c_4 \)-k\(\sqrt{1-c_4^2})\)
●\(B_6\)=\(c_4 \)+k\(\sqrt{1-c_4^2}\)

④R管理図の場合

変数\(D_1\),\(D_2\),\(D_3\),\(D_4\)を導出します。

範囲Rは順序統計量の同時分布で考える

範囲Rは計算しやすいですが、確率密度関数の導出が劇難です。

1. 範囲R は0以上であること
2. 変数Xと同じ1次式(標準偏差sは2次式)

上の2つを満たす確率密度関数の導出が難しくなります。詳細は関連記事をご覧ください。日本で、範囲Rの確率密度関数から\(d_2\),\(d_3\)を端折らず導出できる人はいないかもしれないくらいです。

範囲Rの確率密度関数イメージを描きます。

順序統計量の同時分布における、期待値Eと標準偏差Dは、
●E(R)=\(d_2\)σ
●D(R)=\(d_3\)σ
とします。

ただし、
●\(d_2\)=\(\displaystyle \int_{-\infty}^{\infty} [1-(1-φ(x))^n-(φ(x))^n]dx\)
●\(d_3\)= \(\sqrt{2\displaystyle \int_{-\infty}^{\infty} \displaystyle \int_{-\infty}^{y}f(x,y)dxdy-d_2^2 }\)
( \(f(x,y)=1-φ(y)^n-(1-φ(x))^n+(φ(y)-φ(x))^n\))
と超複雑な式です。

変数\(D_1\),\(D_2\)の導出

変数\(D_1\),\(D_42\)の導出

●σが既知の場合、変数\(D_1\),\(D_2\)を使います。

E(R)±kD(R)
=\(d_2\)σ±k\(d_3\)σ

よって
●\(D_1\)=\(d_2 σ\)-\(k d_3\)σ
●\(D_2\)=\(d_2 σ\)+\(k d_3\)σ

なお、\(D_1\)は0以上としたいので、まとめます。
●\(D_1\)=\(max(0,d_2-k d_3)\)
●\(D_2\)=\(d_2 σ\)+\(k d_3\)

変数\(D_3\),\(D_4\)の導出

●σが未知の場合、変数\(D_3\),\(D_4\)を使いますが、変数\(D_1\),\(D_2\)の導出方法から変わりません。

E(R)±kD(R)
=\(d_2\)σ±k\(d_3\)σ
=\(d_2\)\((1±\frac{k d_3}{d_2})\)σ

よって
●\(D_3\)=\((1-\frac{k d_3}{d_2})\)
●\(D_4\)=\((1+\frac{k d_3}{d_2})\)

なお、\(D_3\)は0以上としたいので、まとめます。
●\(D_3\)=\(max(0, 1-\frac{k d_3}{d_2})\)
●\(D_4\)=\((1+\frac{k d_3}{d_2})\)

シューハートの管理図の管理限界線の係数の導出を解説しました。基本はE(a)±kD(a)ですべて導出できることがわかりました。

まとめ

シューハートの管理図の管理限界線の係数の導出を解説しました。

• ①変数の基本は期待値と標準偏差から導出
• ②\(\bar{X}\)管理図の場合
• ③s管理図の場合
• ④R管理図の場合

管理図

「シューハート管理図の管理線公式をまとめた表が無い？」、「係数表をまとめた表がない」など困っていませんか？

こういう疑問に答えます。

本記事のテーマ

• ①シューハート管理図の管理線公式
• ②シューハート管理図の係数表

記事の信頼性

記事を書いている私は、QC検定®１級合格しましたが、管理図の係数表、群内変動・群間変動の解き方に疑問が残りました。そこで、管理図の理論を研究しました。その成果をブログで解説します。

困ったときは、最新のJIS規格を見ましょう。ちょっと高価ですが。

●Youtube動画でも解説しています。ご覧ください。

①シューハート管理図の管理線公式

管理図の種類、仮定する確率分布、中心、管理限界を一覧にまとめます。

苦手な管理図、見たことが無い管理図が無いように注意しましょう。

②シューハート管理図の係数表

まずは、表から値を使いこなす練習をしましょう。
慣れたら値の導出方法を調べたり、考えたりしてください。

JISZ9020-2(2016)から抜粋

管理図を使った解き方は教科書や他のブログ、webサイトに解説があるので、割愛します。

本記事では、試験前に１枚で全部わかる一覧表を解説しました。

まとめ

シューハートの管理図の公式、係数をまとめた表を解説しました。

• ①シューハート管理図の管理線公式
• ②シューハート管理図の係数表

管理図

「管理図はシューハート管理図だけやればいいのか？」、「解き方はわかるが、理論がいまいちピンとこない」など困っていませんか？

こういう疑問に答えます。

本記事のテーマ

• ①管理図を勉強した人のありがちな考え方
• ②管理図の最新JIS規格は3部しかない
• ③一般指針とシューハートの管理図の違い
• ④管理図は自分で考えて設計すべき

記事の信頼性

記事を書いている私は、QC検定®１級合格しましたが、管理図の係数表、群内変動・群間変動の解き方に疑問が残りました。そこで、管理図の理論を研究しました。その成果をブログで解説します。

困ったときは、最新のJIS規格を見ましょう。ちょっと高価ですが。

●Youtube動画でも解説しています。

①管理図を勉強した人のありがちな考え方

管理図でよく勉強する内容（試験で頻出）は、
(i) シューハートの管理図の種類を覚え、LCL,UCLの導出式を暗記
(ii) 計数値管理図の管理限界は、管理図係数表から導出
(iii) 郡内変動、群間変動の分散の式を暗記し、σの推定値をR/d2で導出

この3つがあれば、QC検定®は全級合格できます。

しかし、次の4つの疑問が残りました。

1. データの特性によらず、管理図係数表の値で管理限界が決まるのは正しいの？
2. 管理図係数表の値はだれがどうやって数学的に正しく導出したのか？
3. 全データを管理しているはずなのに、σの推定値を導出する意味はあるのか？
4. 郡内変動、群間変動の分散の式は正しいのか？

QC検定®１級は一発合格しましたが、管理図の問題は壊滅状態でしたし、今も理解できない解法があります。

QC検定®のための管理図ではなく、自分で考えて工程管理できる管理図を設計してほしいです。

②管理図の最新JIS規格は3部しかない

管理図の最新JIS規格

1. 管理図-一般指針　JISZ9020-1(2016)
2. 管理図-シューハート管理図 JISZ9020-2(2016)
3. 管理図-累積和管理図 JISZ9020-48(2018)

ここで、「累積和管理図 JISZ9020-48(2018)」は使う場面がほぼないので割愛します。

過去には、シューハート管理図 JIS9021 1998がありましたが、廃止され、シューハート管理図 JISZ9020-2(2016)に統合されました。最新版をみると、過去の規格よりシューハート管理図の内容が薄まったように感じます。

それでも、「シューハートの管理図」の存在感は強いです。しかし、
その前に「一般指針」を理解し、
「シューハートの管理図」は「一般指針」の一部でしかない
と理解しましょう。

③一般指針とシューハートの管理図の違い

自分で管理図を考えた場合と、シューハートの管理図を適用した場合では、同じ考え方になるところと、考えが異なるところがあります。

シューハートの管理図で「？」な部分は、自分で理論を作って、対象となる工程管理の妥当性を担保してほしいからです。

シューハート管理図の疑問点

1. データの特性によらず、管理図係数表の値で管理限界が決まるのは正しいの？
2. 管理図係数表の値はだれがどうやって数学的に正しく導出したのか？
3. 全データを管理しているはずなのに、σの推定値を導出する意味はあるのか？
4. 郡内変動、群間変動の分散の式は正しいのか？

一般指針とシューハートの管理図が同じである場合

管理図で調べたいこと、工程の合否など基本的な考え方は同じです。
だから、最初からシューハート管理図に入る傾向が強いと考えます。

一般指針とシューハートの管理図が異なる場合

管理図を自分で考えて作ると、シューハート管理図との違いがはっきりわかるようになります。

④管理図は自分で考えて設計すべき

管理図は手段です。調べたい目的は工程管理の異常を調べることですね。

管理図を使うときに考えること

対象となる工程管理の合否を調べるために必要な事項を挙げます。

• ①検査対象
• ②サンプル数の規模
• ③調べたい特性の値(計数値、計量値)
• ④工程の合否判定基準(何σなのか？)
• ⑤工程が異常だったら何をすべきか？

以上の項目をまず考えるはずです。
シューハートの管理図はまだ、出てきませんね。

シューハートの管理図は、見せ方の１つの手段として活用するべきです。

最初から、「pn管理図使おう」などしないことです。

異常時の対策と改善が実務では最も重要になります。改善提案を考えるときは、自分で考えて作った管理図の方が考えやすいし、協議しやすいはずです。

でも、今はExcelあれば十分ですね。
むしろ、理論や理由、工程異常のリスクと改善化が重要ではないでしょうか？

計算する手法より、頭で考えることがＱＣでは求められます。

まとめ

以上、管理図の作る際に、自分で考えて管理図を作ることの重要性を解説しました。

• ①管理図を勉強した人のありがちな考え方
• ②管理図の最新JIS規格は3部しかない
• ③一般指針とシューハートの管理図の違い
• ④管理図は自分で考えて設計すべき

管理図