<div><span style="font-family: monospace; font-size: large;">Não estou conseguindo acessar este arquivo, não reconhece o comando gdata nem o read.xls, e coloquei a library readxl</span></div>
<div><span style="font-family: monospace; font-size: large;">da <- gdata::read.xls("db_r-list.xlsx",</span><br style="font-family: monospace; font-size: 10.4px;" /><span style="font-family: monospace; font-size: large;"> stringsAsFactors = FALSE)</span></div>
<div><span style="font-size: large;"> </span></div>
<div><br /><br /></div>
<hr style="border-top: 1px solid #ccc;" />
<div><br /><strong>De:</strong> r-br@listas.c3sl.ufpr.br<br /><strong>Enviada:</strong> Quarta-feira, 2 de Outubro de 2019 11:05<br /><strong>Para:</strong> henrique.bdias@yahoo.com.br<br /><strong>Assunto:</strong> [R-br] Ajuda com análise conjunta de dois experimentos com parcelas subdividas<br /><br /></div>
<div dir="ltr">
<div class="gmail_default" style="font-family: trebuchet ms,sans-serif;">
<div class="gmail_default" style="font-family: trebuchet ms,sans-serif;">Henrique,</div>
<div class="gmail_default" style="font-family: trebuchet ms,sans-serif;"> </div>
<div class="gmail_default" style="font-family: trebuchet ms,sans-serif;">Bom que você conferiu reprodutibilidade MAS o código de conduta da lista pede para que não sejam enviados anexos para a lista.</div>
<div class="gmail_default" style="font-family: trebuchet ms,sans-serif;">Existem outras formas de conferir acesso aos dados sem que seja por anexo, ok?</div>
<div class="gmail_default" style="font-family: trebuchet ms,sans-serif;"> </div>
<div class="gmail_default" style="font-family: trebuchet ms,sans-serif;">> As dúvidas, além de saber se estou ou não num caminho lúcido em termos da estatística:<br />> 1) Em alguns casos houve violação de alguma pressuposição (estou testando norm. dos residuos e homogeneidade apenas),</div>
<div class="gmail_default" style="font-family: trebuchet ms,sans-serif;">> como para a variável resposta 3, no local II colhido na idade 8.<br />> - Procedo com alguma transformação ? tipo utilizar boxcox e ver qual melhor forma de transformar os dados ?<br />> - Ou Parto para outro tipo de abordagem ? alguma sugestão ?</div>
<div class="gmail_default" style="font-family: trebuchet ms,sans-serif;"> </div>
<div class="gmail_default" style="font-family: trebuchet ms,sans-serif;">Nem uma coisa nem outra. Segue com a análise usual. Por que? Bem, para começar, você sabe que, se o nível de significância é de 5%, 1 em cada 20 análises (em média) irá rejeitar a hipótese nula mesmo que ela seja verdadeira. É o erro tipo I. Então, o que está acontecendo com você é basicamente isso. De tanto fazer testes, uma hora haverá uma rejeição que pode ser apenas casual. Não rodei seu script para saber se de fato a rejeição deve-se a alguma anomalia grave. Mas de qualquer forma, fica aí o raciocínio. Outro ponto, eu não gosto dessa abordagem "pilha de testes de hipóteses". Prefiro, ainda que seja considerado subjetiva, a análise gráfica dos resíduos porque a leitura conjunta dos gráficos subsidia alguma ação em caso de necessidade, coisa que o teste não faz porque ele diz apenas "sim" e "não". Um último ponto, sendo a mesma variável resposta vista em divisões dos dados (por local, por época, etc), não seria desejável aplicar uma transformação diferente em cada fração. O melhor seria encontrar uma transformação que seja adequada para todas as frações. Você pode examinar isso graficamente, fazer uma boa análise exploratória e depois de ajustar o modelo, fazer uma boa análise dos resíduos.</div>
<div class="gmail_default" style="font-family: trebuchet ms,sans-serif;"> </div>
<div class="gmail_default" style="font-family: trebuchet ms,sans-serif;">A `VariavelResposta1`, por exemplo, mostra um aumento de variabilidade em função da idade, coisa que é comum em dados de crescimento (não se é o caso). No início as unidades experimentais de um mesmo tratamento são mais parecidas, mas com o passar do tempo na exposição aos várias fatores não controlados, as diferenças nas variáveis resposta tendem a aumentar, aumentando portanto a variância do erro. Talvez uma transformação log, raíz, etc, resolva isso. Mas nada impede de modelar a variância conjuntamente com a média (sem querer complicar).</div>
<div class="gmail_default" style="font-family: trebuchet ms,sans-serif;"> </div>
<div class="gmail_default" style="font-family: trebuchet ms,sans-serif;">> 2) Há alguma forma de análise considerando as todas as idades?</div>
<div class="gmail_default" style="font-family: trebuchet ms,sans-serif;"> </div>
<div class="gmail_default" style="font-family: trebuchet ms,sans-serif;">Se você considerar representar efeito do fator idade (fator quantitativo) por alguma função, por exemplo, polinômio, então mesmo não tendo as mesmas idades (mesmo sendo um fatorial incompleto) nas cultivares, você vai analisar tudo conjuntamente com facilidade. Se assumir idade como variável quantitativa (o que não é recomendável pela própria natureza da variável), aí você terá um fatorial incompleto que é coisa complicada de analisar por causa das celas ausentes.</div>
<div class="gmail_default" style="font-family: trebuchet ms,sans-serif;"> </div>
<div class="gmail_default" style="font-family: trebuchet ms,sans-serif;">> 3) Ao menos duas idades e 3 dos cultivares são comuns em ambos os locais, seria muito abusado tentar realizar uma análise conjunta</div>
<div class="gmail_default" style="font-family: trebuchet ms,sans-serif;">> apenas para essas cultivares? (com abordagem similar ao proposto para o caso do José)</div>
<div class="gmail_default" style="font-family: trebuchet ms,sans-serif;"> </div>
<div class="gmail_default" style="font-family: trebuchet ms,sans-serif;">Poderia ser feito. A idade não é preocupação se for usado regressão nela. Mas a análise exploratória que fiz indica que só existe uma época comum aos dois locais. Dessa forma não tem como declara interação local:epoca, pois pela ausência de mais combinações comuns, o efeito de interação não é estimável se época for considerado fator qualitativo. Se for quantitativo, assumindo uma função, então daria para estimar/representar a interação.</div>
<div class="gmail_default" style="font-family: trebuchet ms,sans-serif;"> </div>
<div class="gmail_default" style="font-family: trebuchet ms,sans-serif;">Tudo isso vai do seu interesse e disposição em fazer a análise. Ao meu ver, o experimento parece não ter sido planejado para uma análise conjunta visando estimar interação etc. O jeito mais tranquilo é analisar os locais separados, até mesmo porque, não se tem como inferir sobre o efeito de local (não creio que essa seja a hipótese), e uma vez que estimar a interação local:epoca é complicado neste caso, não sei se compensa o esforço. A interação local:cultivar dá para estimar, mas apenas baseando-se em 3 cultivares comuns aos locais, mas tendo em mente que as épocas de cultivo foram diferentes (apenas uma comum), o que certamente incluirá um complicador para gerar confundimento.</div>
<div class="gmail_default" style="font-family: trebuchet ms,sans-serif;"> </div>
<div class="gmail_default" style="font-family: trebuchet ms,sans-serif;">Já as idades, como em cada experimento elas são as mesmas, e o padrão é bem linear, você pode começar com o modelo linear simples pro efeito de idade e ir aumentando o grau conforme necessidade.</div>
<div> </div>
<div>
<div class="gmail_default" style="font-family: trebuchet ms,sans-serif;">Segue código usado para explorar os dados.</div>
<div class="gmail_default" style="font-family: trebuchet ms,sans-serif;"> </div>
<div class="gmail_default"><span style="font-family: monospace;"># Lê direto da planilha.<br />da <- gdata::read.xls("db_r-list.xlsx",<br /> stringsAsFactors = FALSE)<br />str(da)<br /><br />factor_roman <- function(x) {<br /> u <- sort(as.roman(unique(x)))<br /> factor(x, levels = as.character(u))<br />}<br /><br /># Converte para fator com níveis romanos ordenados.<br />da <- transform(da,<br /> Epoca = factor_roman(Epoca),<br /> Cultivar = factor_roman(Cultivar),<br /> Local = factor_roman(Local),<br /> Bloco = factor(Repeticao))<br /><br /># Tabelas de frequência dos pontos experimentais.<br />xtabs(~Epoca + Local, data = da)<br />xtabs(~Cultivar + Local, data = da)<br />xtabs(~Cultivar + Epoca, data = da)<br />ftable(xtabs(~Local + Idade + Cultivar, data = da))<br /><br />library(lattice)<br /><br /># Diagrama de dispersão.<br />xyplot(VariavelResposta1 ~ Idade | Epoca + Local,<br /> groups = Cultivar,<br /> data = da,<br /> type = c("p", "a"))<br /><br /># Separando por local.<br />gridExtra::grid.arrange(<br /> ncol = 1,<br /> xyplot(VariavelResposta1 ~ Idade | Epoca,<br /> groups = Cultivar,<br /> data = da,<br /> subset = Local == "I",<br /> layout = c(NA, 1),<br /> type = c("p", "a")),<br /> xyplot(VariavelResposta1 ~ Idade | Epoca,<br /> groups = Cultivar,<br /> data = da,<br /> subset = Local == "II",<br /> layout = c(NA, 1),<br /> type = c("p", "a")))<br /><br /># Cultivares comuns aos dois locais.<br />u <- as.array(xtabs(~Cultivar + Local, data = da))<br />i <- apply(u, 1, prod) > 0<br />common <- as.data.frame(u[i, ])<br />common$Freq <- NULL<br />common<br /><br /># Pega apenas as combinações presentes.<br />db <- droplevels(merge(common, da, all.x = TRUE, all.y = FALSE))<br />summary(db)<br /><br /># Diagrama de dispersão para a seleção.<br />xyplot(VariavelResposta1 ~ Idade | Epoca + Local,<br /> groups = Cultivar,<br /> data = db,<br /> type = c("p", "a"))<br /><br />xyplot(VariavelResposta1 ~ Idade | Cultivar + Local,<br /> groups = Epoca,<br /> data = db,<br /> type = c("p", "a"))</span></div>
<div class="gmail_default" style="font-family: trebuchet ms,sans-serif;"> </div>
<div class="gmail_default" style="font-family: trebuchet ms,sans-serif;">À disposição.</div>
<div class="gmail_default" style="font-family: trebuchet ms,sans-serif;">Walmes.</div>
</div>
</div>
</div>
<div>_______________________________________________<br />R-br mailing list<br />R-br@listas.c3sl.ufpr.br<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 /><br /></div>