Al modelar datos de recuento de reclamos en un entorno de seguro, comencé con Poisson pero luego noté una sobredispersión. Un Cuasi-Poisson mejor modeló la mayor relación media-varianza que el Poisson básico, pero noté que los coeficientes eran idénticos en los modelos Poisson y Cuasi-Poisson. Si...
Estoy realizando un modelo de regresión de Poisson con 1 variable de respuesta y 6 covariables. La selección del modelo usando AIC da como resultado un modelo con todas las covariables, así como 6 términos de interacción. Sin embargo, el BIC da como resultado un modelo con solo 2 covariables y sin...
Estoy llevando a cabo una regresión de Poisson con el objetivo final de comparar (y tomar la diferencia de) los recuentos medios pronosticados entre dos niveles de factores en mi modelo: , mientras otras covariables modelo (que son todas binarias) constantes. Me preguntaba si alguien podría...
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...
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 gustaría probar en qué regresión se ajusta mejor a mis datos. Mi variable dependiente es un recuento y tiene muchos ceros. Y necesitaría ayuda para determinar qué modelo y familia usar (poisson o cuasipoisson, o regresión de poisson inflada a cero) y cómo probar los supuestos. Regresión de...
El mgcvpaquete Rtiene dos funciones para ajustar las interacciones del producto tensorial: te()y ti(). Entiendo la división básica del trabajo entre los dos (ajustar una interacción no lineal versus descomponer esta interacción en efectos principales y una interacción). Lo que no entiendo es por...
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é...
Estoy analizando los datos de seguimiento ocular de un experimento diseñado. Una versión simplificada de mis datos se ve así (puede obtener los datos dput () aquí ), head(lookDATA) participant fixationImage fixationCount 1 9 Automobile 81 2 9 Bird 63 3 9 Chair 82 4 9 Dog 64 5 9 Face 90 6 9 Plant...
Estoy tratando de meterme un poco en las estadísticas, pero estoy atascado con algo. Mis datos son los siguientes: Year Number_of_genes 1990 1 1991 1 1993 3 1995 4 Ahora quiero construir un modelo de regresión para poder predecir la cantidad de genes para un año determinado con base en los...
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...
Estoy interesado en estimar una razón de riesgo ajustada, análoga a cómo se estima una razón de probabilidad ajustada usando regresión logística. Parte de la literatura (p. Ej., Esto ) indica que usar la regresión de Poisson con los errores estándar de Huber-White es una forma basada en el modelo...
Zuur 2013 Beginners Guide to GLM & GLMM sugiere validar una regresión de Poisson trazando los residuos de Pearson contra los valores ajustados. Zuur afirma que no deberíamos ver los residuos desplegándose a medida que aumentan los valores ajustados, como el gráfico adjunto (dibujado a...
He visto muchas cosas que discuten si una regresión de Poisson básica es una versión anidada de una regresión de Poisson inflada a cero. Por ejemplo, este sitio argumenta que lo es, ya que este último incluye parámetros adicionales para modelar ceros adicionales, pero por lo demás incluye los...
Tengo estos datos: set.seed(1) predictor <- rnorm(20) set.seed(1) counts <- c(sample(1:1000, 20)) df <- data.frame(counts, predictor) Corrí una regresión de Poisson poisson_counts <- glm(counts ~ predictor, data = df, family = "poisson") Y una regresión binomial...
Quiero calcular las tasas de incidencia para presentar a lo largo de las razones de riesgo para presentar medidas de riesgo relativas y absolutas. Vi en otros estudios que tales tasas de incidencia pueden calcularse utilizando modelos de Poisson con tiempo de seguimiento en el modelo como...
Estoy trabajando en modelos de datos de recuento cero inflados usando el psclpaquete. ¡Me pregunto por qué no hay desarrollo de modelos para modelos de datos de conteo inflados! Además, ¿por qué no hay desarrollo de modelos de datos de conteo bimodales, digamos cero y 2 inflados! Una vez que generé...
Estoy examinando las tendencias (entre 1998 y 2011) en las tasas de mortalidad entre pacientes con enfermedad de Crohn. Cada paciente (caso) se incluyó entre 1998 y 2011. En el momento de la inclusión, cada paciente se comparó con un control saludable con la misma edad y sexo. Estoy analizando las...
En la "Guía de un practicante de modelos lineales generalizados" en el párrafo 1.83 se afirma que: "En el caso particular de un GLM multiplicativo de Poisson, se puede demostrar que el reclamo de modelado cuenta con un término de compensación igual al logaritmo de la exposición que produjo...