He utilizado la imputación múltiple para obtener una serie de conjuntos de datos completos. He utilizado métodos bayesianos en cada uno de los conjuntos de datos completos para obtener distribuciones posteriores para un parámetro (un efecto aleatorio). ¿Cómo puedo combinar / agrupar los...
Estoy tratando con un modelo lineal jerárquico bayesiano , aquí la red que lo describe. YYY representa las ventas diarias de un producto en un supermercado (observado). XXX es una matriz conocida de regresores, que incluye precios, promociones, día de la semana, clima, días festivos. SSS es el...
Sabemos que una prueba t emparejada es solo un caso especial de ANOVA de medidas repetidas unidireccionales (o dentro del sujeto), así como el modelo lineal de efectos mixtos, que se puede demostrar con la función lme () del paquete nlme en R Como se muestra abajo. #response data from 10 subjects...
¿La refundición de un modelo de regresión lineal multivariante como una regresión lineal múltiple es completamente equivalente? No me refiero a simplemente ejecutar regresiones separadas.ttt He leído esto en algunos lugares (Bayesian Data Analysis - Gelman et al., Y Multivariate Old School -...
Digamos que tengo un conjunto de datos de "catéter renal". Estoy tratando de modelar una curva de supervivencia usando un modelo de Cox. Si considero un modelo de Cox: necesito la estimación del riesgo de referencia. Al usar la función integrada del paquete R , puedo hacerlo fácilmente de esta...
Necesito algunos consejos sobre dos dilemas principales en mi investigación, que es un estudio de caso de 3 grandes productos farmacéuticos e innovación. El número de patentes por año es la variable dependiente. Mis preguntas son ¿Cuáles son los criterios más importantes para un buen modelo?...
Suponga que tiene una sola sección transversal de datos donde los individuos se encuentran dentro de los grupos (por ejemplo, estudiantes dentro de las escuelas) y desea estimar un modelo del formulario Y_i = a + B*X_idonde Xes un vector de características de nivel individual ya una constante. En...
Digamos que tenemos que GLMMs mod1 <- glmer(y ~ x + A + (1|g), data = dat) mod2 <- glmer(y ~ x + B + (1|g), data = dat) Estos modelos no están anidados en el sentido habitual de: a <- glmer(y ~ x + A + (1|g), data = dat) b <- glmer(y ~ x + A + B + (1|g), data = dat) por lo que no...
Si arreglo los valores de los nodos observados de un MRF, ¿se convierte en un
Me interesa aprender cómo calificar y clasificar a las personas en un grupo que solo interactúan / compiten en pares (es decir, sistemas como el sistema de clasificación ELO para el ajedrez). ¿Hay algún método de acceso o métodos más precisos y avanzados? ¿Hay algún paquete R que facilite la...
Estoy tratando de producir un modelo para el que tengo una variable de respuesta que es una proporción entre 0 y 1, esto incluye bastantes 0s y 1s pero también muchos valores intermedios. Estoy pensando en intentar una regresión beta. El paquete que he encontrado para R (betareg) solo permite...
Estoy tratando de usar lmeel nlmepaquete para replicar resultados de aovANOVA de medidas repetidas. He hecho esto para un experimento de medidas repetidas de un solo factor y para un experimento de dos factores con un factor entre sujetos y un factor dentro de los sujetos, pero tengo problemas para...
Supongamos que estamos trabajando en un modelo de efectos aleatorios de algunos datos de conteo a lo largo del tiempo, y queremos controlar algunas tendencias. Normalmente, harías algo como: lmer(counts ~ dependent_variable + (1+t+I(t^2)|ID), family="poisson") para incluir una forma cuadrática...
Tengo una hoja de cálculo de algunos puntajes del equipo. El primer equipo en ganar 10 puntos. Hay 2 jugadores en cada equipo. Los jugadores juegan con diferentes compañeros de equipo todo el tiempo, aunque no se eligen perfectamente al azar. No se mantienen puntajes individuales. Básicamente,...
Básicamente, me pregunto cómo se aplican las diferentes estructuras de covarianza y cómo se calculan los valores dentro de estas matrices. Funciones como lme () nos permiten elegir qué estructura nos gustaría, pero me encantaría saber cómo se estiman. Considere el modelo de efectos lineales mixtos...
Hay una distinción que me hace tropezar con modelos mixtos, y me pregunto si podría obtener algo de claridad. Supongamos que tiene un modelo mixto de datos de conteo. Hay una variable que sabe que quiere como efecto fijo (A) y otra variable para el tiempo (T), agrupada por decir una variable de...
Por favor considere estos datos: dt.m <- structure(list(id = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12), occasion = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c("g1", "g2"), class =...
En respuesta a esta pregunta, con respecto a si mi diseño en el que presentaba al azar a los participantes con imágenes de diferentes categorías era un ejemplo en el que debería usar un ANOVA de medidas repetidas, obtuve la respuesta de que debería usar un modelo mixto en su lugar, con uno de los...
Esto parece muy elemental, pero siempre me atoro en este punto ... La mayoría de los datos con los que trato no son normales, y la mayoría de los análisis se basan en una estructura GLM. Para mi análisis actual, tengo una variable de respuesta que es "velocidad de marcha" (metros / minuto). ¡Es...
Tengo una gráfica de valores residuales de un modelo lineal en función de los valores ajustados donde la heterocedasticidad es muy clara. Sin embargo, no estoy seguro de cómo proceder ahora porque, por lo que entiendo, esta heterocedasticidad invalida mi modelo lineal. (¿Está bien?) Use un ajuste...