Esta página aún no se ha traducido para esta versión. Puede ver la versión más reciente de esta página en inglés.

anova2

El análisis bidireccional de la varianza

Descripción

realiza un análisis de varianza de dos vías (ANOVA) con diseños equilibrados.anova2 Para realizar ANOVA de dos vías con diseños no balanceados, vea.anovan

ejemplo

p = anova2(y,reps) Devuelve los-valores de un ANOVA de dos vías equilibrado para comparar los medios de dos o más columnas y dos o más filas de las observaciones.py

es el número de réplicas para cada combinación de grupos de factores, que deben ser constantes, indicando un diseño equilibrado.reps Para diseños desequilibrados, use.anovan La función prueba los efectos principales de los factores de columna y fila y su efecto de interacción.anova2 Para probar el efecto de interacción, debe ser mayor que 1.reps

también muestra la tabla ANOVA estándar.anova2

ejemplo

p = anova2(y,reps,displayopt) habilita la visualización de la tabla ANOVA cuando es (por defecto) y suprime la visualización cuando es.displayopt'on'displayopt'off'

ejemplo

[p,tbl] = anova2(___) Devuelve la tabla ANOVA (incluidas las etiquetas de columna y fila) en la matriz de celdas.tbl Para copiar una versión de texto de la tabla ANOVA en el portapapeles, seleccione el menú.Edit > Copy Text

ejemplo

[p,tbl,stats] = anova2(___) Devuelve una estructura, que puede usar para realizar una.statsprueba de comparación múltiple Una prueba de comparación múltiple le permite determinar qué pares de medios de grupo son significativamente diferentes. Para realizar esta prueba, utilice, proporcionando la estructura como entrada.multcomparestats

Ejemplos

contraer todo

Cargue los datos de ejemplo.

load popcorn popcorn
popcorn = 6×3

    5.5000    4.5000    3.5000
    5.5000    4.5000    4.0000
    6.0000    4.0000    3.0000
    6.5000    5.0000    4.0000
    7.0000    5.5000    5.0000
    7.0000    5.0000    4.5000

Los datos provienen de un estudio de marcas de palomitas de maíz y tipos de Popper (Hogg 1987). Las columnas de la matriz son marcas, gourmet, nacional y genérica, respectivamente.popcorn Las hileras son tipos de Popper, aceite y aire. En el estudio, los investigadores reventaron un lote de cada marca tres veces con cada Popper, es decir, el número de replicaciones es 3. Las tres primeras filas corresponden a la Popper de aceite, y las tres últimas filas corresponden a la Popper de aire. Los valores de respuesta son el rendimiento en tazas de palomitas de maíz reventado.

Realice un ANOVA de dos vías. Guarde la tabla ANOVA en la matriz de celdas para facilitar el acceso a los resultados.tbl

[p,tbl] = anova2(popcorn,3);

La columna muestra los-valores para las tres marcas de palomitas de maíz (0,0000), los dos tipos de Popper (0,0001), y la interacción entre la marca y el tipo de Popper (0,7462).Prob>Fp Estos valores indican que la marca de palomitas de maíz y el tipo de Popper afectan el rendimiento de las palomitas de maíz, pero no hay evidencia de un efecto de interacción de los dos.

Visualice la matriz de celdas que contiene la tabla ANOVA.

tbl
tbl = 6x6 cell array
  Columns 1 through 5

    {'Source'     }    {'SS'     }    {'df'}    {'MS'      }    {'F'       }
    {'Columns'    }    {[15.7500]}    {[ 2]}    {[  7.8750]}    {[ 56.7000]}
    {'Rows'       }    {[ 4.5000]}    {[ 1]}    {[  4.5000]}    {[ 32.4000]}
    {'Interaction'}    {[ 0.0833]}    {[ 2]}    {[  0.0417]}    {[  0.3000]}
    {'Error'      }    {[ 1.6667]}    {[12]}    {[  0.1389]}    {0x0 double}
    {'Total'      }    {[     22]}    {[17]}    {0x0 double}    {0x0 double}

  Column 6

    {'Prob>F'    }
    {[7.6790e-07]}
    {[1.0037e-04]}
    {[    0.7462]}
    {0x0 double  }
    {0x0 double  }

Almacene el-estadístico para los factores y la interacción del factor en variables separadas.F

Fbrands = tbl{2,5}
Fbrands = 56.7000 
Fpoppertype = tbl{3,5}
Fpoppertype = 32.4000 
Finteraction = tbl{4,5}
Finteraction = 0.3000 

Cargue los datos de ejemplo.

load popcorn popcorn
popcorn = 6×3

    5.5000    4.5000    3.5000
    5.5000    4.5000    4.0000
    6.0000    4.0000    3.0000
    6.5000    5.0000    4.0000
    7.0000    5.5000    5.0000
    7.0000    5.0000    4.5000

Los datos provienen de un estudio de marcas de palomitas de maíz y tipos de Popper (Hogg 1987). Las columnas de la matriz son marcas (gourmet, nacional y genérica).popcorn Las hileras son de tipo Popper aceite y aire. En el estudio, los investigadores reventaron un lote de cada marca tres veces con cada Popper. Los valores son el rendimiento en tazas de palomitas de maíz reventado.

Realice un ANOVA de dos vías. Calcule también las estadísticas que necesita para realizar una prueba de comparación múltiple en los efectos principales.

[~,~,stats] = anova2(popcorn,3,'off')
stats = struct with fields:
      source: 'anova2'
     sigmasq: 0.1389
    colmeans: [6.2500 4.7500 4]
        coln: 6
    rowmeans: [4.5000 5.5000]
        rown: 9
       inter: 1
        pval: 0.7462
          df: 12

La estructura incluyestats

  • El error medio cuadrado ()sigmasq

  • Las estimaciones del rendimiento medio para cada marca de palomitas de maíz ()colmeans

  • El número de observaciones para cada marca de palomitas de maíz ()coln

  • La estimación del rendimiento medio para cada tipo de Popper ()rowmeans

  • El número de observaciones para cada tipo de Popper ()rown

  • El número de interacciones ()inter

  • El-valor que muestra el nivel de significancia del término de interacción ()ppval

  • Los grados de error de libertad ().df

Realice una prueba de comparación múltiple para ver si el rendimiento de las palomitas de maíz difiere entre pares de marcas de palomitas de maíz (columnas).

c = multcompare(stats)
Note: Your model includes an interaction term.  A test of main effects can be  difficult to interpret when the model includes interactions. 

c = 3×6

    1.0000    2.0000    0.9260    1.5000    2.0740    0.0000
    1.0000    3.0000    1.6760    2.2500    2.8240    0.0000
    2.0000    3.0000    0.1760    0.7500    1.3240    0.0116

Las dos primeras columnas de mostrar los grupos que se comparan.c La cuarta columna muestra la diferencia entre el medio de grupo estimado. Las columnas tercera y quinta muestran los límites inferior y superior de los intervalos de confianza del 95% para la diferencia media real. La sexta columna contiene el-valor de una prueba de hipótesis de que la diferencia media correspondiente es igual a cero.p Todos los valores (0, 0 y 0,0116) son muy pequeños, lo que indica que el rendimiento de las palomitas de maíz difiere en las tres marcas.p

La figura muestra la comparación múltiple de los medios. De forma predeterminada, la media del grupo 1 se resalta y el intervalo de comparación está en azul. Debido a que los intervalos de comparación para los otros dos grupos no se cruzan con los intervalos para el grupo 1 significan, se resaltan en rojo. Esta falta de intersección indica que ambos medios son diferentes de la media del grupo 1. Seleccione otros medios de grupo para confirmar que todos los medios de grupo son significativamente diferentes entre sí.

Realice una prueba de comparación múltiple para ver el rendimiento de las palomitas de maíz difiere entre los dos tipos de Popper (filas).

c = multcompare(stats,'Estimate','row')
Note: Your model includes an interaction term.  A test of main effects can be  difficult to interpret when the model includes interactions. 

c = 1×6

    1.0000    2.0000   -1.3828   -1.0000   -0.6172    0.0001

El pequeño valor de 0,0001 indica que el rendimiento de las palomitas de maíz difiere entre los dos tipos de Popper (aire y aceite).p La figura muestra los mismos resultados. Los intervalos de comparación disjuntas indican que los medios del grupo son significativamente diferentes entre sí.

Argumentos de entrada

contraer todo

Datos de ejemplo, especificados como una matriz. Las columnas corresponden a grupos de un factor, y las filas corresponden a los grupos del otro factor y las replicaciones. Las replicaciones son las mediciones o observaciones para cada combinación de grupos (niveles) de la fila y el factor de columna. Por ejemplo, en los siguientes datos, el factor de fila tiene tres niveles, el factor de columna tiene dos niveles y hay dos replicaciones ().ABreps = 2 Los subsubgrupos indican fila, columna y replicación, respectivamente.

B=1B=2[y111y121y112y122y211y221y212y222y311y321y312y322]}A=1}A=2}A=3

Tipos de datos: single | double

Número de replicaciones para cada combinación de grupos, especificadas como un número entero. Por ejemplo, los siguientes datos tienen dos replicaciones () para cada combinación de grupo de factor de fila y factor de columna.reps = 2AB

B=1B=2[y111y121y112y122y211y221y212y222y311y321y312y322]}A=1}A=2}A=3

  • Cuando es (valor predeterminado), devuelve dos valores en Vector:reps1anova2pp

    • El-valor de la hipótesis nula de que todas las muestras del factor (es decir, todas las muestras de columna) se dibujan de la misma población.pBy

    • El-valor para la hipótesis nula, que todas las muestras de factor (es decir, todas las muestras de fila en) se dibujan de la misma población.pAy

  • Cuando es mayor que, también devuelve el-valor para la hipótesis nula que los factores y no tienen interacción (es decir, los efectos debido a los factores y son).reps1anova2pABABadditive

Ejemplo: Especifica que cada combinación de grupos (niveles) tiene tres replicaciones.p = anova(y,3)

Tipos de datos: single | double

Indicador para mostrar la tabla ANOVA como una figura, especificada como o.'on''off'

Argumentos de salida

contraer todo

-Value para el-test, devuelto como un valor escalar.pF Un valor pequeño indica que los resultados son estadísticamente significativos.p Los niveles de significancia comunes son 0,05 o 0,01. Por ejemplo:

  • Un valor suficientemente pequeño para la hipótesis nula para los medios de grupo de factor de fila sugiere que al menos una media de la muestra de fila es significativamente diferente de los otros medios de muestra de fila; es decir, hay un efecto principal debido al factorpAA

  • Un valor suficientemente pequeño para la hipótesis nula para los medios de grupo (nivel) de factor de columna sugiere que al menos una media de muestra de columna es significativamente diferente de los otros medios de muestra de columna; es decir, hay un efecto principal debido al factor.pBB

  • Un valor suficientemente pequeño para las combinaciones de grupos (niveles) de factores y sugiere que hay una interacción entre los factores y.pABAB

Tabla ANOVA, devuelta como una matriz de celdas. tiene seis columnas.tbl

Nombre de columnaDefinición
sourceFuente de la variabilidad.
SSSuma de los cuadrados debidos a cada fuente.
dfGrados de libertad asociados a cada fuente.
MSCuadrados de media para cada fuente, que es la relación.SS/df
F-Estadística, que es la relación de los cuadrados medio.F
Prob>F-Value, que es la probabilidad de que la-estadística puede tomar un valor mayor que el valor de la estadística de prueba calculada. deriva esta probabilidad de la CDF de la distribución.pFanova2F

Las filas de la tabla ANOVA muestran la variabilidad en los datos, dividida por la fuente en tres o cuatro partes, dependiendo del valor de.reps

FilaDefinición
ColumnsLa variabilidad debido a las diferencias entre la columna significa
RowsLa variabilidad debido a las diferencias entre la fila significa
Interaction

Variabilidad debida a la interacción entre las filas y las columnas (si es mayor que su valor predeterminado de 1)reps

ErrorLa variabilidad restante no explicada por ninguna fuente sistemática

Tipos de datos: cell

Estadística para, devuelta como una estructura.pruebas de comparaciones múltiples Se usa para realizar varias pruebas de comparación, suministrando como argumento de entrada. tiene nueve campos.multcomparestatsstats

Campo Definición
sourceFuente de la salidastats
sigmasqError cuadrado medio
colmeansLos valores estimados de la columna significan
colnNúmero de observaciones para cada grupo en columnas
rowmeansLos valores estimados de la fila significan
rownNúmero de observaciones para cada grupo en filas
interNúmero de interacciones
pval-valor para el término de interacciónp
dfGrados de libertad de error (— 1) * * donde está el número de replicaciones y y es el número de grupos en los factores, respectivamente.repsrcrepscr

Tipos de datos: struct

Referencias

[1] Hogg, R. V., and J. Ledolter. Engineering Statistics. New York: MacMillan, 1987.

Introducido antes de R2006a