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...
La ocurrencia no tan infrecuente cuando se trata de modelos mixtos máximos complejos (estimar todos los efectos aleatorios posibles para datos y modelo dados) es perfecta (+1 o -1) o una correlación casi perfecta entre algunos efectos aleatorios. Para el propósito de la discusión, observemos el...
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...
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 |...
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
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 tratando de lidiar con un análisis de tiempo hasta el evento utilizando resultados binarios repetidos. Supongamos que el tiempo hasta el evento se mide en días, pero por el momento discretizamos el tiempo en semanas. Quiero aproximar un estimador de Kaplan-Meier (pero permitir covariables)...
Los siguientes injertos se toman de este artículo . Soy novato en bootstrap e intento implementar el bootstrapping paramétrico, semiparamétrico y no paramétrico para el modelo mixto lineal con R bootpaquete. Código R Aquí está mi
Estoy tratando de analizar el efecto del año en el registro variable para un grupo particular de individuos (tengo 3 grupos). El modelo más simple: > fix1 = lm(logInd ~ 0 + Group + Year:Group, data = mydata) > summary(fix1) Call: lm(formula = logInd ~ 0 + Group + Year:Group, data =...
Estoy tratando de producir un modelo mixto lineal, el código R es el siguiente. lme (Promedio.pago ~ Juego + Tipo + Otros.Tipo + Juego: Tipo + Juego: Otros.Tipo + Tipo: Otros.Tipo, aleatorio = ~ 1 | Sujetos, método = "REML", datos = Sujetosm1) -> lme1 El término de respuesta Average.payoff...
En un estudio longitudinal, los resultados de las unidades i se miden repetidamente en los puntos de tiempo t con un total de m ocasiones de medición fija (fija = las mediciones en unidades se toman al mismo tiempo).YitYitY_{it}iiitttmmm Las unidades se asignan aleatoriamente a un tratamiento, , o...
Tengo una pregunta básica sobre los enfoques para promediar modelos usando criterios de TI para ponderar modelos dentro de un conjunto de candidatos. La mayoría de las fuentes que he leído sobre promedios de modelos abogan por promediar las estimaciones de coeficientes de parámetros basadas en...
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
¿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,...
Ejemplos: Tengo una oración en la descripción del trabajo: "Ingeniero senior de Java en el Reino Unido". Quiero usar un modelo de aprendizaje profundo para predecirlo en 2 categorías: English y IT jobs. Si uso el modelo de clasificación tradicional, solo puede predecir 1 etiqueta con...
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...
La distinción entre variables estrictamente aleatorias (que deberían modelarse como tales) y variables no aleatorias que algunos argumentan que podrían modelarse como aleatorias si se trata de un modelo jerárquico / multinivel, me resulta borrosa. Bates y Bolker ejemplifican los efectos aleatorios...
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 <-...