Maurício, Vendo o que você postou... e consultando meu oráculo, posso suspeitar que o "erro" na sua ANOVA seja que ele está interpretando suas variáveis como numérica e não categóricas! Note que isso é um erro de semântica e não sintaxe! Para ver se é esse o caso, reporte-nos a saída do seus
str(couto)
A afirmativa que esses efeitos estejam sendo interpretados como variáveis numéricas, proceda:
couto <- transform(couto, AMB = factor(AMB), ... outras variáveis a transformar ... )
Então proceda de novo sua anova pelo aov(), fica a sugestão de usar o terms('modelo', keep = TRUE) no argumento formula da função... se for o caso de obter somas de quadrados ajustadas! OBS: Em outra ocasião que postar, leia antes o guia de postagem e produza um CMR, isso facilita para quem se disponibiliza a ajudar! sem mais... à disposição (parafrazeando Walmes), FH Faltou você transformar suas variáveis em fator 2012/6/20 Mauricio Couto <couto_estatistica@hotmail.com>
Pessoal.
gostaria da ajuda de vocês para analisar um modelo de seleção recorrente no R, o problema são os efeitos aninhados com interação (R/AS). Necessito estimar os componentes de variância para estimar por exemplo a herdabilidade entre outros.
Modelo estatístico:
Yijkl = µ + Ai + Sj + ASij+R/ASijk +F/Sjl + AS/Sijl + ε ijkl******
μ é a média****
Ai é o efeito do ambiente****
Sj é o efeito de set****
ASij é o efeito da interação ambiente e set****
R/ASijk é o efeito da repetição dentro de ambiente e set****
F/Sjl é o efeito da família dentro de set****
AF/Sijl é o efeito da interação ambiente e família dentro de set****
ɛijkl é o erro experimental
Obs? Vejam como eu tentei fazer:
couto<-read.table("C:/Users/Mauricio/Desktop/SR.txt", h=T)
couto<-aov(CE~AMB*SET+REP/AMB*SET+Famílias/SET+AMB*Famílias/SET,data=couto)
****
summary(couto)**** Obs: só que a ANOVA saiu errada.
Desde de já agradeço.
*Mauricio Farias Couto*
*Doutorando em Genética e Melhoramento de Plantas (UENF-RJ)*
*Mestre em Estatística Aplicada e Biometria (UFV-MG)*
*Engº Agronomo (UESC-BA)*
**
*Cel: (22) 8117-3725*
* *
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.