library(magrittr)
library(data.table)
library(lmerTest)

d <-
  fread(
    "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