# 実験計画法についてわかりやすく解説！Pythonで実際に実装！ > 当サイト【スタビジ】の本記事では、品質管理の分野においてよく用いられる実験計画法について解説していきます！考え方や事例を説明した後、Pythonで実装していきます！ - URL: https://toukei-lab.com/design-of-experiments - 公開: 2023.11.01 / 更新: 2025.08.18 - サイト: スタビジ（https://toukei-lab.com） :::talk{who="umatan"} 当サイト【スタビジ】の本記事では、品質管理の分野においてよく用いられる実験計画法について解説していきます！考え方や事例を説明した後、Pythonで実装していきます！ ::: こんにちは！スタビジ編集部です！ 統計学の応用分野である**実験計画法(Design Of Experiments)**。 近年発展している機械学習手法などと比較するとかなり単純な理論構成ですが、今もなお品質管理などの分野において頻繁に用いられている手法です。 この記事では、そんな実験計画法の特徴とPythonでの実装を見ていきたいと思います。 :::talk{who="umatan"} 概要や事例を説明した後、Pythonで実装していくよ！ ::: 以下の動画でも詳しく解説しているのであわせてチェックしてみてください！ ::youtube{#8Cx6jaQIS1Y} ## 実験計画法とは ![Content](/wp-content/uploads/2017/11/Content-amico.png) 実験計画法とは、**効率のよい実験計画を設計し、結果を適切に解析する**ことを目的とした方法論です。 フィッシャーが20世紀初頭に創始しました。 :::talk{who="umatan"} [最尤法](/maximum_likelihood_method)を編み出したとした人としても知られるあのフィッシャーだね！ ::: 実験計画法には様々な手法がありますが、それらを使っていく上で共通に理解しておくべき**2つのポイント**があります。 :::point **1. 平方和の分解** **2. 交互作用** ::: :::talk{who="robotan"} あれ？両方とも分散分析で聞いたような、、、 ::: :::talk{who="umatan"} そう！実は分散分析は実験計画法で用いられる分析でもあるんだ！ ::: 少しわかりづらいので図で表すとこんな感じ。 ![実験計画法　分散分析　違い](/wp-content/uploads/2023/10/DOE_ANOVA_diff.jpg) つまり、実験計画法という方法論を用いて**効率のよい実験を計画し**、分散分析を用いて**統計的に解析を行う**ということです（いきなり調査データに対して分散分析を行うこともあります）。 ただし実験計画から解析までをセットで行うことが多いので、入門書などでは**分散分析を行うプロセスまで含めて実験計画法**として紹介されていることがほとんどです。 ゆえに、**平方和の分解と交互作用はしっかりと押さえておきましょう。** :::talk{who="umatan"} 分散分析、平方和の分解、交互作用についてはこれらをチェック！ ::: ## 実験計画の3原則 ![stories Data-pana](/wp-content/uploads/2021/01/Data-pana-e1615808087658-1024x875.png) 実験計画法では統計的判断を活用しています。 そしてその基準は**データから推定される誤差**です。 ゆえに、誤差を精度よく求めることが誤った判断をしないために重要であると言えます。 そこで実験計画法の創始者であるフィッシャーは、実験の場を以下の3つの原則に従って管理することを提唱しました。 :::box **1. 反復(replication)**\ **2. 無作為化(randomization)**\ **3. 局所管理(local control)** ::: これらは**実験計画の3原則**と呼ばれています。 ### 1. 反復(replication) 誤差を評価するためには、**同一の条件下で実験を繰り返すこと**が必要です。 また、繰り返しのデータの平均をとることにより、母平均やその差の推定誤差も向上していきます。 ### 2. 無作為化(randomization) 実験で発生する誤差は偶然誤差だけではありません。 例えば、室温を特に制御せずに実験を行うとき、朝から昼にかけて室温は徐々に上昇し、夜に向けては徐々に下降していきます。 このような誤差を**系統誤差**と呼びます。 系統誤差の管理や評価が難しい場合、実験条件を時間的、空間的に無作為（ランダム）にわりつけることによって、**系統誤差を偶然誤差に転化することが可能**になります。 ### 3. 局所管理(local control) 局所管理の原則とは、実験の場を適当なブロックに分け、ブロック内では条件をなるべく均一になるようにし、効果の検出力と比較の精度の向上を図ることです。 これは、**系統誤差を回避するためのもう1つの方法**です。 系統誤差を発生させるおそれのある時間的、空間的要因をブロックに分け、そのなかに比較したい全ての処理を入れることで、系統誤差の大きい部分はブロック間の差異として取り出し、その分誤差のばらつきを減らすことができます。 :::talk{who="umatan"} これらの原則の下で実験計画法が成り立っているということもしっかりと押さえておこう！ ::: ## 実験計画法の種類 ![python](/wp-content/uploads/2019/08/Admin-amico-min-e1613730531228-1024x914.png) 実験計画法には、以下のような手法があります。 :::box **1. 一元配置法\ 2\. 二元配置法\ 3\. 多元配置法\ 4\. 直交表を用いた実験計画\ 5\. 乱塊法\ 6\. 分割法\ 7\. 枝分かれ実験** ::: 様々な実験の条件や目的に合った手法を選択し、解析まで行なっていきます。 :::talk{who="robotan"} 一元配置法とかって分散分析の名前じゃないの？？ ::: :::talk{who="umatan"} 一つの因子を取り上げる（配置する）って意味では同じだね！ただ、それぞれを「一元配置分散分析」と「一元配置実験」って言い方をする場合もあるよ！ ::: :::talk{who="robotan"} なるほど！ ::: ## 適用例(直交表を用いた実験計画) ![Business Plan](/wp-content/uploads/2020/07/Business-Plan-amico-e1613688762889-1024x903.png) ここから、**直交表を用いた実験計画**の適用例を説明していきます。 ### 1. 語句の説明 [分散分析の記事](/analysis-of-variance)でも語句の説明はされていますが、改めてここで説明しておきます。 **因子**は、ばらつきの要因のことを指します。 **水準**は、因子を量的、あるいは質的に変えた具体的な条件です。 また、実験の結果得られるデータのことを**特性値**と呼びます。 以下の図は、それぞれの語句の例を示したものです。 ![実験計画法　語句の説明](/wp-content/uploads/2023/10/DOE_word.jpg) ### 2. 直交表とは 一元配置法や二元配置法、多元配置法は、まとめて**要因配置実験**と呼ばれます。 要因配置実験は、簡単に言うと**全ての水準組み合わせで実験を行う方法**です。 しかし、たくさんの因子を取り上げた要因配置実験では、実験回数がかなり多くなってしまうというデメリットがあります。 例えば7つの因子（各2水準）を取り上げた場合、全ての因子間に交互作用が存在すると仮定すると、27=128回の実験が必要になり、かなり非効率です。 そこで、実験回数を減らしてくれる**直交表の出番**というわけです。 ![実験計画法　L8直交表](/wp-content/uploads/2023/10/DOE_L8array.jpg) 上の表は、2水準系直交表の1つである**L8直交表**です（3水準系直交表などもあります）。 L8直交表には、2水準の因子を最大7つまで（=列の数）適用することができます。 表内の数字は各因子の水準を表しており、**任意の2因子間の水準の組み合わせが均等になるように配置されています。** つまり、L8直交表を用いる場合、たった8回（=行の数）の実験で最適な水準組み合わせを見つけることができるのです。 :::talk{who="robotan"} そんなに少なくていいのか！ ::: ただし、直交表を用いる際、以下のことに注意してください。 :::warn **2因子間の交互作用を「検討すべきもの」と「無視できるもの」にあらかじめ区分する必要がある** ::: 先ほど「全ての因子間に交互作用が存在すると仮定した場合」と述べましたが、現実的には3因子以上の交互作用は存在しないことが多く、直交表を用いた実験においては無視する（誤差としてまとめる）ことがほとんどです。 そのため、考慮すべきなのは2因子間の交互作用になるわけですが、それが全ての因子間で存在すると仮定してしまうと結局実験回数が増えてしまいます。 ゆえに、技術的側面などを考慮して、**あらかじめ取り上げる交互作用とそうでないものを区分しておく**ことで、より少ない回数の実験で解析を進めていくことができます。 ![実験計画法　交互作用の選定](/wp-content/uploads/2023/10/DOE_interaction_selection.jpg) ### 3. 適用例 ある工場では、コンクリート製のタイルを製造しています。 しかし、タイルの強度にはばらつきがあることが判明しました。 そこで、タイルの強度を最も高める製造条件を探すため、実験計画法（実験計画&分散分析）を適用します。 #### 因子と水準の決定 今回の例では、以下の因子と水準を取り上げます。 :::box **・成型圧力(A)\:A1 A2**\ **・成型時間(B)\:B1 B2**\ **・乾燥温度(C)\:C1 C2**\ **・乾燥時間(D)\:D1 D2**\ **・熱処理温度(F)\:F1 F2**\ **・熱処理時間(G)\:G1 G2**\ **・冷却温度(H)\:H1 H2** ::: また、技術的側面を考慮し、A×B、A×C、B×C、C×D、D×Fの5つの交互作用を取り上げることにしました。 さらに、特性値は**タイルの強度**で、大きいほどよいとします。 :::talk{who="robotan"} なんでアルファベットでEは使われてないの？ ::: :::talk{who="umatan"} 統計学の分野では、Eは誤差を表すものだから！ごっちゃにならないようにしているんだ！ ::: #### 因子の割り付け 続いて取り上げた因子と交互作用を直交表に割り付けていきます。 因子数は7、交互作用は5の合わせて12であるので、**L16直交表**を用いることにします。 ![実験計画法　L16直交表](/wp-content/uploads/2023/10/DOE_L16array.jpg) ただし、ここに適当に因子や交互作用を割り付けてはいけません。 **線点図**を用いて割り付けを決定していきます。 線点図とは、**因子を点**で表し、2点を結んだ**線分を交互作用**に対応させたグラフです。 また、実験者が取り上げた因子と考慮する交互作用を表現したものを**必要な線点図**、それぞれの直交表に応じて用意されているものを**用意されている線点図**と呼びます。 今回の事例における、必要な線点図と用意されている線点図は以下のようになります。 また、L16直交表について用意されている線点図は他にもいくつか存在します。 ![実験計画法　必要な線点図](/wp-content/uploads/2023/10/DOE_diagram_needed.jpg) ![実験計画法　用意されている線点図](/wp-content/uploads/2023/10/DOE_diagram_supply.jpg) :::talk{who="robotan"} なんで適当に割り付けちゃダメなの？ ::: :::talk{who="umatan"} 実は任意の2列間の交互作用が現れる列は元々決まっているんだ！だから、その列に他の因子を割り付けてしまうと、その因子の主効果と交互作用が混ざってしまうんだ！ ::: :::talk{who="robotan"} 本当に見たい因子の効果がわからなくなってしまうってことか！ ::: 線点図を用いて割り付けることで、因子の主効果と交互作用が混ざってしまうことを防ぐことができます。 例えば、用意されている線点図の一番下は、2と4が線分6によって結ばれています。 これは、**直交表の第2列と第4列に割り付けた因子の交互作用が第6列に現れる**ことを示しています。 必要な線点図を用意されている線点図に組み込むと、以下のような線点図になります。 ![実験計画法　線点図　割り付け](/wp-content/uploads/2023/10/DOE_diagram_set.jpg) この線点図をもとに**L16直交表**への割り付けとデータの測定（実験）を行った結果、以下のような表が得られました。 なお、実験順序はランダムです（無作為化の原則）。 ![実験計画法　L16直交表　割り付け後](/wp-content/uploads/2023/10/DOE_L16array_set.jpg) また、直交表では、因子が割り付けられておらず、検討すべき交互作用も現れない列を誤差列と呼びます。 #### 平方和の計算 [分散分析](/analysis-of-variance)における一元配置法や二元配置法と同様、平方和を計算していきます。 ただし、**2水準系の場合にはより簡単な形で計算することができます。** それが、直交表特有の**列平方和**と呼ばれるものです。 第\[$k$]列が第$i$水準のデータ和を$T_{[k]_{i}}$、第\[$k$]列が第$i$水準のデータ数を$N_{[k]_{i}}$、第\[$k$]列が第$i$水準のデータの平均を$\bar x_{[k]_{i}}$とすると、第\[$k$]列に対する列平方和は以下のように定義できます。 :::box （第\[$k$]列の列平方和）=$\displaystyle S_{[k]}=\frac{(T_{[k]_{1}}-T_{[k]_{2}})^2}{N}$ ::: 今回の例であれば、因子Cを第1列に割り付けていますから、$T_{C_{1}}=T_{[1]_{1}}$、$T_{C_{2}}=T_{[1]_{2}}$であり、$S_{C}=S_{1}$となります。 つまり、直交表を用いた解析では、**因子の平方和はその因子を割り付けた列の列平方和と対応する**と言えます。 :::talk{who="umatan"} 誤差列の平方和は誤差平方和に対応するよ！ ::: また、列平方和に対して**列自由度**を考えることができます。 それぞれの列には２水準が設定されているため、以下のように定義できます。 :::box （第\[$k$]列の列自由度）=$\phi_{[k]}=2-1=1$ ::: さらに、L16直交表では、以下の2つの式が成り立ちます。 :::box 総平方和： $S_{T}=S_{[1]}+S_{[2]}+S_{[3]}+\cdots+S_{[15]}$ 総自由度： $\phi_{T}=\phi_{[1]}+\phi_{[2]}+\phi_{[3]}+\cdots+\phi_{[15]}$ ::: このことから、２水準系直交表では、**総平方和$S_{T}$を自由度が1ずつの平方和に分解している**と言えます。 :::talk{who="robotan"} ばらつきを分解するっていう基本はやっぱり変わらないんだな！ ::: #### 分散分析表の作成 さて、平方和の算出方法がわかったので、それぞれの主効果および交互作用、誤差の平方和を計算し、**分散分析表**を作成していきます。 作成した**分散分析表**がこちら。 ![実験計画法　分散分析表](/wp-content/uploads/2023/10/DOE_anova_table.jpg) $F(1,3;0.05)=10.1$、$F(1,3;0.01)=34.1$より、**因子A、D**と**交互作用C×D**が有意であることがわかりました。 有意差については以下の記事で解説しています！ #### プーリング 分散分析表を見ると、$F_{0}$が小さい、つまり効果がない因子や交互作用がいくつかあることがわかります。 そこで、それらを誤差と見なし、再度分散分析表を作成していきます。 この作業のことを**プーリング**と呼びます。 プーリングを行うことで、誤差の自由度が増え、誤差分散の推定精度が上がります。 プーリングの絶対的な基準は存在しませんが、慣例的に用いられているルールには以下のようなものがあります。 :::box **① $F_{0}\leq 1$のときは誤差項にプールする** **② $1それぞれの2因子の2元表を作成することで確かめることができます。 例えば因子AとBの2元表は以下のようになります。 ![実験計画法　AB二元表](/wp-content/uploads/2023/10/DOE_ABtable.jpg) 表より、A1B2の平均値が最も大きな値であることがわかります。 よって、**因子A、Bの最適水準はA1B2**となります。 同様にして、**因子C、DはC1D2**、**因子FはF2**が最適水準となります。 ## 実験計画法をPythonでやってみよう ![stories pc](/wp-content/uploads/2021/01/All-the-data-bro-e1613475965440-1024x789.png) ここからは、Pythonで実際に実験計画法をやっていきます（実務では解析ソフトなどが使われることが多いです）！ 適用例でも紹介した、**直交表を用いた実験計画法**をやっていきましょう。 分析の目的は、最もよい水準を見つけること！ 以下の流れで分析を進めていきます。 :::box **1. ライブラリのインストール\ 2\. データの準備\ 3\. 列平方和と自由度の算出\ 4\. 分散分析表の作成\ 5\. 最適水準の決定** ::: ### 1. ライブラリのインストール 今回はscipyのstatsを使って分散分析（$F$値の算出）を行っていきます。 分析で使用する他のモジュール等もここで全てimportしておきます。 ```python import pandas as pd import numpy as np from scipy.stats import f ``` ### 2. データの準備 今回は因子A\~D（各2水準）と交互作用A×Bを取り上げます。 そのため、**L8**直交表を用いることにします。 直交表への割り付けとデータの測定を既に終えたとして、**直交表を表示**します。 なお、今回のデータは値が大きいほど良いものとします。 ```python #因子A~D #因子は全て2水準ずつ #技術的観点から、交互作用はA×Bのみ存在するとした #L8直交表への割り付け A = np.array([1,1,1,1,2,2,2,2]) #因子A B = np.array([1,1,2,2,1,1,2,2]) #因子B A_B = np.array([1,1,2,2,2,2,1,1]) #交互作用A×B C = np.array([1,2,1,2,1,2,1,2]) #因子C D = np.array([1,2,1,2,2,1,2,1]) #因子D E_1 = np.array([1,2,2,1,1,2,2,1]) #誤差列 E_2 = np.array([1,2,2,1,2,1,1,2]) #誤差列 #直交表をデータフレームとして表示 L8_array = pd.DataFrame({ 'A':A, 'B':B, 'A-B':A_B, 'C':C, 'D':D, 'E1':E_1, 'E2':E_2, 'データ':np.array([65,89,76,91,33,37,48,43]) }) #見やすさのため、行番号を1からに変更 L8_array.index = [i for i in range(1,len(L8_array)+1)] L8_array ``` ![実験計画法　L8　割り付け後](/wp-content/uploads/2023/10/DOE_python_L8.jpg) ### 3. 列平方和と自由度の算出 続いて、列平方和と自由度をそれぞれ算出しておきます。 また、この後扱いやすくするため、それぞれリストにまとめておきます。 ```python #各列の平方和を計算 S_A = (sum(L8_array.query('A==1')['データ'].values)-sum(L8_array.query('A==2')['データ'].values))**2/8 S_B = (sum(L8_array.query('B==1')['データ'].values)-sum(L8_array.query('B==2')['データ'].values))**2/8 S_AB = (sum(L8_array.query('A-B==1')['データ'].values)-sum(L8_array.query('A-B==2')['データ'].values))**2/8 S_C = (sum(L8_array.query('C==1')['データ'].values)-sum(L8_array.query('C==2')['データ'].values))**2/8 S_D = (sum(L8_array.query('D==1')['データ'].values)-sum(L8_array.query('D==2')['データ'].values))**2/8 S_E1 = (sum(L8_array.query('E1==1')['データ'].values)-sum(L8_array.query('E1==2')['データ'].values))**2/8 S_E2 = (sum(L8_array.query('E2==1')['データ'].values)-sum(L8_array.query('E2==2')['データ'].values))**2/8 S_E = S_E1 + S_E2 S_T = sum([S_A, S_B, S_C, S_D, S_AB, S_E]) #総平方和 S_list = [S_A, S_B, S_C, S_D, S_AB, S_E, S_T] #リストにまとめておく #各列の自由度をリストにまとめておく deg_list = [1,1,1,1,1,2,7] ``` ### 4. 分散分析表の作成 続いて、算出した列平方和および自由度をもとに、分散分析表を作成していきます。 さらに有意であるかどうかを判断するため、F分布の5%点および1%点を表示しておきます。 ```python #分散の計算 V = [x / y for (x, y) in zip(S_list,deg_list)] #F0の計算 F0 = [x/V[5] for x in V] #分散分析表の作成 ANOVA_table = pd.DataFrame({ 'S':S_list, 'Φ':deg_list, 'V':V, 'F0':F0 }) #行名を要因名に変更 index_name = ['A','B','C','D','A×B','E','T'] ANOVA_table.index = index_name #分散分析表として表示するため、不要のセルには欠損値を格納 ANOVA_table.loc['T','V'] = None ANOVA_table.loc['E','F0'] = None ANOVA_table.loc['T','F0'] = None ANOVA_table ``` ```python #分子、分母の自由度 dfn, dfd = 1, 2 #5%点(上側なので、95% を指定) five = f.ppf(0.95, dfn, dfd) print('上側確率 5%:', five) #1%点(上側なので、99% を指定) one = f.ppf(0.99, dfn, dfd) print('上側確率 1%:', one) ``` ![実験計画法　分散分析表　python](/wp-content/uploads/2023/10/DOE_python_ANOVA_table.jpg) ![実験計画法　F分布](/wp-content/uploads/2023/10/DOE_f_dist.jpg) よって、**因子Aは高度に有意**であり、**交互作用A×Bも有意**であることがわかりました。 なお、今回はプーリングは行わないこととします。 ### 5. 最適水準の決定 最後に、最適水準の決定を行います。 有意であるのは因子Aと交互作用A×Bであったため、因子AとBの最適水準を求めます。 上の適用例で説明したように、AB二元表を作成し、最適水準を決定します。 ```python #優位となったのは因子A、交互作用A×B #交互作用が有意となったので因子Bも取り上げる #各水準組み合わせのときのデータの平均値を算出 Mean_A1B1 = sum(L8_array.query('A==1 & B==1')['データ'].values)/2 #A1とB1 Mean_A1B2 = sum(L8_array.query('A==1 & B==2')['データ'].values)/2 #A1とB2 Mean_A2B1 = sum(L8_array.query('A==2 & B==1')['データ'].values)/2 #A2とB1 Mean_A2B2 = sum(L8_array.query('A==2 & B==2')['データ'].values)/2 #A2とB2 #AB二元表の作成 AB_table = pd.DataFrame( np.array([[Mean_A1B1,Mean_A1B2], [Mean_A2B1,Mean_A2B2]]), index = ['A1','A2'], columns = ['B1','B2'] ) AB_table ``` ![実験計画法　python AB二元表](/wp-content/uploads/2023/10/DOE_python_ABtable.jpg) AB二元表より、A1B2の水準組み合わせのとき、最もデータの平均値が大きくなっていることがわかります。 よって、**最適水準はA1B2**であると言えます。 ## 実験計画法　まとめ ここまでご覧いただきありがとうございました！ 本記事では、実験計画法の概要からPythonでの実装方法まで解説していきました。 実験計画法の主な使用場所は品質管理の現場ですが、この考え方は様々な場所で応用することができます。 複数の要因の効果を効率よく調べたい場合は、ぜひ**実験計画法**を使ってみましょう！ また、今回分析に用いたPythonについて勉強したい方は以下の記事を参考にしてみてください！ さらに詳しくAIやデータサイエンスの勉強がしたい！という方は当サイト「スタビジ」が提供する[スタビジアカデミー](https://toukei-lab.com/achademy/)というサービスで体系的に学ぶことが可能ですので是非参考にしてみてください！ :::box **AIデータサイエンス特化スクール「スタアカ」** ![スタアカ（スタビジアカデミー）のトップページ](/wp-content/uploads/2026/10/achademy-top.jpg) :::
【価格】 ライトプラン：1,280円/月（全コースが受講し放題）
プレミアムプラン：129,800円（通常149,800円・買い切り）
【オススメ度】
【サポート体制】
【受講形式】 オンライン形式
【学習範囲】 統計学・機械学習・ディープラーニングを網羅的に学ぶ
生成AI（LLM・RAG・プロンプトエンジニアリング）を学ぶ
SQLとPythonを組み合わせて実データを使った様々なワークを行う
実践的なビジネスフレームワークを学ぶ
マーケティングの実行プラン策定
データサイエンティストとしての自分の経験をふまえてエッセンスを詰め込んだのがこちらのスタビジアカデミー、略して「[スタアカ](https://toukei-lab.com/achademy/)」！！ 当メディアが運営するスクールです。 月額1,280円のライトプランで、統計学から機械学習・生成AIまで全コースが受講し放題です！ プレミアムプランでは、24時間以内の質問対応と現役データサイエンティストとの1on1（30分×4回）を実施します！ カリキュラム自体は、他のスクールと比較して圧倒的に良い自信があるのでぜひ受講してみてください！ 他のスクールのカリキュラムはPythonでの機械学習実装だけに焦点が当たっているものが多く、実務に即した内容になっていないものが多いです。 そんな課題感に対して、実務で使うことの多いSQLや機械学習のビジネス導入プロセスの理解なども合わせて学べるボリューム満点のコースになっています！ Pythonが初めての人でも学べるようなカリキュラムにしておりますので是非チェックしてみてください！ ウォルマートのデータを使って商品の予測分析をしたり、実務で使うことの多いGoogleプロダクトのBigQueryを使って投球分析をしたり、データサイエンティストに必要なビジネス・マーケティングの基礎を学んでマーケティングプランを作ってもらったりする盛りだくさんの内容になってます！ :::box **・BigQuery上でSQL、Google Colab上でPythonを使い野球の投球分析**\ **・世界最大手小売企業のウォルマートの実データを用いた需要予測**\ **・ビジネス・マーケティングの基礎を学んで実際の企業を題材にしたマーケティングプランの策定**\ **・生成AIのAPIを使ってチャットボットやAIエージェントを作る** ::: ::button[スタアカの受講はこちら]{href="https://toukei-lab.com/achademy/" note="＼今すぐ試す／"}