<div dir="ltr">Ana Paula,<div><br></div><div>Como você já pôde ver nas respostas sobre sua dúvida (do ponto de vista da técnica do uso da GLM e da sintaxe do R), acho que se enseja também uma outra discussão:</div><div><br>
</div><div>Qual é o sentido para sua modelagem de ter todos os cultivares (e tempos) como níveis (removendo o intercepto)?</div><div><br></div><div>O uso de um nível de cada fator como "referência" me parece mais lógico na análise que um modelo onde cada nível seja explicitamente estimado. senão vejamos o seguinte: qual é a resposta que a significância do nível que na primeira abordagem estava na referência?</div>
<div><br></div><div>Ou a outra pergunta: o quê significa (sem trocadilho, por favor) a não significância estatística para o fator tempo na sua regressão?</div><div><br></div><div>Note também na sua regressão os resultados para o intercepto, que geralmente são ignorados nesse tipo de análise, mas neste caso à guisa de ilustração vamos considerar o resultado: o p-valor para o intercepto (que considera T0 e C1 como referência) é <span style="color:rgb(80,0,80);font-family:arial,sans-serif;font-size:13px">0,98628 </span>suficientemente grande para dizer estatisticamente o quê?</div>
<div><span style="color:rgb(80,0,80);font-family:arial,sans-serif;font-size:13px"><br></span></div><div><span style="color:rgb(80,0,80);font-family:arial,sans-serif;font-size:13px">Uma outra coisa que gostaria de saber é o número de casos que você dispõe para fazer essa regressão?</span></div>
<div><span style="color:rgb(80,0,80);font-family:arial,sans-serif;font-size:13px"> </span></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">2013/12/13 ana paula coelho madeira <span dir="ltr"><<a href="mailto:apcmadeira@hotmail.com" target="_blank">apcmadeira@hotmail.com</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div><div dir="ltr">Entendi, muito obrigada! Mas e quanto aos fatores que não aparecem na análise? Tenho 8 tempos e 4 cultivares e quando peço o summary do modelo o tempo 0 e a cultivar C1 não aparecem na análise. O que ocorreu? <br>
Como avaliar a significância desses fatores?<br>Segue a saída:<div class="im"><br>NP1 <- glm (Y ~ T+C, family = poisson, data=dados)<br></div><div class="im">> Outra coisa, quando rodo o NP1 tenho a seguinte saída:<br>
> <br>> Deviance Residuals:<br>> Min 1Q Median 3Q Max <br>> -5.7588 -1.0631 -0.0002 0.7864 4.7855 <br>> Coefficients:<br>> Estimate Std. Error z value Pr(>|z|) <br>
> (Intercept) -1.718e+01 9.994e+02 -0.017 0.98628 <br>> T15 3.932e-09 1.413e+03 0.000 1.00000 <br>> T29 1.951e+01 9.994e+02 0.020 0.98442 <br>> T43 2.088e+01 9.994e+02 0.021 0.98333 <br>
> T57 2.139e+01 9.994e+02 0.021 0.98292 <br>> T71 2.140e+01 9.994e+02 0.021 0.98291 <br>> T85 2.081e+01 9.994e+02 0.021 0.98339 <br>> T99 1.745e+01 9.994e+02 0.017 0.98607 <br>
> C2 -1.468e-01 5.659e-02 -2.594 0.00950 **<br>> C3 -1.433e-01 5.653e-02 -<a href="tel:2.535%C2%A0%200.01123" value="+12535001123" target="_blank">2.535 0.01123</a> *<br>> C4 -1.800e-01 5.710e-02 -<a href="tel:3.152%C2%A0%200.00162" value="+13152000162" target="_blank">3.152 0.00162</a> **<br>
<br></div>Obrigada!<br> <br><div>Date: Fri, 13 Dec 2013 15:35:26 -0200<br>From: <a href="mailto:paulojus@leg.ufpr.br" target="_blank">paulojus@leg.ufpr.br</a><br>To: <a href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">r-br@listas.c3sl.ufpr.br</a><br>
Subject: Re: [R-br] GLM<div><div class="h5"><br><br><pre>Nao!<br> <br>o ponto "." siginifica todas as variaveis no data-frame<br>alem da resposta (Y)<br> <br>Se o data frame contem apenas Y, T, e C entao as especificaqcoes sao <br>
iguais<br> <br>On Fri, 13 Dec 2013, Fátima Lima Paula wrote:<br> <br>> Ana, o modelo com o ponto é o modelo nulo.<br>> O outro está ajustando pelo T e C.<br>> Não aparecem as categorias que são as de referência.<br>
> Abs<br>> <br>> <br>> Em 13 de dezembro de 2013 13:26, ana paula coelho madeira <<a href="mailto:apcmadeira@hotmail.com" target="_blank">apcmadeira@hotmail.com</a>> escreveu:<br>> Boa tarde pessoal.<br>
> <br>> Tenho dados de contagem em função do tempo e da cultivar. Tenho duvida na maneira como declaro o modelo. Qual é a diferença entre os modelos<br>> <br>> NP1 <- glm (Y ~ T+C, family = poisson, data=dados)<br>
> e<br>> NP2 <- glm (Y ~ ., family = poisson, data=dados).<br>> <br>> Outra coisa, quando rodo o NP1 tenho a seguinte saída:<br>> <br>> Deviance Residuals:<br>> Min 1Q Median 3Q Max <br>
> -5.7588 -1.0631 -0.0002 0.7864 4.7855 <br>> Coefficients:<br>> Estimate Std. Error z value Pr(>|z|) <br>> (Intercept) -1.718e+01 9.994e+02 -0.017 0.98628 <br>
> T15 3.932e-09 1.413e+03 0.000 1.00000 <br>> T29 1.951e+01 9.994e+02 0.020 0.98442 <br>> T43 2.088e+01 9.994e+02 0.021 0.98333 <br>> T57 2.139e+01 9.994e+02 0.021 0.98292 <br>
> T71 2.140e+01 9.994e+02 0.021 0.98291 <br>> T85 2.081e+01 9.994e+02 0.021 0.98339 <br>> T99 1.745e+01 9.994e+02 0.017 0.98607 <br>> C2 -1.468e-01 5.659e-02 -2.594 0.00950 **<br>
> C3 -1.433e-01 5.653e-02 -<a href="tel:2.535%C2%A0%200.01123" value="+12535001123" target="_blank">2.535 0.01123</a> *<br>> C4 -1.800e-01 5.710e-02 -<a href="tel:3.152%C2%A0%200.00162" value="+13152000162" target="_blank">3.152 0.00162</a> **<br>
><br>> Pegunto: Não aparece o efeito da cultivar 1 (C1) e do tempo 0 (T0), por que isso ocorre? Como falar sobre o efeito desses dois fatores?<br>> <br>> Desde já agradeço.<br>> <br>
> Ana<br>> <br>> <br>> <br>> _______________________________________________<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.<br>
> <br>> <br>> <br>> <br>> --<br>> "Minha felicidade depende da qualidade dos meus pensamentos"<br>> <br>><br></pre><br></div></div>_______________________________________________
R-br mailing list
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a>
<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>
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></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>