2016-03-25 11 views
1

線形モデルフィットのassign属性は何ですか?それは何とか応答項の位置を提供するはずですが、実際にはモデルのすべての係数を列挙しているようです。 assignはSからのキャリーオーバーであり、それはglm()によってサポートされていないと私は理解しています。私はglmの同等の情報を抽出する必要がありますが、実装がlmのために何をするのか分かりませんし、ソースコードを見つけることもできません。 lm.fitのヘルプファイルはunhelpfully、言う:

null以外のフィットは、このようなsummaryeffects

として抽出機能で使用するための線形フィットに関連するコンポーネント assigneffectsと(要求されていない限り) qrを、持っています

答えて

5

あなたはこれらの値を作成する、help("model.matrix")でこれを見つける:

属性があり、INTEG "を割り当てます"配列内の各列が のエントリを持つ式のベクトルを使用して、 を列に付けた式の項を与えます。値0はインターセプト(存在する場合)に対応し、 という用語は、オブジェクトに対応する用語構造のterm.labels 属性によって与えられた順序で正の値に対応する。

したがって、デザインマトリックスを数式にマッピングします。

+0

簡単で甘い!非常に役立ちます!!!!! – rimorob

+0

これは言及されていないとばかげているのですが、 '?lm'の' all at_ ... 'は' assign 'が最初から出てきた場所を20分追跡しなければなりませんでした。 – MichaelChirico

2

$assignの数字は、対応する予測変数を表します。プレディクタが3つのレベルのカテゴリに分類されている場合は、$assignコールに対応する番号(3-1)が表示されます。例:

data(mpg, package = "ggplot2") 
m = lm(cty ~ hwy + class,data = mpg) 
m$assign 
    [1] 0 1 2 2 2 2 2 2 
# Note how there is six 2's to represent the indicator variables 
# for the various 'class' levels. (class has 7 levels) 

あなたはそれらが設計式で1つの項によって表されるので、定量的な予測因子でのみ、一つの値(上記の例ではhwy)を有することになるでしょう。

+0

このフィールドを計算するコードはどこですか?私はこれまでそれを見つけることができませんでした。 – rimorob

+0

ああ、私は今問題を見る。これは '$ assign'が' lm'や 'glm'のソースコードにないことは奇妙です。私はなぜまだ見つけていないが、私は見続けるだろう。 – cgage

+1

'stats ::: lm.fit'のソースコードを見てください。 – Roland

関連する問題