В моем data.frame (dat) у меня есть следующие переменные:
st.id: идентификатор ученика
sc.id: идентификатор 4 школ, в которых учащиеся в (W/X/Y/Z)
group: определяет, что две школы - это лечение, а две другие - контроль (C/T)
gender: учащихся (M/F)
time: сколько раз тестировался каждый студент (1/ 2/ 3)
y: переменная непрерывного результата
с использованием lme4 или другого соответствующего пакета R, I Интересно, как написать модель для измерения эффекта group в этом проекте, учитывая вложенность учеников (st.id) в time и в школах (sc.id)?
Вот мой текущий код модели, но я не уверен, как включить time в случайную часть ниже?
lme4::lmer(y ~ group * gender * time + (time |sc.id/st.id), data = dat)