2016-03-25 11 views
1

Was ist das assign Attribut eines linearen Modells passen? Es soll irgendwie die Position des Response-Terms angeben, aber in der Praxis scheint es alle Koeffizienten im Modell aufzulisten. Es ist mein Verständnis, dass assign ist ein Übertrag von S und es wird nicht von glm() unterstützt. Ich muss die äquivalente Information für glm extrahieren, aber ich verstehe nicht, was die Implementierung für lm tut und scheint nicht, den Quellcode auch zu finden. Die Hilfedatei für lm.fit sagt wenig hilfreich:

Nicht-Null-fits Komponenten assign haben, effects und qr auf die lineare Anpassung, für die durch Extraktor Funktionen Arzneimittel (es sei denn, nicht angefordert) wie summary und effects

Antwort

5

Sie finden diese in help("model.matrix"), die diese Werte schafft:

Es ist ein Attribut "zuweisen", eine integ Er Vektor mit einem Eintrag für jede Spalte in der Matrix geben den Begriff in der Formel, die Aufstieg zur Spalte gab. Wert 0 entspricht dem Achsenabschnitt (falls vorhanden) und positive Werte zu Termen in der Reihenfolge, die durch den Term angegeben wird. Labels Attribut der dem Objekt entsprechenden Termstruktur.

Also, es bildet die Design-Matrix auf die Formel.

+0

Kurz und süß! Sehr hilfreich!!!!! – rimorob

+0

lächerlich, dass dies nicht _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ ~ _ _ – MichaelChirico

2

Die Zahlen von $assign repräsentieren die entsprechende Prädiktorvariable. Wenn Ihr Prädiktor kategorisch mit 3 Stufen ist, sehen Sie die entsprechende Nummer (3-1) in Ihrem $assign Aufruf. Beispiel:

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) 

Sie werden die quantitativen Prädiktoren nur einen Wert sehen müssen (hwy im Beispiel oben), da sie von einem Begriff in der Design Formel dargestellt werden.

+0

Wo ist der Code, der dieses Feld berechnet? Ich konnte es bisher nicht finden. – rimorob

+0

Ah, ich sehe das Problem jetzt. Das ist seltsam, dass der '$ assign' nicht im Quellcode für' lm' oder 'glm' enthalten ist. Ich habe den Grund noch nicht gefunden, aber ich werde weiter suchen. – cgage

+1

Suchen Sie im Quellcode von 'stats ::: lm.fit'. – Roland