Random Intercept

Creating Random Intercept Model

Jika sebelumnya kita telah membuat model nul tanpa prediktor, dalam model intersep acak kita akan memasukkan prediktor. Misanya kita ingin mengetahui hubungan antara ESCS terhadap capaian matematika, kita dapat dengan mudah membuat model tersebut dengan langkah yang sama seperti ketika kita membuat model nul.

m1 <- lmer(MATH ~ 1 + ESCS + (1 | CNTSCHID), data = pisa)

The result

summary(m1)

ICC result

tab_model(m1)

Prediction plot

pisa$m1 <- predict(m1)

pisa %>% 
  ggplot(aes(ESCS, m1, color = CNTSCHID, group = CNTSCHID)) + 
  geom_smooth(se = F, method = lm) +
  theme_bw() +
  labs(x = "ESCS", 
       y = "Matematika", 
       color = "CNTSCHID")

QQ-plot

qqmath(ranef(m1, condVar = TRUE))