Estoy haciendo un estudio de simulación que requiere estimaciones de arranque obtenidas de un modelo mixto lineal generalizado (en realidad, el producto de dos estimaciones para efectos fijos, una de un GLMM y otra de un LMM). Hacer bien el estudio requeriría alrededor de 1000 simulaciones con 1000...
He buscado en Google y este sitio y todavía estoy confundido acerca de la función lmer en la biblioteca lme4. Tengo algunos datos recopilados de diferentes salas psiquiátricas, que tienen una estructura multinivel. Para simplificar, elegiré dos variables de nivel 2 y dos de nivel 1, aunque en...
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...
Tengo dos modelos de regresión logística en R hechos con glm(). Ambos usan las mismas variables, pero se hicieron usando diferentes subconjuntos de una matriz. ¿Hay una manera fácil de obtener un modelo promedio que proporcione las medias de los coeficientes y luego usar esto con la función predic...
Estoy investigando la interacción entre dos variables ( y x 2 ). Existe una gran correlación lineal entre estas variables con r > 0.9 . Por la naturaleza del problema, no puedo decir nada sobre la causalidad (si x 1 causa x 2 o al revés). Me gustaría estudiar las desviaciones de la línea de...
Si hay 0 en la tabla de contingencia y estamos ajustando modelos poisson / loglineales anidados (usando la glmfunción de R ) para una prueba de razón de probabilidad, ¿necesitamos ajustar los datos antes de ajustar los modelos glm (por ejemplo, agregar 1/2 a todos los recuentos)? Obviamente,...
En R, ¿por qué la configuración predeterminada de qqplot(linear model)uso de los residuos estandarizados en el eje y? ¿Por qué R no usa los residuos
Me han pedido que analice algunos datos de un ensayo clínico que busca dos métodos para medir la presión arterial. Tengo datos de 50 sujetos, cada uno con entre 2 y 57 medidas utilizando cada método. Me pregunto cómo proceder mejor. Obviamente, necesito una solución que tenga en cuenta el hecho...
Suponga que, en un estudio de 15 sujetos, la variable de respuesta (res) se modela con dos variables explicativas, una (nivel) es categórica con 5 niveles y la otra (tiempo de respuesta: RT) es continua. Con lmer en el paquete lme4 de R, tengo: fm1 <- lmer(res ~ level * RT + (level-1 |...
Dada la siguiente configuración experimental: Se toman múltiples muestras de un sujeto y cada muestra se trata de múltiples maneras (incluido un tratamiento de control). Lo que es principalmente interesante es la diferencia entre el control y cada tratamiento. Puedo pensar en dos modelos simples...
Estas preguntas pueden sonar estúpidas, pero ... ¿ es correcto que los efectos aleatorios solo puedan aplicarse a variables categóricas (como identificación individual, identificación de población, ...), por ejemplo, digamos que es una variable categórica:xixix_i yiyiy_i
Wikipedia sugiere que una forma de ver la confiabilidad entre evaluadores es usar un modelo de efectos aleatorios para calcular la correlación intraclase . El ejemplo de correlación intraclase habla de mirar σ2ασ2α+ σ2ϵσα2σα2+σϵ2\frac{\sigma_\alpha^2}{\sigma_\alpha^2+\sigma_\epsilon^2} de un...
Tengo el siguiente tipo de datos. He evaluado a 10 personas cada una repetidas 10 veces. Tengo una matriz de relación de 10x10 (relación entre todas las combinaciones de los individuos). set.seed(1234) mydata <- data.frame (gen = factor(rep(1:10, each = 10)), repl = factor(rep(1:10, 10)), yld...
Mi colega y yo estamos ajustando una gama de modelos de efectos mixtos lineales y no lineales en R. Se nos pide que realicemos una validación cruzada en los modelos ajustados para poder verificar que los efectos observados son relativamente generalizables. Esta es normalmente una tarea trivial,...
En R, ¿cómo especifico el modelo lmer sin efecto fijo global? Por ejemplo, si digo algo como lmer(y ~ (1 | group) + (0 + x | group), data = my_df) el modelo ajustado será yij=a+αi+βixijyij=a+αi+βixijy_{ij} = a + \alpha_{i} + \beta_{i} x_{ij} ¿Cómo encajo modelo? yij=αi+βixijyij=αi+βixijy_{ij}...
Estoy mirando si la abundancia está relacionada con el tamaño. El tamaño es (por supuesto) continuo, sin embargo, la abundancia se registra en una escala tal que A = 0-10 B = 11-25 C = 26-50 D = 51-100 E = 101-250 F = 251-500 G = 501-1000 H = 1001-2500 I = 2501-5000 J = 5001-10,000 etc... A a Q...
Digamos que genero la probabilidad de un resultado basado en cierto factor y trazo la curva de ese resultado. ¿Hay alguna forma de extraer la ecuación para esa curva de R? > mod = glm(winner~our_bid, data=mydat, family=binomial(link="logit")) > summary(mod) Call: glm(formula = winner ~...
¿Cuál es la diferencia entre un modelo GLM (regresión logística) con una variable de respuesta binaria que incluye sujeto y tiempo como covariables y el modelo GEE análogo que tiene en cuenta la correlación entre mediciones en múltiples puntos de tiempo? Mi GLM se ve así: Y(binary) ~ A +...
He creado mi propia versión ligeramente mejorada del termplot que uso en este ejemplo, puedes encontrarla aquí . He publicado anteriormente en SO, pero cuanto más lo pienso, creo que esto probablemente esté más relacionado con la interpretación del modelo de riesgos proporcionales de Cox que con la...
¿Cómo usaría el resultado de un análisis de componentes principales (PCA) en un modelo lineal generalizado (GLM), suponiendo que el PCA se use para la selección de variables para el GLM? Aclaración: Quiero usar PCA para evitar el uso de variables correlacionadas en el GLM. Sin embargo, PCA me da...