Walmes,<div><br></div><div>Obrigado pela prontidão!</div><div><br></div><div>Pois bem, rodei o seu código... eis que: sua implementação confere com os resultados esperados relativos ao ajuste das médias com recuperação da informação interbloco (popmeans, do doBy quando o modelos é ajustado pelo lme()).</div>
<div><br></div><div>No primeiro ajuste,, o do aov, o popmeans retorna as próprias médias aritméticas dos tratamentos, sem qualquer ajuste!</div><div><br></div><div>O teste correto da soma de quadrados de tratamentos ajustados confere com Cochran & Cox!</div>
<div><br></div><div>Infelizmente ainda não consegui recuperar as médias ajustadas intrabloco, mas com certeza o que é posto por você ajuda e muito!</div><div><br></div><div>Obrigado...</div><div><br></div><div>abraço,</div>
<div>FH<br><br><div class="gmail_quote">2012/11/15 Walmes Zeviani <span dir="ltr"><<a href="mailto:walmeszeviani@gmail.com" target="_blank">walmeszeviani@gmail.com</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div><font face="trebuchet ms,sans-serif">Fernando,</font></div><div><font face="trebuchet ms,sans-serif"><br></font></div><div><font face="trebuchet ms, sans-serif">Já rodei os exemplos do Pimentel e do Ramalho et al. Implementação de ambos estão disponíveis nos scripts que levei para o Curso mais recente que dei na Embrapa Arroz e Feijão (que terminei ontem). O link para arquivos do curso é esse</font></div>
<div><font face="trebuchet ms, sans-serif"><br></font></div><div><a href="http://www.leg.ufpr.br/~walmes/cursoR/cnpaf2/" target="_blank">http://www.leg.ufpr.br/~walmes/cursoR/cnpaf2/</a></div><div><font face="trebuchet ms, sans-serif"><br>
</font></div><div><font face="trebuchet ms, sans-serif">Pegue o script </font><font face="courier new, monospace">lati.R</font><font face="trebuchet ms, sans-serif">. Eu nunca olhei a implementação do C&Cox e fiquei curioso e imaginando será que os livros divergem? Meus códigos reproduzem o Pimentel. Bem, com os seus dados a análise que atualmente faço é essa</font></div>
<div><font face="trebuchet ms, sans-serif"><br></font></div><div><div><font face="courier new, monospace">rm(list=ls(all=TRUE)); ls()</font></div><div class="im"><div><font face="courier new, monospace">download <- read.table('<a href="http://dl.dropbox.com/u/38195533/dados_CC.txt" target="_blank">http://dl.dropbox.com/u/38195533/dados_CC.txt</a>',</font></div>
<div><font face="courier new, monospace"> header = TRUE, # com cabecalho</font></div><div><font face="courier new, monospace"> sep = '\t', # separador de celulas <TAB></font></div>
<div><font face="courier new, monospace"> dec = ',', # separador de decimal <, (virgula)></font></div><div><font face="courier new, monospace"> na.string = '.') # indicador de omissao</font></div>
</div><div><font face="courier new, monospace">str(download)</font></div><div class="im"><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">## lendo e transformando variaveis em fator</font></div>
<div><font face="courier new, monospace">dados <- transform(download, </font></div><div><font face="courier new, monospace"> trat = factor(trat), # transforma em fator</font></div><div><font face="courier new, monospace"> rep = factor(rep), # idem</font></div>
<div><font face="courier new, monospace"> bloco = factor(bloco)) # ...</font></div><div><font face="courier new, monospace"><br></font></div></div><div><font face="courier new, monospace">all(complete.cases(dados)) # completo</font></div>
<div><font face="courier new, monospace">xtabs(~rep+trat, dados)</font></div><div><font face="courier new, monospace">xtabs(~rep+bloco, dados)</font></div><div><font face="courier new, monospace">xtabs(~bloco+trat+rep, dados)</font></div>
<div class="im">
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">str(dados) # estrutura da planilha</font></div><div><font face="courier new, monospace"><br></font></div></div><div><font face="courier new, monospace">m0 <- lm(terms(resp~rep/bloco+trat, keep.order=TRUE), data=dados)</font></div>
<div><font face="courier new, monospace">anova(m0)</font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">require(doBy)</font></div><div><font face="courier new, monospace">popMeans(m0, effect="trat")</font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">require(nlme)</font></div><div><font face="courier new, monospace">dados$blocoin <- with(dados, interaction(rep, bloco, drop=TRUE))</font></div>
<div><font face="courier new, monospace">mm0 <- lme(resp~rep+trat, random=~1|blocoin, dados)</font></div><div><font face="courier new, monospace">anova(mm0)</font></div><div><font face="courier new, monospace"><br></font></div>
<div><font face="courier new, monospace">popMeans(mm0, effect="trat")</font></div><div style="font-family:'trebuchet ms',sans-serif"><br></div></div><div style="font-family:'trebuchet ms',sans-serif">
À disposição.</div><div><font face="trebuchet ms, sans-serif">Walmes.</font></div><font face="trebuchet ms,sans-serif"><br clear="all"></font><span style="font-family:trebuchet ms,sans-serif">==========================================================================</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">Walmes Marques Zeviani</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">Departamento de Estatística - Universidade Federal do Paraná</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">fone: <a href="tel:%28%2B55%29%2041%203361%203573" value="+554133613573" target="_blank">(+55) 41 3361 3573</a></span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">VoIP: (3361 3600) 1053 1173</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">e-mail: <a href="mailto:walmes@ufpr.br" target="_blank">walmes@ufpr.br</a><br>
skype: walmeszeviani<br style="font-family:trebuchet ms,sans-serif"></span><span style="font-family:trebuchet ms,sans-serif">twitter: @walmeszeviani</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">homepage: <a href="http://www.leg.ufpr.br/%7Ewalmes" target="_blank">http://www.leg.ufpr.br/~walmes</a></span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">linux user number: 531218</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">==========================================================================</span><br>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br></div>