Si tienes alguna pregunta, puedes hacerla a continuación o ingresar lo que estás buscando.
Construí varios modelos glmer.nb con diferentes combinaciones de intercepciones aleatorias, y para uno de los modelos (con intercepciones anidadas aleatorias, con el AICc más bajo), consistently obtengo: “límite de iteración alcanzado”, sin la advertencia habitual: “En theta.ml(Y, mu, weights = object@resp$weights, limit = limit, :…”. Esto es lo que sé: . . . Read more
Me gustaría trazar las estimaciones de los efectos fijos de un modelo lme4::lmer. Actualmente, la salida de sjPlot::plot_model incluye dos gráficos adicionales que no quiero: SD (Observations) y SD (Intercept). ¿Cómo puedo eliminar estos dos gráficos adicionales? library(lme4) data(iris) lmerMod <- lmer(Sepal.Length ~ Sepal.Width * Petal.Length + (1|Species), data=iris) library(sjPlot) . . . Read more
No estoy seguro de qué hacer en la siguiente situación, donde faltan algunos niveles del efecto fijo (dentro de un efecto aleatorio) – están desequilibrados. Imagina un acuario con 5.000 peces individuales. Forman parte de 100 especies diferentes. Quiero comprobar si hay una relación entre su peso (continuo) y si . . . Read more
Me gustaría ajustar un modelo de efectos mixtos y obtener una lista de los diez principales modelos candidatos clasificados por AIC. Para modelos de efectos fijos, utilizo la función glmulti:glmulti para ajustar el modelo y la función glmulti:weightable para obtener una lista de modelos candidatos clasificados por AIC, lo que . . . Read more
Me gustaría calcular los intervalos de confianza al 95% para la predicción en los GLMM de distribución binomial negativa, utilizando el paquete lme4. Sin embargo, debido a la falta de una opción para calcular los errores estándar de las predicciones por defecto, intento lo siguiente: #Paquete library(lme4) #Crear un conjunto . . . Read more