c.- Anova multifactorial Split-plot
La base de datos OBrienKaiser se genera a partir de un diseño split-plot 3 x 2 x 3 x 5 tratamiento x genero x fase x hora. Los factores tratamiento y genero son intergrupo los cuales declararemos con el objeto fac.inter y fase y hora los dos factores de medidas repetidas (intragrupo) que definiremos con la lista fac.intra.
Por comodidad (nombre mas corto) asignamos a datos dicha base de datos.
datos=OBrienKaiser; head(datos)
Declaramos el diseño mediante el vector fac.inter y la lista fac.intra.
fac.inter= c('treatment','gender')
fac.intra= list(fase = c('antes','despues','seguimiento'), hora= 1:5 )
Observa que tanto los nombres como los niveles de los factores de medidas repetidas son completamente arbitrarios. A diferencia de los factores intergrupo que deben ser denominados con el mismo nombre que tienen en la base de datos que se esta analizando.
Anova.fnc(datos, fac.inter=fac.inter, fac.intra=fac.intra,
col.empieza.mr=3, ylim=c(0,12), tipo=2)
Como tenemos medidas repetidas, la función necesita obligatoriamente saber en que columna empieza las 3 x 5 condiciones de medidas repetidas del experimento. Si observas la figura anterior veras que estas lo hacen en la columna 3 de la base de datos. Lo indicamos con el argumento: col.empieza.mr=3.
Dado que el número de sujetos por niveles de los factores intergrupo no es igual solicitamos una estimación del modelo Anova por sumas de cuadrados tipo II (tipo=2)
El argumento graficas actúa por defecto como verdadero (graficas=T). Como tenemos mas de tres factores y ello daría lugar a un gran número de gráficas la función te avisa al final de los resultados de que dichas gráficas se guardarán en un archivo externo.
*** Warning. Con mas de 3 factores, las graficas se guardan en pdf
*** Se ha creado el archivo pdf graficas_anova.pdf en el directorio activo
*** Si deseas otro nombre para el archivo incluye el argumento nombre ***
*** Ej: nombre='stress1' y se creara el archivo graf_anova_stress1.pdf ***
Ver gráficas IC anova.split.plot
Como puedes ver en el resultado anterior, por defecto el Anova split-plot te dará:
Tabla resumen de ANOVA asumiendo esfericidad
Prueba de esfericidad de Mauchly
Corrección por no esfericidad de Greenhose-Geisser y HF
Tamaño del efecto con eta cuadrado parcial y semiparcial
Potencia observada
Estadísticos descriptivos
Gráfica de barras con los intervalos de confianza según Masson y Loftus(2003) para cada efecto fijo estimado.