I want perform a split plot model on R like I run using SAS.
 PROC GLM DATA=A; BY EXP;
       CLASS REPETICAO FATOR1 FATOR2;
       MODEL &VAR1 = REPETICAO FATOR1 REPETICAO*FATOR1 FATOR2 FATOR1*FATOR2;
       TEST  H=FATOR1 E=REPETICAO*FATOR1;
How I can run this model using R aov function?
I tried some like this:
aov( resp ~ REPETICAO+FATOR1*FATOR2+Error(REPETICAO*FATOR1), data=y)
I need that R generate the correct F test for "FATOR1"
Somebody can help?
 
    