[R-br] Ajuda na análise de seleção recorrente

Mauricio Couto couto_estatistica em hotmail.com
Segunda Junho 25 11:10:18 BRT 2012


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 em gmail.com
To: r-br em listas.c3sl.ufpr.br
Subject: Re: [R-br]	Ajuda na análise de seleção recorrente

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 em 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 em 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 em 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. 		 	   		  
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120625/e03626a2/attachment.html>


Mais detalhes sobre a lista de discussão R-br