<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style></head>
<body class='hmmessage'><div dir='ltr'>
Pessoal.<BR>
<BR>
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.<BR>
<BR>
Modelo estatístico:<BR>
<P style="MARGIN: 0cm 0cm 10pt" class=MsoNormal><SPAN style="FONT-FAMILY: Arial"><FONT size=3>Y<SUB>ijkl</SUB> = µ + A<SUB>i</SUB> + S<SUB>j</SUB> + AS<SUB>ij</SUB>+R/AS<SUB>ijk</SUB> +F/S<SUB>jl</SUB> + AS/S<SUB>ijl</SUB> + ε <SUB>ijkl</SUB><?xml:namespace prefix = o ns = "urn:schemas-microsoft-com:office:office" /><o:p></o:p></FONT></SPAN></P>
<P style="MARGIN: 0cm 0cm 10pt" class=MsoNormal><SPAN style="LINE-HEIGHT: 115%; FONT-FAMILY: 'Times New Roman'; FONT-SIZE: 12pt">μ é a média<o:p></o:p></SPAN></P>
<P style="MARGIN: 0cm 0cm 10pt" class=MsoNormal><SPAN style="LINE-HEIGHT: 115%; FONT-FAMILY: 'Times New Roman'; FONT-SIZE: 12pt">A<SUB>i</SUB> é o efeito do ambiente<o:p></o:p></SPAN></P>
<P style="MARGIN: 0cm 0cm 10pt" class=MsoNormal><SPAN style="LINE-HEIGHT: 115%; FONT-FAMILY: 'Times New Roman'; FONT-SIZE: 12pt">S<SUB>j</SUB> é o efeito de set<o:p></o:p></SPAN></P>
<P style="MARGIN: 0cm 0cm 10pt" class=MsoNormal><SPAN style="LINE-HEIGHT: 115%; FONT-FAMILY: 'Times New Roman'; FONT-SIZE: 12pt">AS<SUB>ij</SUB> é o efeito da interação ambiente e set<o:p></o:p></SPAN></P>
<P style="MARGIN: 0cm 0cm 10pt" class=MsoNormal><SPAN style="LINE-HEIGHT: 115%; FONT-FAMILY: 'Times New Roman'; FONT-SIZE: 12pt">R/AS<SUB>ijk</SUB> é o efeito da repetição dentro de ambiente e set<o:p></o:p></SPAN></P>
<P style="MARGIN: 0cm 0cm 10pt" class=MsoNormal><SPAN style="LINE-HEIGHT: 115%; FONT-FAMILY: 'Times New Roman'; FONT-SIZE: 12pt">F/S<SUB>jl</SUB> é o efeito da família dentro de set<o:p></o:p></SPAN></P>
<P style="MARGIN: 0cm 0cm 10pt" class=MsoNormal><SPAN style="LINE-HEIGHT: 115%; FONT-FAMILY: 'Times New Roman'; FONT-SIZE: 12pt">AF/S<SUB>ijl</SUB> é o efeito da interação ambiente e família dentro de set<o:p></o:p></SPAN></P>
<P style="MARGIN: 0cm 0cm 10pt" class=MsoNormal><SPAN style="LINE-HEIGHT: 115%; FONT-FAMILY: 'MS Mincho'; FONT-SIZE: 12pt; mso-bidi-font-family: 'MS Mincho'">ɛ</SPAN><SUB><SPAN style="LINE-HEIGHT: 115%; FONT-FAMILY: 'Times New Roman'; FONT-SIZE: 12pt">ijkl</SPAN></SUB><SPAN style="LINE-HEIGHT: 115%; FONT-FAMILY: 'Times New Roman'; FONT-SIZE: 12pt"> é o erro experimental</SPAN></P>
<P style="MARGIN: 0cm 0cm 10pt" class=MsoNormal><SPAN style="LINE-HEIGHT: 115%; FONT-FAMILY: 'Times New Roman'; FONT-SIZE: 12pt"></SPAN> </P>
<P style="MARGIN: 0cm 0cm 10pt" class=MsoNormal><SPAN style="LINE-HEIGHT: 115%; FONT-FAMILY: 'Times New Roman'; FONT-SIZE: 12pt">Obs? Vejam como eu tentei fazer:</SPAN></P><SPAN style="LINE-HEIGHT: 115%; FONT-FAMILY: 'Times New Roman'; FONT-SIZE: 12pt">
<P style="MARGIN: 0cm 0cm 10pt" class=MsoNormal><SPAN style="FONT-FAMILY: 'Times New Roman'; mso-ansi-language: EN-US" lang=EN-US>> couto<-read.table("C:/Users/Mauricio/Desktop/SR.txt", h=T)</SPAN></P><SPAN style="FONT-FAMILY: 'Times New Roman'; mso-ansi-language: EN-US" lang=EN-US>
<P style="MARGIN: 0cm 0cm 10pt" class=MsoNormal><SPAN style="FONT-FAMILY: 'Times New Roman'; mso-ansi-language: EN-US" lang=EN-US>>couto<-aov(CE~AMB*SET+REP/AMB*SET+Famílias/SET+AMB*Famílias/SET,data=couto)<o:p></o:p></SPAN></P><o:p></o:p></SPAN></SPAN>
<P style="MARGIN: 0cm 0cm 10pt" class=MsoNormal><SPAN style="FONT-FAMILY: 'Times New Roman'; mso-ansi-language: EN-US" lang=EN-US><FONT size=3>> summary(couto)<o:p></o:p></FONT></SPAN></P>
Obs: só que a ANOVA saiu errada.<BR>
<BR>
Desde de já agradeço.<BR>
<BR>
<BR>
<BR><BR> <BR>
<P align=center><FONT color=#008000><STRONG>Mauricio Farias Couto</STRONG></FONT></P>
<P align=center><STRONG><FONT color=#008000>Doutorando em Genética e Melhoramento de Plantas (UENF-RJ)</FONT></STRONG></P>
<P align=center><STRONG><FONT color=#008000>Mestre em Estatística Aplicada e Biometria (UFV-MG)</FONT></STRONG></P>
<P align=center><FONT color=#008000><STRONG>Engº Agronomo (UESC-BA)</STRONG></FONT></P>
<P align=center><FONT color=#008000><STRONG></STRONG></FONT> </P>
<P align=center><FONT color=#008000><STRONG>Cel: (22) 8117-3725</STRONG></FONT></P>
<P align=center><STRONG><FONT color=#008000> </FONT></STRONG></P>
<BR> </div></body>
</html>