Questão A variação de uma distribuição binomial negativa (NB) é sempre maior que sua média. Quando a média de uma amostra é maior que sua variância, a tentativa de ajustar os parâmetros de um RN com probabilidade máxima ou com estimativa de momento falhará (não há solução com parâmetros...
Eu tenho o seguinte histograma de dados de contagem. E eu gostaria de ajustar uma distribuição discreta a ela. Não tenho certeza de como devo fazer isso. Devo primeiro sobrepor uma distribuição discreta, digamos, Distribuição binomial negativa, no histograma para obter os parâmetros da...
Existe um pacote que preveja a estimativa do modelo de efeitos mistos binomiais negativos inflados a zero em R? Com isso quero dizer: Inflação zero, onde você pode especificar o modelo binomial para inflação zero, como na função zeroinfl no pacote pscl: zeroinfl (y ~ X | Z, dist = "negbin") onde...
Qual é a estratégia apropriada para decidir qual modelo usar com os dados de contagem? Tenho dados de contagem que preciso modelar como modelo multinível e me foi recomendado (neste site) que a melhor maneira de fazer isso é através de bugs ou MCMCglmm. No entanto, ainda estou tentando aprender...
Eu sempre vi o conselho para verificar se um ajuste do modelo de Poisson está ou não disperso em excesso, envolvendo a divisão do desvio residual pelos graus de liberdade. A proporção resultante deve ser "aproximadamente 1". A questão é de qual faixa estamos falando para "aproximado" - qual é a...
A distribuição binomial negativa tornou-se um modelo popular para dados de contagem (especificamente o número esperado de leituras de seqüenciamento dentro de uma determinada região do genoma de um determinado experimento) em bioinformática. As explicações variam: Alguns o explicam como algo que...
Uma distribuição Poisson pode medir eventos por unidade de tempo, e o parâmetro é λλ\lambda . A distribuição exponencial mede o tempo até o próximo evento, com o parâmetro 1λ1λ\frac{1}{\lambda} . Pode-se converter uma distribuição na outra, dependendo se é mais fácil modelar eventos ou...
Qual é o número esperado de vezes que você deve rolar um dado até que cada lado apareça 3 vezes? Esta pergunta foi feita na escola primária da Nova Zelândia e foi resolvida usando simulações. Qual é a solução analítica para esse
Acabei de executar um GLM binomial negativo e esta é a saída: Call: glm.nb(formula = small ~ method + site + depth, data = size.dat, init.theta = 1.080668549, link = log) Deviance Residuals: Min 1Q Median 3Q Max -2.2452 -0.9973 -0.3028 0.3864 1.8727 Coefficients: Estimate Std. Error z...
Eu conduzi um glm.nb por glm1<-glm.nb(x~factor(group)) sendo o grupo uma categoria c e x uma variável métrica. Quando tento obter o resumo dos resultados, obtenho resultados ligeiramente diferentes, dependendo se uso summary()ou nãosummary.glm . summary(glm1)me dá ... Coefficients:...
Então, eu quero ajustar um modelo negativo-binomial de efeitos aleatórios. Para esse modelo, o STATA pode produzir coeficientes exponenciados. De acordo com o arquivo de ajuda, esses coeficientes podem ser interpretados como taxas de incidência. Infelizmente, eu não sou um falante nativo de inglês...
Eu tenho um conjunto de dados muito pequeno sobre a abundância de abelhas solitárias que estou tendo problemas para analisar. São dados de contagem, e quase todas as contagens estão em um tratamento com a maioria dos zeros no outro tratamento. Também existem alguns valores muito altos (um em cada...
Eu me deparei com várias questões práticas ao modelar dados de contagem de pesquisas experimentais usando um experimento dentro do assunto. Descrevo brevemente o experimento, os dados e o que fiz até agora, seguidos pelas minhas perguntas. Quatro filmes diferentes foram exibidos para uma amostra...
Encontrei três propostas para lidar com a super-dispersão em uma variável de resposta de Poisson e um modelo inicial de efeitos fixos: Use um modelo quase; Use GLM binomial negativo; Use um modelo misto com um efeito aleatório no nível de assunto. Mas qual escolher, e por quê? Existe algum...
Eu tenho um conjunto de dados de séries temporais no qual estou tentando ajustar um Modelo de Markov oculto (HMM) para estimar o número de estados latentes nos dados. Meu pseudo-código para fazer isso é o seguinte: for( i in 2 : max_number_of_states ){ ... calculate HMM with i states ......
Estou tentando modelar as intensidades médias de parasitas que afetam um hospedeiro em R usando um modelo binomial negativo. Eu recebo 50 ou mais avisos que dizem: In dpois(y, mu, log = TRUE) : non-integer x = 251.529000 Como posso lidar com isso? Meu código fica assim: mst.nb =...
A questão é a seguinte: Uma amostra aleatória de n valores é coletada de uma distribuição binomial negativa com o parâmetro k = 3. Encontre o estimador de probabilidade máxima do parâmetro π. Encontre uma fórmula assintótica para o erro padrão deste estimador. Explique por que a...
No presente trabalho (disponível gratuitamente via PubMed Central), os autores usam regressão binomial negativa para modelar a pontuação em um instrumento de triagem de 10 itens marcou 0-40. Este procedimento assume dados de contagem, o que claramente não é o caso aqui. Gostaria de suas opiniões...
Corri modelos binomiais negativos e quase-Poisson com base em uma abordagem de teste de hipóteses. Meus modelos finais, usando os dois métodos, têm covariáveis e interações diferentes. Parece que não há padrões quando planto meus resíduos nos dois casos. Assim, eu queria saber qual teste eu...
Eu tenho um GLMM do formulário: lmer(present? ~ factor1 + factor2 + continuous + factor1*continuous + (1 | factor3), family=binomial) Quando uso drop1(model, test="Chi"), obtenho resultados diferentes dos que utilizo Anova(model, type="III")na embalagem do carro ou summary(model). Estes dois...