desdobramento de uma analise fatorial

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

Segue um CMR. data(ex5, package="ExpDes") str(ex5) xtabs(~trat+genero, ex5) m0 <- aov(sabor~bloco+trat/genero, data=ex5) names(coef(m0)) m0$assign gINt.effects <- names(coef(m0))[m0$assign==3] generoINtrat <- lapply(paste0("trat", levels(ex5$trat)), grep, gINt.effects) summary(m0, split=list("trat:genero"=generoINtrat)) m1 <- aov(sabor~bloco+genero/trat, data=ex5) names(coef(m1)) m1$assign tINg.effects <- names(coef(m1))[m1$assign==3] tratINgenero <- lapply(paste0("genero", levels(ex5$genero)), grep, tINg.effects) summary(m1, split=list("genero:trat"=tratINgenero)) À disposição. Walmes.
participantes (2)
-
helcio duarte pereira
-
walmes .