[R-br] Teste de Kruskal Wallis!

Fernando Colugnati fernando em ipti.org.br
Terça Março 13 11:50:07 BRT 2012


Certamente os dois procedimentos estão utilizando estatísticas de teste
diferentes para comparações múltiplas. Repare que no kruskalmc o nível
crítico para alpha=0.05 é fixo, 19.7 aproximadamente, e no kruskal este
nível crítico muda...aparentemente a estatística é outra, talvez utilizando
uma aproximação para a Normal ou t, já que fornece um IC simétrico para as
diferenças entre os postos. O ideal seria saber quais são as estatísticas,
mas no manual do Agricolae não encontrei nada.

Não sei se ajudei muito...

Abs

Em 13 de março de 2012 10:22, Thiago de Paula Protásio <
depaulaprotasio em yahoo.com.br> escreveu:

> Olá Fernando!
>
> Abaixo seguem os resultados.
> A linha de comando é relativamente simples para as funções dos dois
> pacotes.
>
>
> > require(agricolae)> kruskal(O, Clone, alpha = 0.05, group=F, main = "Clone")
> Study: Clone
> Kruskal-Wallis test's
> Ties or no Ties
>
> Value: 11.6483
> degrees of freedom: 6
> Pvalue chisq  : 0.07029077
>
> Clone,  means of the ranks
>
>         O replication
> G008 22.5           5
> G084 16.4           5
> G122  6.8           5
> U030 14.4           5
> U034 20.3           5
> U060 26.8           5
> U083 18.8           5
>
> Comparison between treatments mean of the ranks
>
>             Difference   pvalue sig         LCL      UCL
> G008 - G084        6.1 0.301092      -5.7600575 17.96006
> G008 - G122       15.7 0.011312   *   3.8399425 27.56006
> G008 - U030        8.1 0.172796      -3.7600575 19.96006
> G008 - U034        2.2 0.706832      -9.6600575 14.06006
> U060 - G008        4.3 0.463864      -7.5600575 16.16006
> G008 - U083        3.7 0.527988      -8.1600575 15.56006
> G084 - G122        9.6 0.108468      -2.2600575 21.46006
> G084 - U030        2.0 0.732354      -9.8600575 13.86006
> U034 - G084        3.9 0.506096      -7.9600575 15.76006
> U060 - G084       10.4 0.083256   .  -1.4600575 22.26006
> U083 - G084        2.4 0.681656      -9.4600575 14.26006
> U030 - G122        7.6 0.199964      -4.2600575 19.46006
> U034 - G122       13.5 0.027140   *   1.6399425 25.36006
> U060 - G122       20.0 0.001776  **   8.1399425 31.86006
> U083 - G122       12.0 0.047532   *   0.1399425 23.86006
> U034 - U030        5.9 0.316920      -5.9600575 17.76006
> U060 - U030       12.4 0.041064   *   0.5399425 24.26006
> U083 - U030        4.4 0.453640      -7.4600575 16.26006
> U060 - U034        6.5 0.271122      -5.3600575 18.36006
> U034 - U083        1.5 0.797474     -10.3600575 13.36006
> U060 - U083        8.0 0.177984      -3.8600575 19.86006
>
>
>
> > require(pgirmess)> kruskalmc(O~Clone, probs = 0.05, cont=NULL)
>
> Multiple comparison test after Kruskal-Wallis
> p.value: 0.05
> Comparisons
>           obs.dif critical.dif difference
> G008-G084     6.1     19.68897      FALSE
> G008-G122    15.7     19.68897      FALSE
> G008-U030     8.1     19.68897      FALSE
> G008-U034     2.2     19.68897      FALSE
> G008-U060     4.3     19.68897      FALSE
> G008-U083     3.7     19.68897      FALSE
> G084-G122     9.6     19.68897      FALSE
> G084-U030     2.0     19.68897      FALSE
> G084-U034     3.9     19.68897      FALSE
> G084-U060    10.4     19.68897      FALSE
> G084-U083     2.4     19.68897      FALSE
> G122-U030     7.6     19.68897      FALSE
> G122-U034    13.5     19.68897      FALSE
> G122-U060    20.0     19.68897       TRUE
> G122-U083    12.0     19.68897      FALSE
> U030-U034     5.9     19.68897      FALSE
> U030-U060    12.4     19.68897      FALSE
> U030-U083     4.4     19.68897      FALSE
> U034-U060     6.5     19.68897      FALSE
> U034-U083     1.5     19.68897      FALSE
> U060-U083     8.0     19.68897      FALSE
>
>
>
> Obrigado!
>
>
>
> --
> Thiago de Paula Protásio
> Acadêmico de Engenharia Florestal
> Universidade Federal de Lavras
> Departamento de Ciências Florestais
> Ciência e Tecnologia da Madeira
> (035) 9183 - 2246
>
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>



-- 
Fernando A.B. Colugnati
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120313/c117727a/attachment.html>


Mais detalhes sobre a lista de discussão R-br