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?