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!
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
@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