[論文レビュー] Optimal Designs for Two-Level Factorial Experiments with Binary Response

Abhyuday Mandal, Jie Yang+1 more|arXiv (Cornell University)|Mar 8, 2010

Optimal Experimental Design Methods被引用数 3

ひとこと要約

本稿は、一般化線形モデルを用いた二水準要因実験における二値応答のための局所的D最適設計を開発する。2²設計における特殊な場合について解析的解を導出し、一般の2ᵏ実験において一様設計が非常に効率的かつ最大最小D最適であることを示し、スクリーニング応用におけるロバストネスと実用的有用性を提供する。

ABSTRACT

We consider the problem of obtaining locally D-optimal designs for factorial experiments with qualitative factors at two levels each with binary response. Our focus is primarily on the 2^2 experiment. In this paper, we derive analytic results for some special cases and indicate how to handle the general case. The performance of the uniform design in examined and we show that this design is highly efficient in general. For the general 2^k case we show that the uniform design has a maximin property.

研究の動機と目的

• 従来の線形モデルが適用できない二値応答を伴う二水準要因実験の最適設計理論を確立すること。
• 未知パラメータに依存する情報行列を伴う一般化線形モデルにおけるD最適性の課題に取り組むこと、特に二値応答の下でのD最適性。
• 2²および2ᵏ要因実験のさまざまな設定において、一様設計の性能とロバストネスを評価すること。
• 特殊な場合の解析的解と、一般の場合の計算戦略（例：円筒代数的分解）を提供すること。
• 一般の2ᵏ設定において一様設計が最大最小D最適性を達成することを示すことにより、パrameterの不確実性下でもロバストであることを保証すること。

提案手法

• Chernoff（1953）の局所最適性アプローチを用い、未知パラメータを仮定値に置き換えてD最適性基準を定義する。
• 一般化線形モデルにおけるフィッシャー情報から導かれる漸近的情報行列の行列式を最大化するD基準を適用する。
• 解析的解が得られない一般の2²設計における非凸最適化問題を解くために、円筒代数的分解（CAD）を用いる。
• 一般化線形モデルフレームワーク内で一般的なリンク関数（ログット、プロビット、ログログ、補完対数ログ）を考慮する。
• 最大最小基準を用いて、さまざまなパrameter設定におけるD効率性とロバストネスを評価することで、一様設計を分析する。
• 不等式（例：AM-GM）を用いてD最適性基準の上限を導出し、ラグランジュおよび摂動解析を用いて最適性条件を証明する。

実験結果

リサーチクエスチョン

• RQ1 主効果モデル下での二値応答を伴う2²要因実験における局所的D最適設計の解析的形は何か？
• RQ2 さまざまなパrameter設定において、一様設計のD効率性は他の設計と比較してどのように異なるか？
• RQ3 一般の2ᵏ要因実験設定において、一様設計が最大最小D最適であることを示せるか？
• RQ4 パラメーターベクトルにどのような条件が満たされると、一様設計が最適または近似的に最適となるか？
• RQ5 解析的解が得られない場合、特に2²ケースにおいて最適設計をどのように計算できるか？

主な発見

• 主効果を伴う2²設計において、主効果の分散がバランスしているとき、一様設計のD効率は少なくとも1 - (3/4) × 2^{1/3} ≈ 0.442に達する。
• 主効果の分散v₁ ≥ v₂ + v₃ + v₄のとき、最適設計は非一様であり、p₁ = 0 かつ p₂ = p₃ = p₄ = 1/3 となる。このときL = v₁/27を達成する。
• 一般の2ᵏ実験において、一様設計は最大最小D最適である。これは、すべてのパrameter値に対してD基準の下限を最大にするという意味である。
• 一様設計は一般に非常に効率的であり、真のパラメータ値が仮定値から逸脱しても、ロバストな性能を示す。
• v₁ = v₂ + v₃ + v₄のとき、一様設計のD効率は正確に1 - (3/4) × 2^{1/3} ≈ 0.442であり、最大最小基準下で達成可能な最小値である。
• 摂動解析により、v₁ < v₂ + v₃ + v₄のとき、p₁をわずかに増加させ、p₂, p₃, p₄を比例して減少させるとD基準が向上することが示され、この領域では一様性が最適でないことが示唆される。