7因子を二水準で全組合せにすると128 runです。資源不足だからと任意の16行を選べば、均衡と効果分離が崩れます。部分要因設計は無作為な省略ではなく、どの効果を一緒に見るかを知って選ぶ構造的節約です。
この単元の問いは一つです。
run数を減らすと、どの効果がデータで区別不能になるか。
2³の半分を構造的に選びます
AとBの四組合せを作り C=AB と置けば、4-run 2^(3−1) half fractionになります。このgeneratorの両辺へCを掛けると defining relation I=ABC を得ます。
設計表でC列とAB列は完全に同じです。応答がその列とともに動いても、原因がC主効果かAB交互作用かは区別できません。alias chainは次です。
- I = ABC
- A = BC
- B = AC
- C = AB
「A効果推定」と書かれた数は実際にはAとBCを合わせたcontrastです。BCが0に近いという仮定があって初めて主にAとして解釈できます。
aliasは高い相関より強い同一性です
完全aliasでは二つのmodel columnが同じか符号だけ反対で、同じデータから別々の係数を推定できません。非正規・optimal設計では完全同一でなくても推定量が相関し精度が下がることがあります。
ソフトウェアが係数一つを出しても、すべての候補効果を推定したことにはなりません。実行前にalias matrix、estimability、correlation of estimatesを確認します。
経済性はeffect sparsityと高次交互作用が小さいという仮定に依存します。aliasされた交互作用が大きければ主効果へ誤って帰属します。重要な相互作用候補はfactor assignment、高いresolution、foldover、後続実験で分離してください。
resolutionは最も深刻なalias次数を要約します
defining relationでidentityを除く最短word長がresolutionです。
- Resolution III: 主効果が二因子交互作用とalias
- Resolution IV: 主効果は二因子交互作用と分離、二因子同士はalias
- Resolution V: 主効果と二因子交互作用は分離、二因子は三因子とalias
高いresolutionには一般に多くのrunが必要です。数が高いことが常に最善ではなく、多因子の高速screeningか特定interactionの推定かで必要な構造が変わります。完全要因設計はeffect aliasの意味でresolution infinityとも呼ばれます。block作成では高次interactionとblock effectがconfoundされ得ます。
fraction選択には隠れた仮定があります
実用性は有用だが法則ではない次の原則に依ります。
- sparsity: 多数候補のうち実際に大きい効果は少数
- hierarchy: 高次interactionが重要なら関連する低次項も考慮
- heredity: interactionは関連主効果と現れやすい
主効果が小さくてもcrossover interactionは大きくなり得るため、heredityを盲信しません。機序上重要な組合せは設計でclearにします。
In-Silico Lab: 同じ列を直接確認します
- 四runでCとAB列が全行一致するか見ます。
- defining relationへA、B、Cを掛け、各alias chainを導きます。
- A contrastが8なら、A=8かA+BC=8か答えます。
- 八run完全設計へ広げるとどの列が分離するか考えます。
4-run 반분할의 alias를 직접 읽으세요
2³ 완전설계의 절반에서 C=AB를 generator로 선택하면 어떤 열들이 완전히 같아지는지 확인합니다.
처음이라면: 무엇을 눌러야 하나요?
- 1. 질문을 먼저 읽기Lab 제목에서 이번에 비교할 한 가지를 확인합니다.
- 2. 조건 하나만 바꾸기처음에는 n, 효과, 산포 같은 입력 중 하나만 바꾸십시오.
- 3. 새 합성 표본 누르기새 합성 데이터가 만들어집니다. 같은 조건도 표본에 따라 달라질 수 있습니다.
- 4. 그림과 계산 결과 비교하기바꾸기 전후 무엇이 움직이고 무엇이 그대로인지 한 문장으로 적어보십시오.
막히면 초기화로 돌아가 기본 결과를 본 뒤 조건 하나만 바꾸십시오. 이 Lab은 정답 판정기가 아니라 패턴 관찰 도구입니다.
같은 설정의 합성 관측
| run | A | B | C=AB | ABC |
|---|---|---|---|---|
| 1 | -1 | -1 | 1 | 1 |
| 2 | -1 | 1 | -1 | 1 |
| 3 | 1 | -1 | -1 | 1 |
| 4 | 1 | 1 | 1 | 1 |
계산 결과
I = ABCA = BCB = ACC = ABA 추정값이 커도 원인이 A인지 BC인지 이 네 run만으로 분리할 수 없습니다. 작은 고차 교호작용이라는 가정이 절약의 가격입니다.
교육용 synthetic model · bjs-factorial-sequence-v1. 실제 연구 판단에는 실험단위, 결측, 분포, 다중성, 사전계획과 도메인 기준을 별도로 반영해야 합니다.
JMP Design Evaluationは実行前に読みます
각 효과 추정에 섞이는 다른 효과를 명시합니다.
resolution과 추정 가능성을 실행 전에 확인합니다.
0이 아닌 상관은 효과 분리가 약해졌다는 신호입니다.
Alias Tableは結果解釈後の付録でなく、設計選択前のGateです。予想されるactive interactionが主効果とaliasなら、fractionまたはfactor assignmentを変えます。候補modelがestimableかを確認することが核心です。
部分要因設計は逐次実験の第一波です
screening fractionで有力因子を見つけても、同じデータで最終最適点を確定しません。次の段階は状況により、反対fractionを加えるfoldover、aliasを切る選択run、重要因子だけの完全要因設計、中心点、反応曲面設計です。第一波のnull/active判定に全決定を賭けず、後続実験で曖昧なaliasを解きます。
結果文の例
5因子の
2^(5−1)Resolution V half fractionを用いた。事前候補modelの全主効果と二因子交互作用は相互にclearで、各二因子項は三因子以上とaliasだった。AとABが有力に見えたが、sparsity仮定に依存するため、独立confirmationと選択runで構造を確認するよう計画した。
Resolution III screeningなら「A効果」でなく「A+BC alias contrast」と正直に書きます。
この単元の要点
- fractionは完全設計の構造的な一部です。
- generatorからdefining relationとalias chainが出ます。
- aliasされた効果は同じデータで分離推定できません。
- resolutionは最も深刻な効果混同次数を要約します。
- sparsity・hierarchy・heredityは有用な仮定で保証ではありません。
- 部分要因結果のaliasは後続実験で解消します。
この連続執筆は、要請された自動停止線U20で止まります。screening後続実験、中心点、曲率は候補ですが、別の再開前に執筆してはいけません。
公式補足資料
- NIST/SEMATECH · Constructing a 2^(3−1) design
- NIST/SEMATECH · Confounding and aliasing
- NIST/SEMATECH · Design resolution
この記事の設計とLabは教育用の合成資料です。