Likelihood Ratio Test para significância de fatores fixos em GLMM Poisson e BN
2 participantes
Página 1 de 1
Likelihood Ratio Test para significância de fatores fixos em GLMM Poisson e BN
Boa tarde.
Estou trabalhando em Modelos Lineares Generalizados Mistos para dados de cotnagem (Poisson e Binomial Negativa). Estou analisando o efieto de tratamento, tempo e interação no número de indivíduos e recrutas em uma floresta.
Preciso obter os valores p (ou algo do tipo) para esses fatores fixos do meu modelo. Dentre os métodos existentes, selecionei o Likelihood Ratio test. Para os modelos que ajustei, a interação deu significativa. O problema é que quando uso o comando drop1 para obter os resutlados do LRT, só aparece os valores para a interação, e não para os efeitos fixos (tratamento e tempo recebem 0) separadamente. Como faço para obter esses valores para os efeitos fixos separadamente?
Meu código:
> drop1(GLMM.R, scope = ~ treatment + year + treatment*year, test = "Chisq")
Single term deletions
Model:
recruits ~ treatment + year + treatment * year + (1 | plot)
npar AIC LRT Pr(Chi)
<none> 111.19
treatment 0 111.19 0.000
year 0 111.19 0.000
treatment:year 4 196.60 93.418 < 2.2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Estou trabalhando em Modelos Lineares Generalizados Mistos para dados de cotnagem (Poisson e Binomial Negativa). Estou analisando o efieto de tratamento, tempo e interação no número de indivíduos e recrutas em uma floresta.
Preciso obter os valores p (ou algo do tipo) para esses fatores fixos do meu modelo. Dentre os métodos existentes, selecionei o Likelihood Ratio test. Para os modelos que ajustei, a interação deu significativa. O problema é que quando uso o comando drop1 para obter os resutlados do LRT, só aparece os valores para a interação, e não para os efeitos fixos (tratamento e tempo recebem 0) separadamente. Como faço para obter esses valores para os efeitos fixos separadamente?
Meu código:
> drop1(GLMM.R, scope = ~ treatment + year + treatment*year, test = "Chisq")
Single term deletions
Model:
recruits ~ treatment + year + treatment * year + (1 | plot)
npar AIC LRT Pr(Chi)
<none> 111.19
treatment 0 111.19 0.000
year 0 111.19 0.000
treatment:year 4 196.60 93.418 < 2.2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Carlos Delano C. Olvieira- Mensagens : 6
Data de inscrição : 30/03/2020
Valeria Andrade- Mensagens : 94
Data de inscrição : 12/03/2018
Tópicos semelhantes
» Pacote que para função do envelope quase-poisson
» GLMM - Fixed effect para todas as variáveis
» Calcular Total, Média e SD para mais de uma variável e fatores
» Módulo IV dúvida de terminologia para o test t
» como calcular sobredispersão - GLMM
» GLMM - Fixed effect para todas as variáveis
» Calcular Total, Média e SD para mais de uma variável e fatores
» Módulo IV dúvida de terminologia para o test t
» como calcular sobredispersão - GLMM
Página 1 de 1
Permissões neste sub-fórum
Não podes responder a tópicos
|
|