Hace poco comencé a aprender sobre los modelos mixtos lineales generalizados y estaba usando R para explorar la diferencia que hace al tratar la pertenencia a grupos como un efecto fijo o aleatorio. En particular, estoy mirando el conjunto de datos de ejemplo discutido
¿Debo incluir efectos aleatorios en un modelo incluso si no son estadísticamente significativos? Tengo un diseño experimental de medidas repetidas, en el que cada individuo experimenta tres tratamientos diferentes en orden aleatorio. Me gustaría controlar cualquier efecto individual y de orden,...
Espero que alguien pueda ayudar con lo que creo que es una pregunta relativamente simple, y creo que sé la respuesta, pero sin confirmación se ha convertido en algo de lo que no puedo estar seguro. Tengo algunos datos de conteo como una variable de respuesta y quiero medir cómo cambia esa...
Tengo algunas preguntas con respecto a la especificación e interpretación de GLMM. 3 preguntas son definitivamente estadísticas y 2 son más específicamente sobre R. Estoy publicando aquí porque, en última instancia, creo que el problema es la interpretación de los resultados de GLMM. Actualmente...
Introducción Tengo participantes que están tocando repetidamente superficies contaminadas con E. coli en dos condiciones ( A = usando guantes, B = sin guantes). Quiero saber si hay una diferencia entre la cantidad de bacterias en la punta de los dedos con y sin guantes, pero también entre la...
Primero simulemos algunos datos para una regresión logística con partes fijas y aleatorias: set.seed(1) n <- 100 x <- runif(n) z <- sample(c(0,1), n, replace=TRUE) b <- rnorm(2) beta <- c(0.4, 0.8) X <- model.matrix(~x) Z <- cbind(z, 1-z) eta <- X%*%beta + Z%*%b pr <-...
Mi diseño es el siguiente. yyy es la respuesta de Bernoulli X1X1x_1 es una variable continua X2X2x_2 es una variable categórica (factor) con dos niveles El experimento está completamente dentro de los sujetos. Es decir, cada sujeto recibe cada combinación deX1X1x_1 y X2X2x_2. Esta es una...
He leído los excelentes comentarios sobre cómo lidiar con los valores perdidos antes de aplicar SVD, pero me gustaría saber cómo funciona con un ejemplo simple: Movie1 Movie2 Movie3 User1 5 4 User2 2 5 5 User3 3 4 User4 1 5 User5 5 1 5 Dada la matriz anterior, si elimino los valores de NA,...
Tengo datos de recolección a largo plazo y me gustaría probar si el número de animales recolectados está influenciado por los efectos del clima. Mi modelo se ve a continuación: glmer(SumOfCatch ~ I(pc.act.1^2) +I(pc.act.2^2) + I(pc.may.1^2) + I(pc.may.2^2) + SampSize + as.factor(samp.prog) +...
Considere el siguiente ejemplo sencillo: library( rms ) library( lme4 ) params <- structure(list(Ns = c(181L, 191L, 147L, 190L, 243L, 164L, 83L, 383L, 134L, 238L, 528L, 288L, 214L, 502L, 307L, 302L, 199L, 156L, 183L), means = c(0.09, 0.05, 0.03, 0.06, 0.07, 0.07, 0.1, 0.1, 0.06, 0.11,...
Estoy usando el lme4paquete en R para hacer algunos modelos logísticos de efectos mixtos. Entiendo que la suma de cada efecto aleatorio debe ser cero. Cuando hago modelos mixtos lineales de juguete usando lmer, los efectos aleatorios generalmente son < confirmando mi creencia de que el But en...
Me está costando llegar a la misma página que mi supervisor cuando se trata de validar mi modelo. He analizado los residuos (observados contra los valores ajustados) y utilicé esto como argumento para discutir los resultados obtenidos por mi modelo, sin embargo, mi supervisor insiste en que la...