[R-br] GLM model

Cesar Rabak cesar.rabak em gmail.com
Sex Nov 9 16:21:17 -02 2018


A resposta está no 1º parágrafo da *vignette* para o pacote quasi:  «
Computing “quasi-AIC” (QAIC), in R is a minor pain, because the R Core team
(or at least the ones who wrote glm, glmmPQL, etc.) are purists and don’t
believe that quasi- models should report a likelihood. »


HTH
--
Cesar Rabak


On Thu, Nov 8, 2018 at 9:34 PM Lucas Stempkowski por (R-br) <
r-br em listas.c3sl.ufpr.br> wrote:

> Olá boa noite,
>
> Estou utilizando um GLM para análise de incidência de doença. Utilizei o
> seguinte comando:
>
> Call:
>
> *glm(formula = incid ~ factor(bloco) + factor(cod) + cultivar + *
> *    Amostragem, family = quasibinomial("probit"), data = dados)*
>
> e obtive os resultados a seguir:
>
> Minha dúvida fica em relação ao* AIC: NA, *há algo errado no modelo usado?
>
> Deviance Residuals:
>      Min        1Q    Median        3Q       Max
> -0.62644  -0.17630  -0.07311   0.07820   0.85815
>
> Coefficients:
>                  Estimate Std. Error t value Pr(>|t|)
> (Intercept)     -0.717767   0.107812  -6.658 1.40e-10 ***
> factor(bloco)2  -0.002541   0.087289  -0.029   0.9768
> factor(bloco)3  -0.093526   0.089447  -1.046   0.2966
> factor(bloco)4  -0.051161   0.088405  -0.579   0.5632
> factor(cod)2    -0.153266   0.083405  -1.838   0.0671 .
> factor(cod)3     0.027784   0.079826   0.348   0.7280
> factor(cod)4    -1.132366   0.135880  -8.334 3.21e-15 ***
> factor(cod)5    -1.266976   0.151660  -8.354 2.79e-15 ***
> cultivarParrudo -0.752089   0.084474  -8.903  < 2e-16 ***
> cultivarReponte -0.331743   0.071247  -4.656 4.92e-06 ***
> Amostragem      -0.038367   0.022256  -1.724   0.0858 .
> ---
> Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
>
> (Dispersion parameter for quasibinomial family taken to be 0.06977337)
>
>     Null deviance: 39.898  on 299  degrees of freedom
> Residual deviance: 18.381  on 289  degrees of freedom
> AIC: NA
>
> Number of Fisher Scoring iterations: 7
>
>
> Minha dúvida fica em relação ao* AIC: NA, *há algo errado no modelo
> usado?
>
>
> *--*
>
> *-----------------------------------------------------*
>
> Lucas Antonio Stempkowski
>
> Eng. Agronômo - UFFS
>
> Mestrando em Produção Vegetal CAV/UDESC
>
> +55 54 99153-0124/ 98448-4633
>
> --------------------------------------------------
>
>
>
>
>
> <http://www.avg.com/email-signature?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail> Livre
> de vírus. www.avg.com
> <http://www.avg.com/email-signature?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail>.
> <#m_-8074017122499481837_DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>
> _______________________________________________
> 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.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20181109/f7d28db5/attachment.html>


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