管理図係数計算

　品質管理の上でも重要な管理図係数は、JIS(Z9020-2)に表形式で掲載されています。
　PCレベルによる標準偏差計算が一般的となっている現在では、この数表(シューハート管理図係数表：Shewhart Control Chart Constant)を用いる必然性はないとも言えますが、ここではあくまでも計算技術としての視点で解説していきます(PCを誰もが容易に使えない時代(1930年代)に、管理図係数を数表形式で提示したことには敬意を表したいと思います)。

　この管理図の意味や使い方などは別途解説予定ですが、ここでは管理図係数の中でも累積分布関数の積分を含んだ算出難易度の高い2つの係数を、Excel VBAを用いて計算する方法について概要を説明します。
　係数の定義式自体は公開されていますが、これをプログラムで数値計算するためには、本サイトで解説している数値積分法(シンプソン法)を駆使する必要があります。
　この係数計算を、RやMaximaといった言語で実施した例はありますが、VBAで実装した例があまりないため、参考になればと思います。

　なお、ここで紹介した手法を用いたExcel計算シートも掲載しておきます。
　(数式保護のためシートロックをかけていますがパスワードは設定していないので、興味のある方はロック解除して下さい)
　また、指定した管理図係数の値を返すExcel VBA用のコードも添付します。

　参考ですが、本文で説明している関数f1とf2のグラフィカルなイメージは次のようなものです(群の大きさ=5の場合)。
　どちらの関数も、変数xやyが±6程度でほぼゼロに収束していることに着目して下さい。