a.- Anova multifactorial intergrupo

A partir de la base de datos OBrienKaiser declaramos los 2 factores intergrupo del experimento con el objeto fac.inter. Por comodidad (nombre mas corto) asignamos a datos dicha base de datos.

datos=OBrienKaiser; head(datos)

Observa que los factores intergrupo indicados en fac.inter, deben coincidir con sus respectivos nombres en la base de datos.

fac.inter= c('treatment','gender')

Anova.fnc(datos, fac.inter=fac.inter, vd='pre.1')

Los argumentos obligatorios para este ejemplo serán obviamente el nombre de la base de datos (datos) , los nombres de los factores intergrupo (fac.inter) y la variable numérica (vd) sobre la que deseamos observar los efectos fijos (pre.1).

Por defecto para este tipo de Anova la función te dará:

1. Tabla resumen de ANOVA

2. Prueba de Levene

3. Tamaño del efecto

4. Potencia observada

5. Estadísticos descriptivos

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

Veamos otros argumentos posibles:

Por defecto las gráficas con los intervalos de confianza se crean en tonos de grises (argumento color=F). Si incluyes el argumento color=T, creará las gráficas en color.

Anova.fnc(datos, fac.inter=fac.inter, vd='pre.1', color=T)

En el supuesto de que los colores no fueran de tu agrado, puedes hacer uso del argumento que.colores con los códigos en versión hexadecimal que puedes buscar en http://www.colorhexa.com. Veamos un ejemplo en el que deseamos que los tres colores de la gráfica izquierda sean en tres tonos diferentes de azul (2 figuras inferiores). Nos trasladamos a la página web de colorhexa y buscamos blue en el cuadro de búsqueda del principio de la página. Descendemos hasta Color Schemes y pasamos el cursor del ratón por encima de los 3 tonos que deseamos. De esa manera podremos ver los 3 códigos correspondientes: #0080ff, #0000ff y #8000ff. Ahora debemos incluir esos tres códigos (entre apóstrofes) en el argumento que.colores y volver a ejecutar la función Anova.fnc.

Anova.fnc(datos, fac.inter=fac.inter, vd='pre.1',

que.colores=c('#0080ff','#0000ff' ,'#8000ff'))

Si queremos que esas gráficas sean guardadas en un archivo pdf, debemos incluir el argumento to.pdf=T. Se guardará un archivo pdf en tu directorio de trabajo de nombre graf_anova.pdf con todas las gráficas de barras de los efectos estimados.

Anova.fnc(datos, fac.inter=fac.inter, vd='pre.1', to.pdf=T)

Si queremos que esas gráficas sean guardadas una a una en un archivo tiff (sin compresión) para poder editar posteriormente en aplicaciones de gráficos vectoriales y por capas como gimp, etc. Incluiremos el argumento to.tif=T. Se guardará en tu carpeta de trabajo tantas gráficas de barras como efectos principales e interacciones (sólo dobles) tengas en la tabla resumen.

Anova.fnc(datos, fac.inter=fac.inter, vd='pre.1', to.tif=T)