2017-07-16 1 views
0

Ich versuche herauszufinden, warum meine lm() Schätzungen sind anders als die geom_smooth für die gleichen Daten und Formel. Insbesondere ist die Steigung für meine Gruppierungsvariable "cat" level 5> 0 in dem lm() Ausgang, aber < 0 in dem geom_smooth (und daher scheint das Diagramm nicht die Zusammenfassungstabelle zu reflektieren).Verschiedene geom_smooth und lm() Schätzungen in R: vergessen, Basis Steigung und Interaktion Steigung zusammen

Hier sind the data. (Einfacher als kommen mit Beispieldaten, die sich ähnlich verhalten.)

Das Modell: summary(lm(data=df, y~x*cat))

Hinweis für die Steigung x:cat5> 0 ist.

Das Grundstück:

library(ggplot2) 
plt <- ggplot(df, aes(x=x, y=y, group=cat)) + 
    geom_smooth(method="lm", show.legend=FALSE) + 
    facet_wrap(~cat, nrow=1) + 
    geom_point(aes(color=color) 

die geom_smooth Schätzungen abrufen (nach @Pedro Aphalo Antwort here):

library(ggpmisc)  
my.formula <- y~x 
plt + stat_poly_eq(formula = my.formula, 
      aes(label = paste(..eq.label.., ..rr.label.., sep = "~~~")), 
      parse = TRUE) 

Beachten Sie die Steigung in Facette 5 < 0. Sind lm() und geom_smooth ein mit verschiedene Summe von Quadraten oder etwas? Welche Version melde ich in einer Zeitung? Wenn möglich, möchte ich die beiden dazu bringen, dem zuzustimmen, damit ich das Diagramm mit geom_smooth und die Zusammenfassungstabelle von lm() in der Zeitung verwenden kann. Vielen Dank!

+0

Können Sie die Ausgabe aus 'Zusammenfassung' für die Daten (und das Modell) hinzufügen, so können wir sicher sein, dass wir die Daten richtig gelesen haben. – Spacedman

+0

@Spacedman, Ihre Ausgabe bestätigt, dass Sie die Daten richtig gelesen haben. Ich werde die Frage nicht mit der Ausgabe aktualisieren, da Sie die folgende Frage gelöst haben. – BonnieM

Antwort

3

Es sieht alles gut für mich. Die Zusammenfassung Linien für CAT5 sind:

Coefficients: 
      Estimate Std. Error t value Pr(>|t|)  
(Intercept) 1.932248 0.053131 36.368 < 2e-16 *** 
x   -0.006651 0.001962 -3.389 0.000721 *** 
... 
cat5  -1.080554 0.075138 -14.381 < 2e-16 *** 
... 
x:cat5  0.005602 0.002775 2.019 0.043720 * 

der die Steigung für CAT5 bedeutet die Gesamtneigung für x plus die Steigung für die x: CAT5-Interaktion:

> -0.006651+0.005602 
[1] -0.001049 

und auf dem Grundstück Ich sehe - 0,00105

der Schnitt als 0.852, die

ist
> 1.932248+(-1.080554) 
[1] 0.851694 

soweit ich, die beiden sehen also gezeigt die Dinge stimmen überein.

+0

Oh mein Gott, @ Spacedman. (Selbststirnschlag!) Danke, dass du meinen Fehler begreifst. Ich habe vergessen, die x-Steigung zur Interaktions-Steigung hinzuzufügen! Danke vielmals!!! – BonnieM

Verwandte Themen