<div dir="ltr"><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Minha intervenção vai contra o contexto das mensagens de erro. Você tá aplicando um SK para nÃveis do fator tempo que contÃnuo. Embora não esteja errado do ponto de vista estatÃstico, do ponto de vista pÅ•atico/cientÃfico expressar o efeito do tempo por alguma função (polinômio ou modelo não linear) é considerado mais adequado por ser um fator qualitativo, no caso com 7 nÃveis. Outro ponto que não estou seguro é se o número de nÃveis deve ser o mesmo quando tem um fator só em DIC, como afirma o Ivan, mas no meu teste parece que é isso mesmo. Estranho, pois uma vez me atrevi a programar o teste e não me recordo se era requisito ter o mesmo número de observações por nÃvel.<br>
<br><span style="font-family:courier new,monospace">> require(ScottKnott)<br>> <br>> da <- data.frame(trat=gl(4,5))<br>> da$y <- rnorm(nrow(da))<br>> table(da$trat)<br><br>1 2 3 4 <br>5 5 5 5 <br>> <br>
> summary(SK(aov(y~trat, da), which="trat"))<br> Levels     Means SK(5%)<br>     3 0.2190466     a<br>     1 0.1650521     a<br>     2 -0.1431724     a<br>     4 -0.2933151     a<br>> <br>> db <- data.frame(trat=gl(2,3))<br>
> db$y <- rnorm(nrow(db))<br>> dc <- rbind(da, db)<br>> table(dc$trat)<br><br>1 2 3 4 <br>8 8 5 5 <br>> <br>> summary(SK(aov(y~trat, dc), which="trat"))<br> Levels     Means SK(5%)<br>     3 0.2190466     a<br>
     1 0.0773068     a<br>     4 -0.2933151     a<br>     2 -0.3164204     a<br>Mensagens de aviso perdidas:<br>1: In if ((lam < valchisq) | (ord1 == k)) { :<br> a condição tem comprimento > 1 e somente o primeiro elemento será usado<br>
2: In if (lam > valchisq) { :<br> a condição tem comprimento > 1 e somente o primeiro elemento será usado<br>3: In if (lam < valchisq) { :<br> a condição tem comprimento > 1 e somente o primeiro elemento será usado<br>
> <br></span><br>À disposição.<br>Walmes.<br></div></div>