1

私は、7 年間の時間傾向変数 (つまり、2000=1,2001=2,...,2006=7) と、6 年間のダミー変数 (それぞれに 1 つのバイナリ変数) を含むモデルを持っています。 2000年を除く年）。R にこの線形モデルを適合させると、次のようになります。

olsmodel=lm(lnyield ~ lnx1+ lnx2+ lnx3+ lnx4+ lnx5+ x6+ x7+ x8+ timetrend+ 
                     yeardummy2001+ yeardummy2002+ yeardummy2003+ yeardummy2004+ 
                     yeardummy2005+ yeardummy2006)

モデル サマリーの最後のダミー変数に対して生成された NA を取得します。次の「係数: (特異点のため 1 は定義されていません)」とともに。

すべての x_i 変数が連続的であり、ダミーのサブセットがなく、時間傾向が互いの線形結合であるため、なぜこれが起こっているのかわかりません。

なぜこれが起こっているのかについての助けをいただければ幸いです!

4

1 に答える 1

1

問題は、年の傾向を に設定し、1:n各年のダミー変数も含めると、たまたま非フル列ランクの共変量行列が生成されることです。

r1、r2、r3の3つのカテゴリしかない場合、モデルは次のとおりでy ~ trend + c2 + c3あり、共変量行列は次のとおりです。

> mat
     int trend c2 c3
[1,]   1     1  0  0
[2,]   1     1  0  0
[3,]   1     2  1  0
[4,]   1     2  1  0
[5,]   1     3  0  1
[6,]   1     3  0  1

matまた、共変量行列の列ランクは、推定に必要な係数の数 (4) ではなく、わずか 3 であることがわかります。つまりt(mat)%*%mat、特異です。そのため、エラーが発生する可能性があります。

于 2012-11-19T21:21:50.640 に答える