Preguntas etiquetadas con mixed-model

Supongamos que tengo alguna variable de respuesta yijyijy_{ij} que se midió a partir del hermano jjj en la familia iii . Además, se recopilaron algunos datos de comportamiento xijxijx_{ij} al mismo tiempo de cada sujeto. Estoy tratando de analizar la situación con el siguiente modelo lineal de efectos mixtos: yij=α0+α1xij+δ1ixij+εijyij=α0+α1xij+δ1ixij+εijy_{ij} …

14 mixed-model  lme4-nlme  covariance-matrix  non-independent 

¿Puede alguien decirme cómo hacer que R calcule el punto de ruptura en un modelo lineal por partes (como un parámetro fijo o aleatorio), cuando también necesito estimar otros efectos aleatorios? He incluido un ejemplo de juguete a continuación que se ajusta a una regresión de palo de hockey / …

14 r  mixed-model  lme4-nlme  change-point  piecewise-linear 

Douglas Bates afirma que los siguientes modelos son equivalentes "si la matriz de varianza-covarianza para los efectos aleatorios con valores vectoriales tiene una forma especial, llamada simetría compuesta" ( diapositiva 91 en esta presentación ): m1 <- lmer(y ~ factor + (0 + factor|group), data) m2 <- lmer(y ~ factor …

13 r  anova  mixed-model  repeated-measures  lme4-nlme 

Veo mencionado en varios lugares que ANOVA hace su estimación utilizando el método de los momentos. Esa afirmación me confunde porque, aunque no estoy familiarizado con el método de los momentos, entiendo que es algo diferente y no equivalente al método de máxima probabilidad; Por otro lado, ANOVA puede verse …

13 anova  mixed-model  maximum-likelihood  method-of-moments 

Aprendí en las estadísticas elementales que, con un modelo lineal general, para que las inferencias sean válidas, las observaciones deben ser independientes. Cuando se produce la agrupación, la independencia ya no puede conducir a una inferencia no válida a menos que esto se tenga en cuenta. Una forma de dar …

13 r  mixed-model  inference  independence 

Consideremos este conjunto de datos hipotético: set.seed(12345) num.subjects <- 10 dose <- rep(c(1,10,50,100), num.subjects) subject <- rep(1:num.subjects, each=4) group <- rep(1:2, each=num.subjects/2*4) response <- dose*dose/10 * group + rnorm(length(dose), 50, 30) df <- data.frame(dose=dose, response=response, subject=subject, group=group) podemos usar lmepara modelar la respuesta con un modelo de efectos aleatorios: require(nlme) …

13 mixed-model  linear-model  prediction 

Supongamos que he medido alguna variable en hermanos, que están anidados dentro de las familias. La estructura de datos se ve así: valor de hermano familiar ------ ------- ----- 1 1 y_11 1 2 y_12 2 1 y_21 2 2 y_22 2 3 y_23 ... ... ... Quiero saber la …

13 mixed-model  intraclass-correlation 

¿Alguien puede decirme la diferencia entre usar aov()y lme()analizar datos longitudinales y cómo interpretar los resultados de estos dos métodos? A continuación, se analiza el mismo conjunto de datos usando aov()y lme()y tengo 2 resultados diferentes. Con aov(), obtuve un resultado significativo en la interacción tiempo por tratamiento, pero ajustando …

13 r  mixed-model  repeated-measures  panel-data 

He escuchado comúnmente que los modelos LME son más sólidos en el análisis de datos de precisión (es decir, en experimentos de psicología), ya que pueden trabajar con distribuciones binomiales y otras distribuciones no normales que los enfoques tradicionales (por ejemplo, ANOVA) no pueden. ¿Cuál es la base matemática de …

13 mixed-model 

He estado leyendo sobre el cálculo valores de en modelos mixtos y después de leer las preguntas frecuentes de R-sig, otras publicaciones en este foro (vincularía algunas pero no tengo suficiente reputación) y varias otras referencias entiendo que usando valores en el contexto de modelos mixtos son complicados.R 2R2R2R^2R2R2R^2 Sin …

13 r  mixed-model  r-squared  lme4-nlme 

He estado analizando algunos datos usando el modelado de efectos mixtos lineales en R. Estoy planeando hacer un póster con los resultados y me preguntaba si alguien con experiencia en modelos de efectos mixtos podría sugerir qué gráficos usar para ilustrar los resultados de modelo. Estaba pensando en gráficos residuales, …

13 r  data-visualization  mixed-model 

Me gustaría obtener intervalos de confianza del 95% en las predicciones de un nlmemodelo mixto no lineal . Como no se proporciona nada estándar para hacer esto dentro nlme, me preguntaba si es correcto usar el método de "intervalos de predicción de población", como se describe en el capítulo del …

12 r  mixed-model  confidence-interval  lme4-nlme 

En el modelado multinivel, la correlación intraclase a menudo se calcula a partir de un ANOVA de efectos aleatorios yyo j= γ00+ uj+ eyo jyyoj=γ00+tuj+miyoj y_{ij} = \gamma_{00} + u_j + e_{ij} donde son los residuos de nivel 2 y son los residuos de nivel 1. Luego obtenemos estimaciones, y …

12 correlation  mixed-model  intraclass-correlation 

TL; DR: la lme4optimización parece ser lineal en la cantidad de parámetros del modelo por defecto, y es mucho más lenta que un glmmodelo equivalente con variables ficticias para grupos. ¿Hay algo que pueda hacer para acelerarlo? Estoy tratando de ajustar un modelo logit jerárquico bastante grande (~ 50k filas, …

12 r  mixed-model  optimization  lme4-nlme  numerics 

Resumen de la pregunta Advertencia: esta pregunta requiere mucha configuración. Por favor, tenga paciencia conmigo. Un colega mío y yo estamos trabajando en un diseño de experimento. El diseño debe evitar una gran cantidad de restricciones, que enumeraré a continuación. Desarrollé un diseño que satisface las limitaciones y que nos …

12 mixed-model  experiment-design  simulation  confounding