<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/">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><font face="courier new, monospace">download <- read.table('<a href="http://dl.dropbox.com/u/38195533/dados_CC.txt">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><font face="courier new, monospace">str(download)</font></div><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><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><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><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: (+55) 41 3361 3573</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>