[R-br] desdobramento de uma analise fatorial

[helcio duarte pereira]

Pessoal, estou tentando desdobrar uma interação pela função aov mas só dá erro. Estou trabalhando com o objeto "ex5" da função "ExpDes", cuja estrutura basica é:

str(ex5)
'data.frame':160 obs. of  4 variables:
 $ trat  : Factor w/ 4 levels "10g","15g","15t",..: 1 1 1 1 1 
 $ genero: Factor w/ 2 levels "F","M": 1 1 1 1 1 1 1 1 1 1 ...
 $ bloco : Factor w/ 20 levels "1","2","3","4",..: 1 2 3 4 5 6 7
 $ sabor : int  8 9 7 9 8 8 8 5 6 7 ...

o modelo de analise que usei foi:
anava<-aov(formula = sabor~trat*genero + bloco, data = ex5)

Response: sabor
             Df Sum Sq Mean Sq F value  Pr(>F)  
trat          3  19.37  6.4562  1.9452 0.12537  
genero        1   7.66  7.6563  2.3068 0.13118  
bloco        19  97.82  5.1484  1.5512 0.07832 .
trat:genero   3   3.92  1.3062  0.3936 0.75783  
Residuals   133 441.43  3.3190    

Eu sei que a interação foi não significativa, mas estou tentando desdobra-la para aprender, já que irei precisar fazer uma analise assim em breve e que provavelmente vai dar interação significativa

O primeiro desdobramento(genero dentro de tratamento) rodou sem problemas:

desd1<-aov(formula = ex5$sabor~ex5$trat/ex5$genero + bloco, 
           data = ex5)
summary.aov(object = desd1, split = list("ex5$trat:ex5$genero"=
                  list("10g"=1, "15g"=2, "15t"=3, "20t"=4)))

Porém o segundo:

desd2<-aov(formula = ex5$sabor~ex5$genero/ex5$trat + bloco, 
           data = ex5)
summary.aov(object = desd2, split = list("ex5$genero:ex5$trat"=
                    list("F"=1, "M"=2)))

Não roda de maneira alguma. Usei o mesmo procedimento do anterior, assim, alguém sabe no que estou errando?

Att.
Helcio
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140927/68028ee0/attachment.html>