Toledo,
segui suas sugestões e em parte deu certo, contudo as intereções dentro de outro fator, tipo AMB*Familias/SET ou REP/AMB*SET, tá errado como vc pode ver abaixo:
couto<-read.table("C:/Users/Mauricio/Desktop/SR.txt",h=T)
> couto<-transform(couto,AMB=factor(AMB),SET=factor(SET),Famílias=factor(Famílias),REP=factor(REP))
> couto<-aov(CE~AMB*SET+REP/AMB*SET+Famílias/SET+ AMB*Famílias/SET,data=couto)
> summary(couto)
Df Sum Sq Mean Sq F value Pr(>F)
AMB 1 10 10.46 1.287 0.25726
SET 7 403 57.52 7.077 5.85e-08 ***
REP 1 57 56.56 6.959 0.00868 **
Famílias 192 13551 70.58 8.684 < 2e-16 ***
AMB:SET 7 98 13.97 1.719 0.10298
AMB:REP 1 8 7.84 0.965 0.32661
SET:REP 7 53 7.58 0.933 0.48060
AMB:Famílias 192 2012 10.48 1.289 0.01911 *
AMB:SET:REP 7 43 6.11 0.751 0.62845
Residuals 384 3121 8.13
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Ps: desculpe o incomodo.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
Date: Wed, 20 Jun 2012 15:55:10 -0300
From: fernandohtoledo@gmail.com
To: r-br@listas.c3sl.ufpr.br
Subject: Re: [R-br] Ajuda na análise de seleção recorrente_______________________________________________ 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.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),FHFaltou você transformar suas variáveis em fator2012/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.
_______________________________________________
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.