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)

Anova medidas repetidas

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á:

  1. Tabla resumen de ANOVA asumiendo esfericidad

  2. Prueba de esfericidad de Mauchly

  3. Corrección por no esfericidad de Greenhose-Geisser y HF

  4. Tamaño del efecto con eta cuadrado parcial y semiparcial

  5. Potencia observada

  6. Estadísticos descriptivos

  7. Gráfica de barras con los intervalos de confianza según Masson y Loftus(2003) para cada efecto fijo estimado.