Deutschsprachiges Forum zur Statistikumgebung R
http://forum.r-statistik.de/
Aus deinem Model brauchst du nur die Predictions (auf grp level) in ein data.frame packen (geht z.b. mit broom sehr angenehm: https://cran.r-project.org/web/packages/broom.mixed/vignettes/broom_mixed_intro.html, oder halt selbst mit predict() rechnen). Danach kann man das mit ggplot2 rechnen...irgendwie gibt es diese Package nicht (mehr)? was gibt es für ALternativen wenn ich eine Hierarchische Lin Reg bunt darstellen will.
Code: Alles auswählen
library("lme4")
library("broom.mixed")
library("ggplot2")
mod <- lmer(Reaction ~ Days + (1 | Subject), sleepstudy)
pdat <- augment(mod)
head(pdat)
ggplot(pdat, aes(x = Days, y = .fitted, col = Subject)) +
geom_line() +
geom_line(aes(y = .fixed), size = 2, col = "black")