Tag: MIXED-MODELS
En Tutz y Schmid (2016), ejemplo 9.4 (Cáncer de mama) comparan varios modelos de supervivencia. Entre ellos, hay un modelo de “efectos fijos” para peligros heterogéneos, donde los efectos fijos están limitados por una penalización que toma esta forma El breve ejemplo menciona que utilizan mgcv::gam con la penalización en . . . Read more
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 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