Eu gostaria de executar anova de medida repetida em R usando modelos de regressão, em vez de uma 'Análise de Variância' (AOV) função.

Aqui está um exemplo do meu código AOV para 3 fatores dentro do assunto:

m.aov<-aov(measure~(task*region*actiontype) + Error(subject/(task*region*actiontype)),data)  

Alguém pode me dar a sintaxe exata para executar a mesma análise usando modelos de regressão? Eu quero ter certeza de respeitar a independência dos resíduos, ou seja, usar termos específicos de erro como o AOV.

Em um post anterior eu li uma resposta do tipo:

lmer(DV ~ 1 + IV1*IV2*IV3 + (IV1*IV2*IV3|Subject), dataset))

Eu realmente não tenho certeza sobre essa solução, uma vez que ela ainda trata as variáveis ​​entre os assuntos, e eu não entendo como adicionar fatores aleatórios mudaria isso.

Alguém sabe executar anova de medida repetida com lm ​​/ lmer levando em conta a independência residual?

Muito obrigado, Solene