library(magrittr)
library(data.table)
library(lmerTest)
d <-
"x1 x2 x3 x4 y1 y2 y3
-1 -1 -1 -1 28 27.1 26.2
-1 -1 1 1 46.3 43.5 46.5
-1 1 -1 1 21.9 21 20.1
-1 1 1 -1 52.9 53.7 52
1 -1 -1 1 75 73.1 74.3
1 -1 1 -1 127.7 126.9 128.7
1 1 -1 -1 86.2 86.5 87
1 1 1 1 195 195.9 195.7"
) %>%
.[, subject := as.character(seq(1, nrow(d)))]
d.melt <-
melt(d,
id.vars = c("x1", "x2", "x3", "x4", "subject"),
value.name = "y")
m <- lmer(y ~ x1 + x2 + x3 + x4 + (1 | subject), data = d.melt)
summary(m) %>% print
fixef(m) %>% print
ranef(m) %>% print
resid(m) %>% print
bGlicmFyeShtYWdyaXR0cikKbGlicmFyeShkYXRhLnRhYmxlKQpsaWJyYXJ5KGxtZXJUZXN0KQoKZCA8LQogIGZyZWFkKAogICAgIngxICAgICAgeDIgICAgICB4MyAgICAgIHg0ICAgICAgeTEgICAgICB5MiAgICAgIHkzCiAgICAtMSAgICAgIC0xICAgICAgLTEgICAgICAtMSAgICAgIDI4ICAgICAgMjcuMSAgICAyNi4yCiAgICAtMSAgICAgIC0xICAgICAgMSAgICAgICAxICAgICAgIDQ2LjMgICAgNDMuNSAgICA0Ni41CiAgICAtMSAgICAgIDEgICAgICAgLTEgICAgICAxICAgICAgIDIxLjkgICAgMjEgICAgICAyMC4xCiAgICAtMSAgICAgIDEgICAgICAgMSAgICAgICAtMSAgICAgIDUyLjkgICAgNTMuNyAgICA1MgogICAgMSAgICAgICAtMSAgICAgIC0xICAgICAgMSAgICAgICA3NSAgICAgIDczLjEgICAgNzQuMwogICAgMSAgICAgICAtMSAgICAgIDEgICAgICAgLTEgICAgICAxMjcuNyAgIDEyNi45ICAgMTI4LjcKICAgIDEgICAgICAgMSAgICAgICAtMSAgICAgIC0xICAgICAgODYuMiAgICA4Ni41ICAgIDg3CiAgICAxICAgICAgIDEgICAgICAgMSAgICAgICAxICAgICAgIDE5NSAgICAgMTk1LjkgICAxOTUuNyIKICApICU+JQogIC5bLCBzdWJqZWN0IDo9IGFzLmNoYXJhY3RlcihzZXEoMSwgbnJvdyhkKSkpXQpkLm1lbHQgPC0KICBtZWx0KGQsCiAgICAgICBpZC52YXJzID0gYygieDEiLCAieDIiLCAieDMiLCAieDQiLCAic3ViamVjdCIpLAogICAgICAgdmFsdWUubmFtZSA9ICJ5IikKCm0gPC0gbG1lcih5IH4geDEgKyB4MiArIHgzICsgeDQgKyAoMSB8IHN1YmplY3QpLCBkYXRhID0gZC5tZWx0KQpzdW1tYXJ5KG0pICU+JSBwcmludApmaXhlZihtKSAlPiUgcHJpbnQKcmFuZWYobSkgJT4lIHByaW50CnJlc2lkKG0pICU+JSBwcmludAo=