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).
Salida Anova intergrupo
Por defecto para este tipo de Anova la función te dará:
Tabla resumen de ANOVA
Prueba de Levene
Tamaño del efecto
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.
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)