Main Content

N-Way ANOVA

Introduction to N-Way ANOVA

You can use the function anovan to perform N-way ANOVA. Use N-way ANOVA to determine if the means in a set of data differ with respect to groups (levels) of multiple factors. By default, anovan treats all grouping variables as fixed effects. For an example of ANOVA with random effects, see ANOVA with Random Effects. For repeated measures, see fitrm and ranova.

N-way ANOVA is a generalization of two-way ANOVA. For three factors, for example, the model can be written as

yijkr=μ+αi+βj+γk+(αβ)ij+(αγ)ik+(βγ)jk+(αβγ)ijk+εijkr,

where

  • yijkr is an observation of the response variable. i represents group i of factor A, i = 1, 2, ..., I, j represents group j of factor B, j = 1, 2, ..., J, k represents group k of factor C, and r represents the replication number, r = 1, 2, ..., R. For constant R, there are a total of N = I*J*K*R observations, but the number of observations does not have to be the same for each combination of groups of factors.

  • μ is the overall mean.

  • αi are the deviations of groups of factor A from the overall mean μ due to factor A. The values of αi sum to 0.

    i=1Iαi=0.

  • βj are the deviations of groups in factor B from the overall mean μ due to factor B. The values of βj sum to 0.

    j=1Jβj=0.

  • γk are the deviations of groups in factor C from the overall mean μ due to factor C. The values of γk sum to 0.

    k=1Kγk=0.

  • (αβ)ij is the interaction term between factors A and B. (αβ)ij sum to 0 over either index.

    i=1I(αβ)ij=j=1J(αβ)ij=0.

  • (αγ)ik is the interaction term between factors A and C. The values of (αγ)ik sum to 0 over either index.

    i=1I(αγ)ik=k=1K(αγ)ik=0.

  • (βγ)jk is the interaction term between factors B and C. The values of (βγ)jk sum to 0 over either index.

    j=1J(βγ)jk=k=1K(βγ)jk=0.

  • (αβγ)ijk is the three-way interaction term between factors A, B, and C. The values of (αβγ)ijk sum to 0 over any index.

    i=1I(αβγ)ijk=j=1J(αβγ)ijk=k=1K(αβγ)ijk=0.

  • εijkr are the random disturbances. They are assumed to be independent, normally distributed, and have constant variance.

Three-way ANOVA tests hypotheses about the effects of factors A, B, C, and their interactions on the response variable y. The hypotheses about the equality of the mean responses for groups of factor A are

H0:α1=α2=αIH1: at least one αi is different, i=1, 2, ..., I.

The hypotheses about the equality of the mean response for groups of factor B are

H0:β1=β2==βJH1: at least one βj is different,  j=1, 2, ..., J.

The hypotheses about the equality of the mean response for groups of factor C are

H0:γ1=γ2==γKH1: at least one γk is different, k=1, 2, ..., K.

The hypotheses about the interaction of the factors are

H0:(αβ)ij=0H1:at least one (αβ)ij0

H0:(αγ)ik=0H1:at least one (αγ)ik0H0:(βγ)jk=0H1:at least one (βγ)jk0H0:(αβγ)ijk=0H1:at least one (αβγ)ijk0

In this notation parameters with two subscripts, such as (αβ)ij, represent the interaction effect of two factors. The parameter (αβγ)ijk represents the three-way interaction. An ANOVA model can have the full set of parameters or any subset, but conventionally it does not include complex interaction terms unless it also includes all simpler terms for those factors. For example, one would generally not include the three-way interaction without also including all two-way interactions.

Prepare Data for N-Way ANOVA

Unlike anova1 and anova2, anovan does not expect data in a tabular form. Instead, it expects a vector of response measurements and a separate vector (or text array) containing the values corresponding to each factor. This input data format is more convenient than matrices when there are more than two factors or when the number of measurements per factor combination is not constant.

y=[y1,y2,y3,y4,y5,,yN]g1={'A','A','C','B','B',,'D'}g2=[12131,2]g3={'hi','mid','low','mid','hi',,'low'}

Perform N-Way ANOVA

This example shows how to perform N-way ANOVA on car data with mileage and other information on 406 cars made between 1970 and 1982.

Load the sample data.

load carbig

The example focusses on four variables. MPG is the number of miles per gallon for each of 406 cars (though some have missing values coded as NaN). The other three variables are factors: cyl4 (four-cylinder car or not), org (car originated in Europe, Japan, or the USA), and when (car was built early in the period, in the middle of the period, or late in the period).

Fit the full model, requesting up to three-way interactions and Type 3 sums-of-squares.

varnames = {'Origin';'4Cyl';'MfgDate'};
anovan(MPG,{org cyl4 when},3,3,varnames);

Note that many terms are marked by a # symbol as not having full rank, and one of them has zero degrees of freedom and is missing a p-value. This can happen when there are missing factor combinations and the model has higher-order terms. In this case, the cross-tabulation below shows that there are no cars made in Europe during the early part of the period with other than four cylinders, as indicated by the 0 in tbl(2,1,1).

[tbl,chi2,p,factorvals] = crosstab(org,when,cyl4)
tbl = 
tbl(:,:,1) =

    82    75    25
     0     4     3
     3     3     4


tbl(:,:,2) =

    12    22    38
    23    26    17
    12    25    32

chi2 = 207.7689
p = 8.0973e-38
factorvals=3×3 cell array
    {'USA'   }    {'Early'}    {'Other'   }
    {'Europe'}    {'Mid'  }    {'Four'    }
    {'Japan' }    {'Late' }    {0x0 double}

Consequently it is impossible to estimate the three-way interaction effects, and including the three-way interaction term in the model makes the fit singular.

Using even the limited information available in the ANOVA table, you can see that the three-way interaction has a p-value of 0.699, so it is not significant.

Examine only two-way interactions.

[p,tbl2,stats,terms] = anovan(MPG,{org cyl4 when},2,3,varnames);

terms
terms = 6×3

     1     0     0
     0     1     0
     0     0     1
     1     1     0
     1     0     1
     0     1     1

Now all terms are estimable. The p-values for interaction term 4 (Origin*4Cyl) and interaction term 6 (4Cyl*MfgDate) are much larger than a typical cutoff value of 0.05, indicating these terms are not significant. You could choose to omit these terms and pool their effects into the error term. The output terms variable returns a matrix of codes, each of which is a bit pattern representing a term.

Omit terms from the model by deleting their entries from terms.

terms([4 6],:) = []
terms = 4×3

     1     0     0
     0     1     0
     0     0     1
     1     0     1

Run anovan again, this time supplying the resulting vector as the model argument. Also return the statistics required for multiple comparisons of factors.

[~,~,stats] = anovan(MPG,{org cyl4 when},terms,3,varnames)

stats = struct with fields:
         source: 'anovan'
          resid: [3.1235 0.1235 3.1235 1.1235 2.1235 0.1235 -0.8765 -0.8765 -0.8765 0.1235 NaN NaN NaN NaN NaN 0.1235 -0.8765 NaN 0.1235 -0.8765 -2.3832 7.1235 3.1235 6.1235 0.6168 1.2857 0.2857 -0.7143 0.2857 1.2857 6.1235 -4.8765 ... ] (1x406 double)
         coeffs: [18x1 double]
            Rtr: [10x10 double]
       rowbasis: [10x18 double]
            dfe: 388
            mse: 14.1056
    nullproject: [18x10 double]
          terms: [4x3 double]
        nlevels: [3x1 double]
     continuous: [0 0 0]
         vmeans: [3x1 double]
       termcols: [5x1 double]
     coeffnames: {18x1 cell}
           vars: [18x3 double]
       varnames: {3x1 cell}
       grpnames: {3x1 cell}
        vnested: []
            ems: [5x5 double]
          denom: []
        dfdenom: []
        msdenom: []
         varest: []
          varci: []
       txtdenom: []
         txtems: []
        rtnames: []

Now you have a more parsimonious model indicating that the mileage of these cars seems to be related to all three factors, and that the effect of the manufacturing date depends on where the car was made.

Perform multiple comparisons for Origin and Cylinder.

[results,~,~,gnames] = multcompare(stats,'Dimension',[1,2]);

Display the multiple comparison results and the corresponding group names in a table.

tbl = array2table(results,"VariableNames", ...
    ["Group A","Group B","Lower Limit","A-B","Upper Limit","P-value"]);
tbl.("Group A") = gnames(tbl.("Group A"));
tbl.("Group B") = gnames(tbl.("Group B"))
tbl=15×6 table
              Group A                         Group B               Lower Limit      A-B      Upper Limit     P-value  
    ____________________________    ____________________________    ___________    _______    ___________    __________

    {'Origin=USA,4Cyl=Other'   }    {'Origin=Japan,4Cyl=Other' }      -5.4891      -3.8412      -2.1932      4.2334e-10
    {'Origin=USA,4Cyl=Other'   }    {'Origin=Europe,4Cyl=Other'}      -4.4146      -2.7251      -1.0356      6.2974e-05
    {'Origin=USA,4Cyl=Other'   }    {'Origin=USA,4Cyl=Four'    }      -9.9992      -8.5828      -7.1664               0
    {'Origin=USA,4Cyl=Other'   }    {'Origin=Japan,4Cyl=Four'  }      -14.024      -12.424      -10.824               0
    {'Origin=USA,4Cyl=Other'   }    {'Origin=Europe,4Cyl=Four' }      -12.898      -11.308       -9.718               0
    {'Origin=Japan,4Cyl=Other' }    {'Origin=Europe,4Cyl=Other'}     -0.71714        1.116       2.9492          0.5085
    {'Origin=Japan,4Cyl=Other' }    {'Origin=USA,4Cyl=Four'    }      -7.3655      -4.7417      -2.1179      3.8678e-06
    {'Origin=Japan,4Cyl=Other' }    {'Origin=Japan,4Cyl=Four'  }      -9.9992      -8.5828      -7.1664               0
    {'Origin=Japan,4Cyl=Other' }    {'Origin=Europe,4Cyl=Four' }      -9.7464      -7.4668      -5.1872      1.4557e-20
    {'Origin=Europe,4Cyl=Other'}    {'Origin=USA,4Cyl=Four'    }      -8.5396      -5.8577      -3.1757      6.9888e-09
    {'Origin=Europe,4Cyl=Other'}    {'Origin=Japan,4Cyl=Four'  }      -12.052      -9.6988      -7.3459               0
    {'Origin=Europe,4Cyl=Other'}    {'Origin=Europe,4Cyl=Four' }      -9.9992      -8.5828      -7.1664               0
    {'Origin=USA,4Cyl=Four'    }    {'Origin=Japan,4Cyl=Four'  }      -5.4891      -3.8412      -2.1932      4.2334e-10
    {'Origin=USA,4Cyl=Four'    }    {'Origin=Europe,4Cyl=Four' }      -4.4146      -2.7251      -1.0356      6.2974e-05
    {'Origin=Japan,4Cyl=Four'  }    {'Origin=Europe,4Cyl=Four' }     -0.71714        1.116       2.9492          0.5085

See Also

| | | |

Related Examples

More About