¿Existe una razón intuitiva para que los efectos aleatorios se reduzcan a su valor esperado en el modelo mixto lineal
Estructura de datos > str(data) 'data.frame': 6138 obs. of 10 variables: $ RT : int 484 391 422 516 563 531 406 500 516 578 ... $ ASCORE : num 5.1 4 3.8 2.6 2.7 6.5 4.9 2.9 2.6 7.2 ... $ HSCORE : num 6 2.1 7.9 1 6.9 8.9 8.2 3.6 1.7 8.6 ... $ MVMNT : Factor w/ 2 levels "_Withd","Appr": 2 2...
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 tratando de adaptar un modelo mixto multivariado (es decir, respuesta múltiple) R. Aparte de los paquetes ASReml-ry SabreR(que requieren software externo), parece que esto solo es posible en MCMCglmm. En el documento que acompaña al MCMCglmmpaquete (págs. 6), Jarrod Hadfield describe el...
Dadas tres variables, yy x, que son positivas continuas, y z, que es categórica, tengo dos modelos candidatos dados por: fit.me <- lmer( y ~ 1 + x + ( 1 + x | factor(z) ) ) y fit.fe <- lm( y ~ 1 + x ) Espero comparar estos modelos para determinar qué modelo es más apropiado. Me parece que...
¿Alguien sabe cómo calcular (o extraer) el apalancamiento y las distancias de Cook para un merobjeto de clase (obtenido a través del lme4paquete)? Me gustaría trazar estos para un análisis de
Entiendo el efecto de un efecto aleatorio categórico en un modelo de efectos mixtos en el sentido de que realiza una agrupación parcial de las observaciones por nivel en el efecto aleatorio, suponiendo efectivamente que las observaciones no son independientes en sí mismas, sino solo sus...
Estoy trabajando en un conjunto de datos. Después de usar algunas técnicas de identificación de modelos, obtuve un modelo ARIMA (0,2,1). Utilicé la detectIOfunción en el paquete TSAen R para detectar un valor atípico innovador (IO) en la observación número 48 de mi conjunto de datos...
Siempre me han enseñado que los efectos aleatorios solo influyen en la varianza (error), y que los efectos fijos solo influyen en la media. Pero he encontrado un ejemplo en el que los efectos aleatorios también influyen en la media: la estimación del coeficiente: require(nlme) set.seed(128) n...
Hay varios hilos en R-sig-ME sobre cómo obtener intervalos de confianza para las predicciones usando lme4y nlmeen R. Por ejemplo, aquí y aquí en 2010, incluidos algunos comentarios de Dougals Bates, uno de los autores de ambos paquetes. Dudo en citarlo literalmente, por temor a que sean sacados de...
Esta pregunta / tema surgió en una discusión con un colega y estaba buscando algunas opiniones sobre esto: Estoy modelando algunos datos usando una regresión logística de efectos aleatorios, más precisamente una regresión logística de intercepción aleatoria. Para los efectos fijos tengo 9...
En el modelo mixto, suponemos que los efectos aleatorios (parámetros) son variables aleatorias que siguen distribuciones normales. Se ve muy similar al método bayesiano, en el que se supone que todos los parámetros son aleatorios. Entonces, ¿es el modelo de efectos aleatorios un tipo de caso...
Estoy analizando los resultados de un experimento de tiempo de reacción en R. Ejecuté un ANOVA de medidas repetidas (1 factor dentro del sujeto con 2 niveles y 1 factor entre sujetos con 2 niveles). Ejecuté un modelo lineal mixto similar y quería resumir los resultados de lmer en forma de tabla...
¿Cuál es la estructura de varianza-covarianza predeterminada para efectos aleatorios en glmero lmeren lme4paquete? ¿Cómo se especifica otra estructura de varianza-covarianza para efectos aleatorios en el código? No pude encontrar ninguna información al respecto en la
Estoy analizando algunos datos de comportamiento utilizando lme4en R, sobre todo después de excelentes tutoriales de Bodo invierno , pero no entiendo si yo estoy manejando adecuadamente las interacciones. Peor aún, nadie más involucrado en esta investigación utiliza modelos mixtos, por lo que estoy...
Tengo el siguiente modelo m_plotequipado con lme4::lmerefectos aleatorios cruzados para los participantes ( lfdn) y los elementos ( content): Random effects: Groups Name Variance Std.Dev. Corr lfdn (Intercept) 172.173 13.121 role1 62.351 7.896 0.03 inference1 24.640 4.964 0.08 -0.30...
Antes de hacer esta pregunta, busqué en nuestro sitio y encontré muchas preguntas similares (como aquí , aquí y aquí ). Pero creo que esas preguntas relacionadas no fueron bien respondidas o discutidas, por lo que me gustaría plantear esta pregunta nuevamente. Siento que debería haber una gran...
Sinceramente espero haber formulado esta pregunta de tal manera que pueda responderse definitivamente; si no, ¡hágamelo saber y volveré a intentarlo! También debería adivinar que utilizaré R para estos análisis. Tengo varias medidas de las plant performance (Ys)que sospecho que fueron...
Estoy confundido por las declaraciones en una página web de la UCLA sobre la regresión logística de efectos mixtos. Muestran una tabla de coeficientes de efectos fijos al ajustar dicho modelo y el primer párrafo a continuación parece interpretar los coeficientes exactamente como una regresión...
Considere un experimento con múltiples participantes humanos, cada uno medido varias veces en dos condiciones. Se puede formular un modelo de efectos mixtos (usando la sintaxis lme4 ) como: fit = lmer( formula = measure ~ (1|participant) + condition ) Ahora, digamos que quiero generar...