<div dir="ltr">Lucas,<div><br></div><div>Conforme informei, o ks.test não resolveu.</div><div>Luiz Roberto</div></div><div class="gmail_extra"><br clear="all"><div><div>Luiz Roberto Martins Pinto<br>Prof. Pleno/DCET/UESC</div>
<div>Laboratório de Estatística Computacional</div><div>Universidade Estadual de Santa Cruz</div><div>Ilhéus-Bahia<br><br><a href="mailto:luizroberto.uesc@gmail.com" target="_blank">luizroberto.uesc@gmail.com</a><br>skype: lrmpinto</div>
<div><a href="http://lattes.cnpq.br/2732314327604831" target="_blank">http://lattes.cnpq.br/2732314327604831</a> <br><br><br></div></div>
<br><br><div class="gmail_quote">Em 23 de julho de 2014 16:31, <span dir="ltr"><<a href="mailto:luc.souza@usp.br" target="_blank">luc.souza@usp.br</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div><div style="font-family:times new roman,new york,times,serif;font-size:12pt;color:#000000">O teste de Kolmogorov-Smirnov (K-S) não resolve?<br><br><br><div><span name="x"></span>Lucas F. de Souza<br>Eng. Agrícola - CONAB<br>
M.Sc. Produção Vegetal - UENF<br>Doutor em Agrometeorologia - ESALQ/USP<br>tel: (19) 98160-7068<span name="x"></span><br></div><br><hr><blockquote style="border-left:2px solid rgb(16,16,255);margin-left:5px;padding-left:5px;color:#000;font-weight:normal;font-style:normal;text-decoration:none;font-family:Helvetica,Arial,sans-serif;font-size:12pt">
<b>De: </b>"Luiz Roberto Martins Pinto" <<a href="mailto:luizroberto.uesc@gmail.com" target="_blank">luizroberto.uesc@gmail.com</a>><br><b>Para: </b><a href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">r-br@listas.c3sl.ufpr.br</a><br>
<b>Enviadas: </b>Quarta-feira, 23 de Julho de 2014 10:40:53<br><b>Assunto: </b>[R-br] Teste para comparar distribuições<div><div class="h5"><br><br><div dir="ltr"><div>Caros amigos,</div><div><br></div><div>Fiz uma pesquisa no Google, mas não consegui uma solução satisfatória. Então peço ajuda.</div>
<div><br></div><div>Preciso de encontrar um teste que <b>quantifique</b> a semelhança entre dois vetores para checar se eles têm o mesmo tipo de distribuição, por meio de p.value.</div>
<div>O interesse principal é comparar os vetores x e ab.</div><div><br></div><div>Data set:<br></div><div><a href="http://www.datafilehost.com/d/f3757310" target="_blank">http://www.datafilehost.com/d/f3757310</a><br></div>
<div><br></div>
<div>library(Hmisc)</div><div>library(grDevices)</div><div>load('NNC D_d3kk2.RData')</div><div><br></div><div># Avaliação qualitativa (gráfica):</div><div><br></div><div>dev.set(which=1) </div><div>bpplot(d3)</div>
<div>boxplot(d3)</div><div><br></div><div><div># Aparentemente os vetores x, abC e ab são muitos semelhantes</div><div># A minha expectativa é encontrar um teste que quantifique esta semelhança por meio de p.value</div><div>
# Então fiz os testes abaixo</div></div><div><br></div><div># Avaliação quantitativa:</div><div><br></div><div>v2=d3$ab</div><div><br></div><div># Comparando o vetor x com o vetor ab:</div><div><br></div><div>v1=d3$x</div>
<div>round(ks.test(v1, v2)$p.value,2)</div><div><br></div><div># Out put: p.value=0. Então o vetor x é estatisticamente DIFERENTE do vetor ab.</div><div><br></div><div># Comparando o vetor abC com o vetor ab:</div><div><br>
</div><div>v1=d3$abC</div><div>round(ks.test(v1, v2)$p.value,2)<br></div><div><br></div><div># Out put: p.value=0.65. Então o vetor x é estatisticamente IGUAL do vetor abC.<br></div><div><br></div><div># Comparando o vetor abC com o vetor x:</div>
<div><br></div><div>v2=d3$x<br></div><div>round(ks.test(v1, v2)$p.value,2)<br></div><div><br></div><div># Out put: p.value=0. Então o vetor abC é estatisticamente DIFERENTE do vetor x.<br></div><div><br></div><div>### <b>Testes alternativos</b>, para comparar o vetor x com ab:</div>
<div><br></div><div><div>v2=d3$ab</div><div>v1=d3$x</div><div><br></div><div>t1=round(ks.test(v1, v2)$p.value,2);t1 # Out put: p.value=0</div><div>t2=round(ks.test(v1, v2, alternative = "l")$p.value,2);t2 # Out put: p.value=0<br>
</div><div>t3=round(ks.test(v1, v2, alternative = "g")$p.value,2);t3 # Out put: p.value=0.85</div></div><div><br></div><div># Então, para este Data set (NNC D_d3kk2.RData), e para o t3=0.85. Então, por meio deste teste, o vetor x é estatisticamente IGUAL do vetor ab.</div>
<div><br></div><div># <b>No entanto</b>, fazendo os testes com conjuntos de dados diferentes (10000 Data set) (mas, sempre com mesmo tamanho de vetor), quando comparo o vetor x com o vetor ab:</div><div><br></div><div> ## em aproximadamente 60% das comparações (ou em 6000 Data set) t1>0.05</div>
<div> ## em aproximadamente 80% das comparações (ou em 8000 Data set) t2>0.05</div><div> ## em aproximadamente 80% das comparações (ou em 8000 Data set) t3>0.05</div><div><br></div><div> ## em aproximadamente 85% das comparações (ou em 8500 Data set) t1 <b><font color="#ff0000">OU</font></b> t2 >0.05 </div>
<div> ## em aproximadamente 85% das comparações (ou em 8500 Data set) t1 <b><font color="#ff0000">OU</font></b> t3 >0.05 </div><div> </div><div> ## em 0% das comparações (ou em 0 Data set) t2 <font color="#ff0000"><b>E</b></font> t3 >0.05 </div>
<div> ## em 100% das comparações (ou em 10000 Data set) t2 <b><font color="#ff0000"><span>OU</span> </font></b>t3 >0.05 </div><div><br></div><div><br></div><div><div>Luiz Roberto Martins Pinto<br>
Prof. Pleno/DCET/UESC</div><div>Laboratório de Estatística Computacional</div><div>Universidade Estadual de Santa Cruz</div><div>Ilhéus-Bahia<br><br><a href="mailto:luizroberto.uesc@gmail.com" target="_blank">luizroberto.uesc@gmail.com</a><br>
skype: lrmpinto</div>
<div><a href="http://lattes.cnpq.br/2732314327604831" target="_blank">http://lattes.cnpq.br/2732314327604831</a> <br><br><br></div></div>
</div>
<br></div></div>_______________________________________________<br>R-br mailing list<br><a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">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.</blockquote><br></div></div><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>