He leído que la distribución para el estimador de efectos aleatorios en lme4 está muy sesgada y por esta razón no se informan los errores estándar. Me pregunto si alguien puede proporcionar una referencia para esto. Tengo acceso al libro de Bates y Pinherio, pero no de Raudenbush y Bryk. Un...
¿Existe un método de Box-Cox para modelos gaussianos lineales con efectos
Tengo un conjunto de datos muy grande donde he repetido mediciones a lo largo del tiempo para ubicaciones individuales. Algunas ubicaciones pueden tener 10 puntos de datos y algunas ubicaciones tienen solo 1 punto de datos. Encajo un modelo mixto y uso ubicaciones como efectos aleatorios. Mi...
Tengo un conjunto de números que se supone que provienen de una distribución de Poisson. El conjunto también tiene algunos valores atípicos y, debido a eso, las estimaciones de máxima probabilidad se ven gravemente afectadas. Escuché que los procedimientos de estimación robustos pueden ayudar en...
La regresión logística binomial tiene asíntotas superior e inferior de 1 y 0 respectivamente. Sin embargo, los datos de precisión (solo como ejemplo) pueden tener asíntotas superiores e inferiores muy diferentes a 1 y / o 0. Puedo ver tres posibles soluciones a esto: No se preocupe si obtiene...
Supongamos que A es un factor fijo y B es un factor aleatorio. A y B están cruzados. Otro factor aleatorio C está anidado en B. Teniendo en cuenta las interacciones AB y AC, ¿cómo debo ajustar el modelo usando la función lme en el paquete nlme?
Estoy tratando de obtener predicciones para observaciones de un objeto lme. Se supone que esto es bastante sencillo. Sin embargo, dado que recibo diferentes tipos de errores para diferentes ensayos, me parece que me falta algo. Mi modelo es el siguiente: model <- lme(log(child_mortality) ~...
Me gustaría construir predicciones para un modelo mixto (logística a través de glmer) en un nuevo conjunto de datos usando solo los efectos fijos, manteniendo los efectos aleatorios en 0. Pero tengo problemas para configurar la matriz del modelo para poder calcularlos. Dado que la clase mer no...
Me gustaría obtener las simulaciones posteriores de los componentes de varianza de un modelo lmer () con la función mcmcsamp (). Cómo hacer ? Por ejemplo, a continuación se muestra el resultado de un ajuste lmer (): > fit Linear mixed model fit by REML Formula: y ~ 1 + (1 | Part) + (1 |...
Espero que alguien pueda ayudarme con qué modelo (fragilidad, estratos o clúster) debería usar para mis datos. He emparejado datos, así que tengo que tenerlos en cuenta al modelar el Cox PH y no estoy seguro de qué modelo me dará un resultado más preciso. Mi estudio analizó el tiempo que le tomó a...
Bien, este es mi problema: Tengo 12 participantes. Cada participante pasó 3 noches en mi laboratorio haciendo una tarea de tiempo de reacción en cuatro puntos de tiempo durante la noche (12, 1, 2 y 3 en punto), con una semana entre cada una de esas noches. Cada noche, cada participante estuvo...
Estoy calculando algunas probabilidades condicionales y los intervalos de confianza del 95% asociados. Para muchos de mis casos, tengo recuentos directos de xéxitos fuera de los nensayos (de una tabla de contingencia), por lo que puedo usar un intervalo de confianza binomial, como se proporciona...
Digamos, por ejemplo, que tenemos un modelo de regresión logística que genera la probabilidad de que un paciente desarrolle una enfermedad particular basada en muchas covariables. Podemos tener una idea de la magnitud y dirección del efecto de cada covariable en general examinando los...
Trabajo para una organización que supervisa alrededor de 180 escuelas en todo el país. Regularmente recopilamos datos de resultados de estas escuelas y los clasificamos. Me han pedido que haga una encuesta en la oficina nacional (alrededor de 100 personas) para determinar cuáles son las mejores y...
Tengo un problema de regresión lineal. En resumen, tengo un conjunto de datos, lo dividí en dos subconjuntos. Un subconjunto se usa para encontrar la regresión lineal (subconjunto de entrenamiento), otro se usa para evaluarlo (subconjunto de evaluación). Mi pregunta es ¿cómo evaluar el resultado de...
Recientemente completé un estudio mediante el cual asigné al azar a los participantes a uno de los dos grupos de tratamiento. Probé a los participantes al inicio del estudio, inmediatamente después de la intervención, 1 mes y 4 meses en un número bastante grande de variables de resultado. Estaba...
Tengo dos grupos de 10 participantes que fueron evaluados tres veces durante un experimento. Para probar las diferencias entre los grupos y entre las tres evaluaciones, ejecuté un ANOVA de diseño mixto 2x3 con group(control, experimental), time(primero, segundo, tres) y group x time. Ambos timey...
Cerrado. Esta pregunta está fuera de tema . Actualmente no está aceptando respuestas. ¿Quieres mejorar esta pregunta? Actualice la pregunta para que esté en el tema de Cross Validated. Cerrado hace 8 meses . Actualmente estoy en el proceso de...
He estado investigando el paquete de ratones , y aún no he descubierto una forma de usar las imputaciones múltiples para hacer un modelo de Cox, luego validar ese modelo con la función del paquete rmsvalidate() . Aquí hay un código de muestra de lo que tengo hasta ahora, usando el conjunto de datos...
Estoy tratando de realizar un análisis de crecimiento de clase latente (LCGA) y / o modelos de mezcla de crecimiento (GMM) en R. Los datos que estoy usando es un número cada vez mayor de tenedores de repositorios git (variable discreta, no categórica), como puedes ver en este conjunto de datos...