Preguntas etiquetadas con lme4-nlme

Por un simple ejemplo, suponga que hay dos modelos de regresión lineal Modelo 1 tiene tres predictores, x1a, x2b, yx2c El modelo 2 tiene tres predictores del modelo 1 y dos predictores adicionales x2ayx2b Hay una ecuación de regresión poblacional donde la varianza poblacional explicada es para...

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 tratando de calcular las predicciones de efectos aleatorios de un modelo mixto lineal a mano, y usando la notación proporcionada por Wood en modelos aditivos generalizados: una introducción con R (pg 294 / pg 307 de pdf), me estoy confundiendo sobre cada parámetro representa. A continuación...

La fórmula que uno debe especificar para entrenar un modelo multinivel (usando lmerdesde la lme4 Rbiblioteca) siempre me da. He leído innumerables libros de texto y tutoriales, pero nunca lo he entendido correctamente. Así que aquí hay un ejemplo de este tutorial que me gustaría ver formulado en...

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...

Estoy tratando de reproducir varias pruebas de interacción entre ambas lmy lmeren medidas repetidas (2x2x2). La razón por la que quiero comparar ambos métodos es porque el GLM de SPSS para medidas repetidas produce exactamente los mismos resultados que el lmenfoque presentado aquí, por lo que al...

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...

¿Cómo puedo obtener efectos aleatorios agrupados para lmer después de una imputación múltiple? Estoy usando ratones para imputar múltiples un marco de datos. Y lme4 para un modelo mixto con intercepción aleatoria y pendiente aleatoria. Agrupar lmer funciona bien, excepto que no agrupa los efectos...

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...

Estoy realizando pruebas post-hoc en un modelo lineal de efectos mixtos en R( lme4paquete). Estoy usando el multcomppaquete ( glht()función) para realizar las pruebas post-hoc. Mi diseño experimental es de medidas repetidas, con un efecto de bloqueo aleatorio. Los modelos se especifican...

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...

En la lmerfunción dentro lme4de Rhay una llamada para la construcción de una matriz de modelo de efectos aleatorios, , como se explica aquí , páginas 7 - 9.ZZZ El cálculo de implica productos de KhatriRao y / o Kronecker de dos matrices, y . ZZZJyoJyoJ_iXyoXyoX_i La matriz es un bocado: "matriz...

Como contraparte de esta publicación , trabajé en la simulación de datos con variables continuas, prestándome a intercepciones y pendientes correlacionadas. Aunque hay excelentes publicaciones sobre este tema en el sitio y fuera del sitio , tuve dificultades para encontrar un ejemplo de principio...

Entonces, tengo 16 ensayos en los que estoy tratando de autenticar a una persona de un rasgo biométrico usando Hamming Distance. Mi umbral está establecido en 3.5. Mis datos están a continuación y solo la prueba 1 es un verdadero positivo: Trial Hamming Distance 1 0.34 2 0.37 3 0.34 4 0.29 5...

Estoy ajustando un glmermodelo en el lme4paquete R. Estoy buscando una tabla anova con el valor p que se muestra allí, pero no puedo encontrar ningún paquete que se ajuste a ella. ¿Es posible hacerlo en R? El modelo que estoy ajustando es de la

Tengo algunos problemas para obtener resultados equivalentes entre un aovmodelo de medidas repetidas entre dentro y un lmermodelo mixto. Mis datos y script se ven de la siguiente manera data=read.csv("https://www.dropbox.com/s/zgle45tpyv5t781/fitness.csv?dl=1") data$id=factor(data$id) data id...

En un modelo mixto lineal generalizado logístico (familia = binomio), no sé cómo interpretar la varianza de los efectos aleatorios: Random effects: Groups Name Variance Std.Dev. HOSPITAL (Intercept) 0.4295 0.6554 Number of obs: 2275, groups: HOSPITAL, 14 ¿Cómo interpreto este resultado...

Estoy simulando ensayos de Bernoulli con un entre grupos y luego el modelo correspondiente con el paquete:logitθ ∼ N( logitθ0 0,12)logitθ∼norte(logitθ0 0,12)\text{logit}\, \theta \sim {\cal N}(\text{logit}\, \theta_0, 1^2)lme4 library(lme4) library(data.table) I <- 30 # number of groups J <-...

Estoy usando lmer en R para verificar el efecto de la condición ( cond) en algún resultado. Aquí están algunos hechos de seguridad de datos, donde s es el identificador de objeto y a, by cson condiciones. library("tidyr") library("dplyr") set.seed(123) temp <- data.frame(s = paste0("S", 1:30),...