Manova.fnc
Copy, Paste and Adapts
Manova.fnc(dat, variables=2:6, between.factor=c('zone','gender'), poshoc=c('zone','zone:gender') )
Objetive
Produce la estimación de modelos de Análisis Multivariado de la Varianza (MANOVA).
Manova
Utilizaremos la base de datos ejemplo.manova.Rdata que debes guardar en el mismo directorio donde tengas el toolbox.
stress = read.file.fnc('ejemplo.manova.Rdata')
frequency.fnc(stress, 'escuela:sexo')
Disponemos de cuatro variables dependientes (Y1 hasta Y4) medidas en 60 sujetos con dos factores intergrupo: escuela y sexo (10 sujetos por condición). Queremos estimar un modelo MANOVA del efecto fijo del diseño escuela x sexo sobre combinaciones lineales (variables canónicas) de las 4 medidas diferentes de stress (variables dependientes).
Manova.fnc(stress, variables=3:6,
between.factor=c('escuela','sexo') )
Por defecto la estimación Manova dará como resultado:
Tabla resumen Manova basada en el test de Wilks.
SCPC de todos los efectos y de error.
Funciones discriminantes (canónicas) con coeficientes típicos y estructura para cada efecto fijo del modelo estimado.
Estadísticos descriptivos univariados y Anova univariado para cada variable dependiente.
Gráficas de proyección de centroides en las funciones discriminantes para cada efecto fijo.
La primera figura nos muestra las dos funciones canónicas significativas del efecto del factor escuela. En ella podemos ver que la primera función explica el 64.7% de la varianza conjunta de las variables dependientes frente al 35.5% de la segunda, sobre el 100% de varianza explicada por el factor escuela que es del 1-0.12=88% (1-wilks). La primera función separa a las escuelas medias y elemental de la high, mientras que la segunda separa a la elemental del resto. Asimismo vemos que las variables dependientes Y1 e Y2, tienen la mayor carga en la segunda función, y Y3 e Y4 en la primera.
POSHOC CONTRASTS
Los contrastes poshoc (par a par y ortogonales definidos por el usuario o de tendencia) para cada uno de los efectos se llevan a cabo de la misma forma que lo realizábamos en los poshoc de Anova.
Dado que el factor escuela (con 3 niveles) ha sido significativo, podemos plantearnos contrastes poshoc sobre las dos funciones canónicas significativas.
Manova.fnc(stress, variables=3:6,
between.factor=c('escuela','sexo'),
poshoc='escuela')
Los contrastes poshoc realizados para ambas funciones canónicas confirman lo comentado anteriormente en la gráfica correspondiente al efecto del factor escuela.
En la salida del Manova hemos podido ver que no hay interacción de escuela:sexo (p > 0.05). Sin embargo de haber resultado significativa podríamos solicitar los efectos simples para este efecto mediante:
Manova.fnc(stress, variables=3:6,
between.factor=c('escuela','sexo'),
poshoc='escuela:sexo')
Si el usuario desea que la gráfica de poshoc se presente sexo:escuela (al contrario) debe modificar el orden de dichos factores en between.factor y repetir el procedimiento Manova.