<div dir="ltr"><div class="gmail_default" style="font-family:georgia,serif;font-size:small">Oi, Benilton,<br>Obrigada.<br>Já venho estudando sobre GLMs e continuarei buscando informações sobre essas características dos meus dados para conseguir montar um modelo apropriado.<br>
Valeu pelas informações!<br>Polliana.</div><div class="gmail_extra"><br><br><div class="gmail_quote">Em 30 de maio de 2013 14:56, Benilton Carvalho <span dir="ltr"><<a href="mailto:beniltoncarvalho@gmail.com" target="_blank">beniltoncarvalho@gmail.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Para que alguem possa reproduzir o que vc quer dizer (se vc nao se<br>
importar com confidencialidade dos dados), eu sugiro que vc poste o<br>
resultado de:<br>
<br>
dput(alelo)<br>
<br>
Se nao quiser colocar seus dados aqui, crie um conjunto de dados de<br>
exemplo (q "finja" ser os seus dados originais) e poste o resultado de<br>
dput, conforme acima.<br>
<br>
Definitivamente (pelo menos da forma q vc esta' utilizando),<br>
Kolmogorov-Smirnov nao e' adequado... Se vc ler o manual do comando,<br>
vera' que ele so' aceita 2 grupos por vez. E, apesar de vc nao dizer<br>
explicitamente, eu imagino que vc queira testar as diferentes colunas<br>
do data.frame alelo (control, ph5.7, v100, v200, etc).<br>
<br>
Certamente, havera' pessoas interessadas em te dar sugestoes. Apesar<br>
de eu nao recomendar um metodo em particular, acho valido vc olhar as<br>
tecnicas que descrevi no meu email anterior. Para usa-las, voce<br>
precisara' "melt" os dados (o exemplo abaixo e' o mais simples, para<br>
casos nos quais vc nao tenha a variavel 'tempo', entao sera'<br>
necessario estudar o funcionamento do "melt" para faze-lo funcionar<br>
direitinho)...<br>
<br>
library(reshape2)<br>
novosDados = melt(alelo)<br>
<br>
<br>
E, entao, usar 'novosDados' para a modelagem... Voce basicamente faria<br>
um modelo parecido com:<br>
<br>
germinados ~ hora + tratamento<br>
<br>
ou ate'<br>
<br>
germinados ~ hora * tratamento<br>
<br>
Se vc estiver fazendo isso sozinha, e' definitivamente importante que<br>
voce estude a teoria destes modelos (ao inves de simplesmente usar o<br>
que postarem aqui, pois cada caso e' um caso). Se vc estiver fazendo<br>
isso com o acompanhamento de um estatistico (recomendado!), consulte-o<br>
a respeito das alternativas recomendadas por ele e suas implementacoes<br>
em R. Mas os modelos que mencionei anteriormente (Poisson, Binomial<br>
Negativa... e, possivelmente ambos com inflacao no zero) sao<br>
estrategias comumente consideradas em casos similares.<br>
<br>
b<br>
<br>
b<br>
<br>
<br>
<br>
<br>
Em 30 de maio de 2013 14:36, Polliana Zocche <<a href="mailto:farbby@gmail.com">farbby@gmail.com</a>> escreveu:<br>
<div class="HOEnZb"><div class="h5">> Olá, Benilton,<br>
> Desculpe a falta de clareza. Vou tentar explicar com mais informações.<br>
> Observei o efeito de 8 extratos aquosos de uma planta sobre a germinação de<br>
> sementes de alface. Então a coleta dos dados foi a contagem do total de<br>
> sementes germinadas em cada tratamento por dia.<br>
> _______<br>
> Exemplo do formato dos dados<br>
>> alelo$time #variável X (tempo em horas)<br>
> [1] 24 48 72 96 120 144 168 192 216 240 264 288 312 336 360<br>
> alelo$control #exemplo dos dados<br>
> [1] 91 95 96 98 98 99 99 99 99 99 100 100 100 100 100<br>
>> alelo$v200<br>
> [1] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0<br>
> alelo$s100<br>
> [1] 38 45 45 45 45 46 46 45 46 47 48 49 49 51 55 56 56 56 56 58<br>
>> alelo$v200ph<br>
> [1] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 5 5 5 5 5 5<br>
> .<br>
> .<br>
> .<br>
> ks.test(alelo$control, alelo$ph5.7, alelo$v100, alelo$v200,<br>
> alelo$v200ph, alelo$s100, alelo$s200, alelo$s200ph)<br>
> _________<br>
><br>
> Procuro um teste que seja apropriado para esses tipos de dados. Você sugere<br>
> um GLM?<br>
> Espero que tenha ficado mais claro agora.<br>
> Obrigada pela ajuda.<br>
> Polliana.<br>
><br>
><br>
><br>
> Em 30 de maio de 2013 14:16, Benilton Carvalho <<a href="mailto:beniltoncarvalho@gmail.com">beniltoncarvalho@gmail.com</a>><br>
> escreveu:<br>
><br>
>> Nem todo mundo entende que o que voce esta' olhando e' a influencia de<br>
>> um organismo noutro e mensurando este efeito com o numero de<br>
>> germinacoes por dia.<br>
>><br>
>> Para que voce obtenha uma ajuda mais util, voce precisa descrever os<br>
>> dados. Apresentar, pelo menos, uma amostra deles (mas nao envie<br>
>> arquivos anexados para a lista) e explicar quais sao as comparacoes<br>
>> que voce deseja fazer.<br>
>><br>
>> Eu tenho la' minhas duvidas sobre a adequacidade de Kolmogorov-Smirnov<br>
>> para isso... principalmente se vc tiver mais que 2 grupos. Talvez,<br>
>> voce queira olhar a familia de GLM e suas extensoes para<br>
>> "zero-inflated cases" (Poisson, Binomial Negativa, etc)... Mas sem a<br>
>> apresentacao precisa dos seus dados (e tambem esquema de amostragem e<br>
>> coleta), e' dificil alguem poder te ajudar.<br>
>><br>
>> b<br>
>><br>
>> Em 30 de maio de 2013 13:21, Polliana Zocche <<a href="mailto:farbby@gmail.com">farbby@gmail.com</a>> escreveu:<br>
>> > Olá,<br>
>> ><br>
>> > Eu estou procurando que teste seria apropriado para dados que contém<br>
>> > zeros<br>
>> > em seus vetores, vetores com valores repetidos e vetores com valores<br>
>> > crescentes até 100.<br>
>> ><br>
>> > Eu usei o ks.test() para verificar as diferentes entre número de<br>
>> > germinações/dia entre diferentes tratamentos avaliando alelopatia, mas o<br>
>> > teste retorna um warning por causa dos valores repetidos (ties):<br>
>> ><br>
>> > "Mensagens de aviso perdidas:<br>
>> > In ks.test(alelo$control, alelo$ph5.7, alelo$v100, alelo$v200,<br>
>> > alelo$v200ph,<br>
>> > :<br>
>> > cannot compute exact p-values with ties"<br>
>> ><br>
>> > Alguém tem uma sugestão de que teste seria apropriado para fazer<br>
>> > comparações<br>
>> > com dados dessa natureza?<br>
>> > Muito obrigada!<br>
>> ><br>
>> > --<br>
>> > Polliana Zocche de Souza<br>
>> > Bióloga/Mestre em Ecologia<br>
>> > Doutoranda em Ecologia<br>
>> > Departamento de Biologia Vegetal, IB, UNICAMP<br>
>> ><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<br>
>> > código<br>
>> > mínimo reproduzível.<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<br>
>> código mínimo reproduzível.<br>
><br>
><br>
><br>
><br>
> --<br>
> Polliana Zocche de Souza<br>
> Bióloga/Mestre em Ecologia<br>
> Doutoranda em Ecologia<br>
> Departamento de Biologia Vegetal, IB, UNICAMP<br>
><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<br>
> mínimo reproduzível.<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.</div></div></blockquote></div>
</div></div>