Al especificar un modelo de efectos mixtos cruzados, estoy tratando de incluir interacciones. Sin embargo, recibo el siguiente mensaje de error: Error in lme.formula(rate ~ nozzle, random = ~nozzle | operator, data = Flow) : nlminb problem, convergence error code = 1 message = iteration limit...
Esta pregunta no pertenece específicamente R, pero elegí usarla Rpara ilustrarla. Considere el código para producir bandas de confianza alrededor de una línea qq (normal): library(car) library(MASS) b0<-lm(deaths~.,data=road) qqPlot(b0$resid,pch=16,line="robust") Estoy buscando una...
Acabo de ejecutar un GLM binomial negativo y esta es la salida: Call: glm.nb(formula = small ~ method + site + depth, data = size.dat, init.theta = 1.080668549, link = log) Deviance Residuals: Min 1Q Median 3Q Max -2.2452 -0.9973 -0.3028 0.3864 1.8727 Coefficients: Estimate Std. Error z...
Tengo un conjunto de datos de medidas repetidas desequilibradas para analizar, y he leído que la forma en que la mayoría de los paquetes estadísticos manejan esto con ANOVA (es decir, la suma de cuadrados del tipo III) es incorrecta. Por lo tanto, me gustaría usar un modelo de efectos mixtos para...
Tengo un conjunto de datos en el que la variable que me gustaría usar como efecto aleatorio solo tiene una observación única para algunos niveles. Según las respuestas a las preguntas anteriores, he deducido que, en principio, esto puede estar bien. ¿Puedo ajustar un modelo mixto con sujetos que...
Sorprendentemente, no pude encontrar una respuesta a la siguiente pregunta usando Google: Tengo algunos datos biológicos de varias personas que muestran un comportamiento de crecimiento más o menos sigmoide en el tiempo. Por lo tanto, deseo modelarlo utilizando un crecimiento logístico...
¿Por qué la regresión lineal y el modelo generalizado tienen supuestos inconsistentes? En la regresión lineal, suponemos que el residuo proviene de Gauss En otra regresión (regresión logística, regresión por envenenamiento), asumimos que la respuesta viene de alguna distribución (binomial,...
Actualmente tengo un problema para entender la sintaxis de R para ajustar un GLM usando la distribución Gamma. Tengo un conjunto de datos, donde cada fila contiene 3 covariables ( ), una variable de respuesta ( ) y un parámetro de forma ( ). Quiero modelar la escala de la distribución Gamma como...
Sé cómo realizar una regresión lineal en un conjunto de puntos. Es decir, sé cómo ajustar un polinomio de mi elección, a un conjunto de datos dado (en el sentido LSE). Sin embargo, lo que no sé es cómo forzar a mi solución a pasar por algunos puntos particulares de mi elección. He visto que esto se...
Me preguntaba si la regresión de Poisson tiene un término de error. ¿Puede una regresión de Poisson tener efectos aleatorios y un término de error? Estoy confundido sobre este punto. En la regresión logística, no hay término de error porque su variable de resultado es binaria. ¿Es ese el único...
He estado leyendo algunos documentos, y no estoy claro acerca de las definiciones específicas de Efecto de tratamiento promedio (ATE) y Efecto de tratamiento marginal (MTE). ¿Son lo mismo? De acuerdo con Austin ... Un efecto condicional es el efecto promedio, a nivel de sujeto, de mover un...
Quiero incluir el tiempo dedicado a hacer algo (semanas de lactancia materna, por ejemplo) como una variable independiente en un modelo lineal. Sin embargo, algunas observaciones no se involucran en el comportamiento en absoluto. Codificarlos como 0 no es realmente correcto, porque 0 es...
El tratamiento clásico de la inferencia estadística se basa en el supuesto de que existe una estadística correctamente especificada. Es decir, la distribución P∗(Y)P∗(Y)\mathbb{P}^*(Y) que generó los datos observados yyy es parte del modelo estadístico MM\mathcal{M} :
Estoy estudiando modelos gráficos probabilísticos , un libro para el autoaprendizaje. ¿Los bordes en un gráfico acíclico dirigido (DAG) representan relaciones causales? ¿Qué sucede si quiero construir una red bayesiana , pero no estoy seguro de la dirección de las flechas en ella? Todos los datos...
He estado usando el ajuste del modelo caret, pero luego volví a ejecutar el modelo usando el gbmpaquete. Entiendo que el caretpaquete usa gbmy el resultado debe ser el mismo. Sin embargo, solo una ejecución de prueba rápida usando data(iris)muestra una discrepancia en el modelo de aproximadamente...
Una pregunta probablemente muy básica sobre ANOVA multifactorial. Supongamos un diseño bidireccional donde probamos los efectos principales A, B y la interacción A: B. Cuando se prueba el efecto principal para A con SS tipo I, el efecto SS se calcula como la diferencia R SS( 1 ) - R SS( A...
He estado leyendo sobre el cálculo valores de en modelos mixtos y después de leer las preguntas frecuentes de R-sig, otras publicaciones en este foro (vincularía algunas pero no tengo suficiente reputación) y varias otras referencias entiendo que usando valores en el contexto de modelos mixtos son...
He visto referencias al aprendizaje de la teoría de probabilidad bayesiana en R, y me preguntaba si hay más como esta, ¿quizás específicamente en Python? ¿Orientado hacia el aprendizaje de la teoría de probabilidad bayesiana, inferencia, estimación de máxima verosimilitud, modelos gráficos y el...
Soy relativamente nuevo en estadísticas y agradecería ayuda para entender esto mejor. En mi campo hay un modelo de forma comúnmente utilizado: PAGt= Po( Vt)αPAGt=PAGo(Vt)αP_t = P_o(V_t)^\alpha Cuando las personas ajustan el modelo a los datos, generalmente lo linealizan y ajustan lo...
No tengo ningún fondo en matemáticas, pero entiendo cómo funciona el simple Perceptron y creo que entiendo el concepto de un hiperplano (lo imagino geométricamente como un plano en el espacio 3D que separa dos nubes de puntos, al igual que una línea se separa) dos nubes de puntos en el espacio...