¿Por qué es negativa la variación de mi efecto aleatorio?

Question

Tengo dos variables, response y group, y he ajustado el modelo en el que group es un efecto aleatorio.

proc mixed data=myData nobound;
   class group;
   model response =;
   random group;
run;

Sin embargo, SAS devuelve una desviación negativa para el efecto aleatorio...

Me he dado cuenta de que, inesperadamente, la variabilidad de la response dentro group puede ser mayor que la variabilidad total de response. Supongo que esta es la razón por la SAS devuelve una desviación negativa. Es su opinión? Si es así, ¿cómo podría ser esto visto desde un punto de vista matemático?

DATOS

group   response
1   70.8
1   64.2
1   67.8
1   62
1   60 
1   66.6
1   67.8
1   60
1   70
1   63.3
1   37
1   53.8
2   82.7
2   65.5
2   67.7
2   51.6
2   53.8
2   64.3
2   72.4
2   48.3
2   54.8
2   64.8
3   88.5
3   69
3   60.9
3   71.4
3   52
3   58.3
3   63.6
3   51.9
3   52.4
3   52.2
3   65.4 

Answer 1

Cuando se utiliza ANOVA o MINQUE método de estimación de componentes de varianza de la estimación de la varianza de un factor aleatorio de hecho puede que en ocasiones sea negativo. Puede haber varias causas, y no es fácil detectar la real. Ver Hocking R. R. (1985), El análisis de modelos lineales. Por ejemplo, los valores atípicos y el pequeño tamaño de la muestra puede ser culpable. Trate de eliminar un par de valores extremos y/o ampliar la muestra y ver lo que se obtiene.

Además. Algunas de las posibles razones:

• Variación de obsevations puede ser demasiado grande para el tamaño de la muestra. Más la muestra que se necesita.
• Los valores atípicos
• La verdadera varianza del factor aleatorio es 0
• Grupos en el azar o el factor fijo es demasiado desequilibrada (desiguales de tamaño)
• Componentes de varianza no es un modelo adecuado para la estructura de covarianza de los datos

Answer 2

Cuándo usted borrar nobound de código de SAS la varianza del grupo = 0. También al utilizar el código siguiente puede obtener los mismos resultados; PROC VARCOMP MÉTODO = REML; Grupo de clases; MODELO respuesta = grupo/fijo = 0; ejecutar; Creo que estos resultados confirman que la varianza del grupo = 0 a pesar de algunos métodos dieron a valor negativo como Mivque.

Answer 3

Cuando duplicó sus datos y realizar el análisis se obtendrá también valor negativo pero menos. Para cada aumento de datos tamaño de la disminución de la varianza negativa. El segundo factor es que las diferencias entre los tres medios de los grupos estaba cerca. Estos dos factores están responsables del valor negativo.