JISZ9003計量抜取検査(標準偏差既知)で上限合格判定値が既知の抜取方式
「計量抜取検査(標準偏差既知) (JISZ9003)がよくわからない」、「合格判定値と規格値の違いで求め方がどう変わるのかがわからない」など困っていませんか?
こういう疑問に答えます。
本記事のテーマ
JISZ9003計量抜取検査(標準偏差既知)で上限合格判定値が既知の抜取方式
- ①上限合格判定値についての関係式を導出
- ②上限合格判定値と上限規格値の導出方法の違い
- ③演習問題
- ④OC曲線を描く
下限規格値については、関連記事で確認ください。
JISZ9003計量抜取検査(標準偏差既知)で下限合格判定値が既知の抜取方式 JISZ9003計量抜取検査(標準偏差既知)で下限合格判定値が既知の抜取方式について解説します。 サンプル数n、下限合格判定値の関係式の導出やOC曲線の描き方を解説します。計量抜取検査をマスターしたい方は必見です。 |
①上限合格判定値についての関係式を導出
関係式を導出するためのモデル図を作成
次のような計量抜取検査を考えます。
● 上限合格判定値\(\bar{X_U}\)未満のロットは合格
● 上限合格判定値\(\bar{X_U}\)以上のロットは不合格
とする。前者はできるだけ合格させたいが、後者はできるだけ不合格にさせたいような抜取検査を考えたい。
モデル図を下図のように作ります。
このモデル図がしっかり作りこむことが意外と重要です。よく眺めてください。
できるだけ合格させたいp0はα=0.05(生産者危険)
できるだけ不合格にさせたいp1はβ=0.1(消費者危険)
の確率になるような抜取方式を検討します。
関係式を導出
モデル図から次の式が導かれます。見たらわかりますね。
計量抜取検査の理論は、モデル図から式を導出します
①\(σ_{\bar{x}}\)=\(\frac{σ}{\sqrt{n}}\)
②上限合格判定値\(\bar{X_U}\)の関係式を作ります。
ここから
・\(μ_0\)=\(\bar{X_U}\)-\(K_{α}σ_{\bar{x}}\)
・\(μ_1\)=\(\bar{X_U}\)+\(K_{β}σ_{\bar{x}}\)
サンプル数nの導出
式変形します。
\(μ_1\)-\(μ_0\)=\(K_{α}σ_{\bar{x}}\)+\(K_{β}σ_{\bar{x}}\)
よって、
n=\((\frac{K_α+K_β}{μ_1-μ_0})^2 σ^2\)
下限合格判定値が与えられた場合と同じ式になります。
上限合格判定値\(\bar{X_U}\)の導出
\(\bar{X_U}\)=\(μ_0\)+\(K_{α}σ_{\bar{x}}\)
=\(μ_0\)+\(K_{α} \frac{μ_1-μ_0}{ K_{α}+ K_{β}}\)
=\(\frac{ K_{β}μ_0+ K_{α}μ_1}{ K_{α}+ K_{β}}\)
なお、サンプル数nと、合格判定値の式は
下限合格判定値が与えられた場合と同じ式になります。
②上限合格判定値と上限規格値の導出方法の違い
上限規格値の場合
●合格判定係数k=k\((K_{α}, K_{β}, K_{p0}, K_{p1})\)
を導出
上限合格判定値の場合
●上限合格判定値\(\bar{X_U}\)=\((K_{α}, K_{β}, μ_0,μ_1 \))
を導出
サンプル数nを表現する変数が変わることと
合格判定係数kではなく、直接上限合格判定値を使うこと
の2点の違いがあります。
③演習問題
サンプル数nと上限合格判定値\(\bar{X_U}\)を使って、計量抜取検査のOC曲線が描けます。その前に演習問題を出して考えましょう。
あるスナック菓子のロットの平均値が150g以上ある場合は、できるだけ合格させたいが、152g以下の平均値をもつロットはできるだけ不合格にしたい。ただし、ロットの標準偏差は5gとわかっている。このとき、第1種の誤りであるα=0.05,第2種の誤りであるβ=0.10とした場合の抜取方式を決めよ。
うん、難しそう。。。でも1つずつ見ていきましょう。
まず、検査は抜取検査をやろうとしていますね。
次に、扱う変数は厚さという計量値を検査しようとしていますね。
最後に、上限合格判定値が決まっていますね。
サンプル数nと上限合格判定値\(\bar{X_U}\)を導出した公式から求めましょう。
まず、確率から\(K_{α}\)、\(K_{β}\)
がわかります。正規分布表を活用します。
\(K_{α}\)=1.645 (α=0.05のときのK値)
\(K_{β}\)=1.282(β=0.10のときのK値)
正規分布表に苦手意識があれば関連記事で復習しましょう。
【簡単】正規分布は怖くない!正規分布表や確率計算の求め方がすぐわかる 「正規分布とは何か?」、「正規分布の難解な式が理解できない」、「正規分布表の意味がわからない」など困っていませんか?本記事では、教科書やwebサイトより正規分布の基本やポイントをわかりやすく解説します。最も重要な正規分布を理解したい方は必見です。 |
サンプル数nは
n=\((\frac{K_α+K_β}{μ_0-μ_1})^2 σ^2\)
=\((\frac{1.645+1.282}{152-150})^2 5^2\)
=54
\(\bar{X_U}\)
=\(\frac{ K_{β}μ_0+ K_{α}μ_1}{ K_{α}+ K_{β}}\)
=\(\frac{ 1.282×150+ 1.645×152}{ 1.645+ 1.282}\)
=151.12
まとめると
(n, \(\bar{X_U}\))=(54,151.12)の値で、
平均値が151.12g以上ならロット合格、未満ならロット不合格
となります。
④OC曲線を描く
上の演習問題の結果をOC曲線で描きます。
OC曲線を描くための準備
なお、OC曲線を描くために、の関係式を再度書きます。
・\(μ_1\)=\(\bar{X_U}\)+\(K_{β} σ/\sqrt{n}\)
ここで、\(μ_1\),βを一般化して、
\(μ_1\)⇒μ
β⇒L(μ)
に変えます。慣れないとここの変化は無理矢理感がありますけど。
μ=\(\bar{X_U}\)+\(K_{ L(μ)} σ/\sqrt{n}\)
\(\frac{μ_\bar{X_U}}{σ} \sqrt{n}\)=\(K_{ L(μ)}\)
L(p)の作り方
- μを変数として値を振る。
- \(\frac{μ-\bar{X_U}}{σ} \sqrt{n}\)を計算する。
- μとL(μ)の関係からOC曲線を描く。
では、実際にやってみましょう。表にまとめます。
μ | μ-\(\bar{X_U}\) | \(\frac{μ-\bar{X_L}}{σ}\) | \(\frac{μ-\bar{X_L}}{σ} \sqrt{n}\) | L(μ) |
149 | -2.12 | -0.424 | -3.116 | 0.999 |
149.5 | -1.62 | -0.324 | -2.381 | 0.991 |
150 | -1.12 | -0.224 | -1.646 | 0.95 |
150.5 | -0.62 | -0.124 | -0.911 | 0.819 |
151 | -0.12 | -0.024 | -0.176 | 0.57 |
151.5 | 0.38 | 0.076 | 0.558 | 0.288 |
152 | 0.88 | 0.176 | 1.293 | 0.098 |
152.5 | 1.38 | 0.276 | 2.028 | 0.021 |
153 | 1.88 | 0.376 | 2.763 | 0.003 |
ここで、表の計算式をExcelの式を使って表現しています。
L(p)=1- (NORM.DIST(\(K_{L(p)}\)の値,0,1,TRUE))
OC曲線を描く
OC曲線です。計数抜取検査と似たような曲線になります。
計量抜取検査は式変形が多いですが、慣れましょう。
まとめ
JISZ9003計量抜取検査(標準偏差既知)で上限合格判定値が既知の抜取方式について、解説しました。
- ①上限合格判定値についての関係式を導出
- ②上限合格判定値と上限規格値の導出方法の違い
- ③演習問題
- ④OC曲線を描く
Warning: count(): Parameter must be an array or an object that implements Countable in /home/qcplanets/qcplanets.com/public_html/wp-content/themes/m_theme/sns.php on line 119