Antes de hacer esta pregunta, busqué en nuestro sitio y encontré muchas preguntas similares (como aquí , aquí y aquí ). Pero creo que esas preguntas relacionadas no fueron bien respondidas o discutidas, por lo que me gustaría plantear esta pregunta nuevamente. Siento que debería haber una gran...
Tengo el siguiente modelo lineal: Para abordar la heteroscedasticidad residual, he tratado de aplicar una transformación logarítmica en la variable dependiente como Iniciar sesión( Y+ 1 )Iniciar sesión(Y+1)\log(Y + 1) pero todavía veo el mismo efecto de abanico en los residuos. Los valores de...
Tengo una serie de distribuciones de cola izquierda / sesgada que me gustaría mostrar. Hay 42 distribuciones a través de tres factores (etiquetados como A, By Cmás adelante). Además, la variación se está reduciendo a través del factor B. El problema que tengo es que las distribuciones son...
Tengo datos del siguiente diseño experimental: mis observaciones son recuentos del número de éxitos ( K) del número correspondiente de ensayos ( N), medidos para dos grupos, cada uno compuesto de Iindividuos, de Ttratamientos, donde en cada combinación de factores hay Rréplicas . Por lo tanto, en...
Una de las suposiciones de la regresión lineal es que debe haber una variación constante en los términos de error y que los intervalos de confianza y las pruebas de hipótesis asociadas con el modelo se basan en esta suposición. ¿Qué sucede exactamente cuando los términos de error no tienen una...
Supongamos que observo vectores variables independientes y y la variable dependiente . Me gustaría ajustar un modelo de la forma: donde es una función dos veces diferenciable de valor positivo, \ sigma es un parámetro de escala desconocido y \ epsilon es una variable aleatoria gaussiana de media...
Usando R en algunos datos e intentando ver si mis datos son o no heteroscedásticos, he encontrado dos implementaciones de la prueba Breusch-Pagan, bptest (package lmtest) y ncvTest (package car). Sin embargo, estos producen resultados diferentes. ¿Cuál es la diferencia entre los dos? ¿Cuándo...
Estoy tratando de simular un conjunto de datos que coincida con los datos empíricos que tengo, pero no estoy seguro de cómo estimar los errores en los datos originales. Los datos empíricos incluyen la heterocedasticidad, pero no estoy interesado en transformarlos, sino en utilizar un modelo lineal...
En este artículo , ( Inferencia bayesiana para componentes de varianza usando solo contrastes de error , Harville, 1974), el autor afirma es un "bien conocido relación ", para una regresión lineal
Quiero ejecutar correlaciones en una serie de mediciones donde se usaron escalas Likert. Al observar los diagramas de dispersión, parece que los supuestos de linealidad y homocedasticidad pueden haber sido violados. Dado que parece haber cierto debate sobre la calificación del nivel ordinal que...
Sé que OLS es imparcial pero no eficiente bajo heteroscedasticidad en una configuración de regresión lineal. En Wikipedia http://en.wikipedia.org/wiki/Minimum_mean_square_error El estimador MMSE es asintóticamente imparcial y converge en distribución a la distribución normal: , donde I (x) es...
En el método de mínimos cuadrados queremos estimar los parámetros desconocidos en el modelo: Yj=α+βxj+εj(j=1...n)Yj=α+βxj+εj(j=1...n)Y_j = \alpha + \beta x_j + \varepsilon_j \enspace (j=1...n) Una vez que hayamos hecho eso (para algunos valores observados), obtenemos la línea de regresión...
En Rpuedo realizar una prueba de Breusch-Pagan para la heterocedasticidad utilizando la ncvTestfunción del carpaquete. Una prueba de Breusch-Pagan es un tipo de prueba de ji cuadrado. ¿Cómo interpreto estos resultados? > require(car) > set.seed(100) > x1 = runif(100, -1, 1) > x2 =...
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...
Tengo problemas para distinguir entre los conceptos de escedasticidad y estacionariedad. Según tengo entendido, la heterocedasticidad es una variabilidad diferente en las subpoblaciones y la no estacionariedad es una media / varianza cambiante a lo largo del tiempo. Si esta es una comprensión...
Tengo un conjunto de datos de ejemplo de la siguiente manera: Volume <- seq(1,20,0.1) var1 <- 100 x2 <- 1000000 x3 <- 30 x4 = sqrt(x2/pi) H = x3 - Volume r = (x4*H)/(H + Volume) Power = (var1*x2)/(100*(pi*Volume/3)*(x4*x4 + x4*r + r*r)) Power <- jitter(Power, factor = 1, amount...
Estoy buscando ayuda, consejos o sugerencias sobre cómo explicar la heterogeneidad / heterocedasticidad a los biólogos de mi departamento. En particular, quiero explicar por qué es importante buscarlo y tratarlo si existe, estaba buscando opiniones sobre las siguientes preguntas. ¿La...
Estoy buscando la terminología estadística correcta para describir el siguiente problema. Quiero caracterizar un dispositivo electrónico que tenga una respuesta lineal Y=β0 0+β1X+ ϵY=β0+β1X+ϵY = \beta_0 + \beta_1 X + \epsilon dónde ϵ ∼ N( 0 ,σ2r o)ϵ∼N(0,σro2)\epsilon \sim N(0,\sigma^2_{ro})es un...
En estos días estoy trabajando con Breusch-Pagan para probar la homocedasticidad. He probado los precios de dos acciones con este método. Este es el resultado: > mod <- lm(prices[,1] ~ prices[,2]) > bp <- bptest(mod) > bp studentized Breusch-Pagan test data: prices[, 1] ~...
Estoy tratando de crear un modelo de predicción usando la regresión. Este es el diagrama de diagnóstico para el modelo que obtengo al usar lm () en R: Lo que leí en el gráfico QQ es que los residuos tienen una distribución de cola pesada, y el gráfico Residuals vs Fitted parece sugerir que la...