2017-07-24 14 views
0

Wie können Sie die Referenzkategorie drucken, die verwendet wird, wenn eine kategoriale/nominale Variable in ein lineares Modell eingegeben wird? Hier ein Beispiel:Drucken (anzeigen) Referenzkategorie in R's lm Zusammenfassung?

summary(lm(data = iris, Sepal.Length ~ Species))

Coefficients: 
        Estimate Std. Error t value Pr(>|t|)  
(Intercept)   5.0060  0.0728 68.76 < 2e-16 *** 
Speciesversicolor 0.9300  0.1030 9.03 8.8e-16 *** 
Speciesvirginica 1.5820  0.1030 15.37 < 2e-16 *** 

Hier ist, was Ich mag würde:

    Estimate Std. Error t value Pr(>|t|)  
(Intercept)   5.0060  0.0728 68.76 < 2e-16 *** 
**Reference: Speciessetosa** 
Speciesversicolor 0.9300  0.1030 9.03 8.8e-16 *** 
Speciesvirginica 1.5820  0.1030 15.37 < 2e-16 *** 

Wenn es einen Weg gibt, diese Arbeit zu machen im allgemeinen (wenn es mehrere kategorische Prädiktoren, dann wird jeder Referenzgruppe ist leicht identifizierbar), das wäre am besten. Und wenn es eine Möglichkeit gibt, die Formatierung besonders deutlich zu machen, wäre das doppelt hervorragend (ich bin nicht der obigen Beispielformatierung verpflichtet).

Antwort

0

Sie können angeben, dass Ihr Achsenabschnitt Null ist.

summary(lm(Sepal.Length ~ Species + 0, data = iris)) 
#Coefficients: 
#     Estimate Std. Error t value Pr(>|t|)  
#Speciessetosa  5.0060  0.0728 68.76 <2e-16 *** 
#Speciesversicolor 5.9360  0.0728 81.54 <2e-16 *** 
#Speciesvirginica 6.5880  0.0728 90.49 <2e-16 *** 
+0

Danke, aber das ändert die Ergebnisse. Ich möchte die Referenzgruppe NAME anzeigen, nicht die Ergebnisse ohne einen Intercept. – dca

Verwandte Themen