質問: 分散分析 4因子3水準で実験する場合、l9直交表だと誤差列がないのですが、 分散分析は可能ですか? この場合の対処法教えて下さい。 そもそもl9ではダメでしょうか?
回答: l9直交表は4列から成り、各列に4因子を割付ければ各因子の主効果に相当する分散を求めることができます。すなわち4つの因子a,b,c,dの相対的な効果の大きさを見積もることができます。これを分散分析と呼ぶかどうか? 一般的に分散分析は要因の分散と誤差の分散との比(分散比)を求め、この値をf検定して有意判定を行い、その結果をもとに変数(要因)の選択を行うことを目的としています。その意味では分散分析とは呼べないかもしれません。しかし、4因子の効果の最も小さいものを誤差とみなせば分散比を求めることも可能です。これが対処法の1つでしょう。一般的には誤差のプーリングと呼ばれます。 ここからが「そもそも論」になりますが、l9の実験回数は9回で自由度は9-1=8です。ご存じのように全平均を計算するために自由度1を使用します。そして4因子を評価することから誤差の自由度は8-4=4とすることができるはずです。しかしl9で4因子を評価すると誤差列がなく誤差が評価できない!何故でしょうか? それは3水準で実験することで各因子の自由度が2となるためです。すなわち主効果(1次の効果)と非線形効果(2次の効果)を評価することを目的として3水準とするのです。3水準の実験を計画することは2次の効果を評価することが前提となります。 ということで「そもそも」2次の効果を評価する目的があるか否か、がl9(3水準)選択の分かれ道になりそうです。もちろん3水準でも1次の効果だけ評価することは可能ですが、実験計画の基本からは外れているかもしれません。