Pregunta La varianza de una distribución binomial negativa (NB) es siempre mayor que su media. Cuando la media de una muestra es mayor que su varianza, el intento de ajustar los parámetros de un NB con la máxima probabilidad o con la estimación de momento fallará (no hay solución con parámetros...
Tengo el siguiente histograma de datos de conteo. Y me gustaría ajustarle una distribución discreta. No estoy seguro de cómo debo hacer esto. ¿Debería superponer primero una distribución discreta, digamos Distribución binomial negativa, en el histograma para obtener los parámetros de la...
¿Existe tal paquete que proporcione la estimación del modelo de efectos mixtos binomiales negativos inflados a cero en R? Con eso quiero decir: Inflación cero donde puede especificar el modelo binomial para inflación cero, como en la función zeroinfl en el paquete pscl: zeroinfl (y ~ X | Z, dist...
¿Cuál es la estrategia adecuada para decidir qué modelo usar con los datos de conteo? Tengo datos de recuento que necesito modelar como modelo multinivel y me recomendaron (en este sitio) que la mejor manera de hacerlo es a través de errores o MCMCglmm. Sin embargo, todavía estoy tratando de...
A menudo he visto consejos para verificar si un ajuste del modelo de Poisson se dispersa en exceso o no, lo que implica dividir la desviación residual entre los grados de libertad. La relación resultante debe ser "aproximadamente 1". La pregunta es de qué rango estamos hablando para "aproximados":...
La distribución binomial negativa se ha convertido en un modelo popular para los datos de recuento (específicamente el número esperado de lecturas de secuenciación dentro de una región dada del genoma de un experimento dado) en bioinformática. Las explicaciones varían: Algunos lo explican como...
Una distribución de Poisson puede medir eventos por unidad de tiempo, y el parámetro es λλ\lambda . La distribución exponencial mide el tiempo hasta el próximo evento, con el parámetro 1λ1λ\frac{1}{\lambda} . Uno puede convertir una distribución en la otra, dependiendo de si es más fácil modelar...
¿Cuál es el número esperado de veces que debes tirar un dado hasta que cada lado haya aparecido 3 veces? Esta pregunta se hizo en la escuela primaria de Nueva Zelanda y se resolvió mediante simulaciones. ¿Cuál es la solución analítica para este
Acabo de ejecutar un GLM binomial negativo y esta es la salida: 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...
Realicé un glm.nb por glm1<-glm.nb(x~factor(group)) siendo el grupo una categoría y x una variable métrica. Cuando intento obtener el resumen de los resultados, obtengo resultados ligeramente diferentes, dependiendo de si uso summary()o summary.glm. summary(glm1)me da ... Coefficients:...
Entonces, quiero ajustar un modelo binomial negativo de efectos aleatorios. Para tal modelo, STATA puede producir coeficientes exponenciados. Según el archivo de ayuda, dichos coeficientes se pueden interpretar como relaciones de tasa de incidencia. Desafortunadamente, no soy un hablante nativo de...
Tengo un conjunto de datos muy pequeño sobre la abundancia de abejas solitarias que estoy teniendo problemas para analizar. Son datos de recuento, y casi todos los recuentos están en un tratamiento con la mayoría de los ceros en el otro tratamiento. También hay un par de valores muy altos (uno en...
Me he encontrado con una serie de preguntas prácticas al modelar datos de conteo de investigaciones experimentales utilizando un experimento dentro del sujeto. Describo brevemente el experimento, los datos y lo que he hecho hasta ahora, seguido de mis preguntas. Se mostraron cuatro películas...
Me he encontrado con tres propuestas para tratar la sobredispersión en una variable de respuesta de Poisson y un modelo de inicio de efectos fijos: Use un cuasi modelo; Use GLM binomial negativo; Use un modelo mixto con un efecto aleatorio a nivel de sujeto. ¿Pero cuál elegir realmente y por...
Tengo un conjunto de datos de series temporales en el que estoy tratando de ajustar un Modelo de Markov Oculto (HMM) para estimar el número de estados latentes en los datos. Mi pseudo código para hacer esto es el siguiente: for( i in 2 : max_number_of_states ){ ... calculate HMM with i states...
Estoy tratando de modelar las intensidades medias de parásitos que afectan a un huésped en R usando un modelo binomial negativo. Sigo recibiendo 50 o más advertencias que dicen: In dpois(y, mu, log = TRUE) : non-integer x = 251.529000 ¿Cómo puedo lidiar con esto? Mi código se ve así:...
La pregunta es la siguiente: Se recoge una muestra aleatoria de n valores de una distribución binomial negativa con el parámetro k = 3. Encuentre el estimador de máxima verosimilitud del parámetro π. Encuentre una fórmula asintótica para el error estándar de este estimador. Explique por...
En este documento (disponible gratuitamente a través de PubMed central), los autores utilizan la regresión binomial negativa para modelar la puntuación en un instrumento de detección de 10 ítems con una puntuación de 0-40. Este procedimiento asume datos de recuento, que claramente no es el caso...
He ejecutado modelos binomiales negativos y cuasi-Poisson basados en un enfoque de prueba de hipótesis. Mis modelos finales que usan ambos métodos tienen diferentes covariables e interacciones. Parece que no hay patrones cuando trazo mis residuos en ambos casos. Por lo tanto, me preguntaba qué...
Tengo un GLMM de la forma: lmer(present? ~ factor1 + factor2 + continuous + factor1*continuous + (1 | factor3), family=binomial) Cuando lo uso drop1(model, test="Chi"), obtengo resultados diferentes a los que uso Anova(model, type="III")del paquete del automóvil o summary(model). Estos dos...