Формула для многоуровневого lme () с повторными измерениями, 1 между субъектами, 1 внутри-субъектным фактором

Question

В RI хотелось бы проанализировать дизайн смешанных повторных измерений как многоуровневый glm с lme ().

Мои данные выглядят так:

• Две группы (между субъектами переменная = "Обработка")
• Каждый участник решал два вида задач (в пределах предмета переменная = "задача")

• Два раза измерения (повторное измерение = "время"), до обработки (= 1), после обработки (= 2)

• переменная результата RT

Таким образом, мой кадр данных выглядит примерно так:

dat

id     time   treatment    task    RT
1       1       a            x     0.2
1       1       a            y     0.3
1       2       a            x     0.1
1       2       a            y     0.2
2       1       b            x     0.2
2       1       b            y     0.3
2       2       b            x     0.01
2       2       b            y     0.02
...

Я хотел бы запустить многоуровневый glm, однако мне было интересно узнать о формуле.

Это

rtmodel<- lme( RT~ time + treatment + task + time:treatment + time:task + task:treatment + time:treatment:task, random= ~1|id/time/treatment/task, data= dat, method= "ML")

или это (потому что лечение это коэффициент между подсчетами):

rtmodel<- lme( RT~ time + treatment + task + time:treatment + time:task + task:treatment + time:treatment:task, random= ~1|id/time/task, data= dat, method= "ML")

?

Или я совершенно не прав насчет моей формулы?