Load the sample data.

load fisheriris

The column vector `species`consists of iris
flowers of three different species: setosa, versicolor, and virginica.
The double matrix `meas` consists of four types of
measurements on the flowers: the length and width of sepals and petals
in centimeters, respectively.

Store the data in a table array.

t = table(species,meas(:,1),meas(:,2),meas(:,3),meas(:,4),...
'VariableNames',{'species','meas1','meas2','meas3','meas4'});
Meas = dataset([1 2 3 4]','VarNames',{'Measurements'});

Fit a repeated measures model, where the measurements
are the responses and the species is the predictor variable.

rm = fitrm(t,'meas1-meas4~species','WithinDesign',Meas);

Perform Mauchly's test to assess the sphericity
assumption.

mauchly(rm)

ans =
W ChiStat DF pValue
_______ _______ __ __________
0.55814 84.976 5 1.1102e-16

The small *p*-value (in the `pValue` field)
indicates that the sphericity, hence the compound symmetry assumption,
does not hold. You should use epsilon corrections to compute the *p*-values
for a repeated measures anova. You can compute the epsilon corrections
using the `epsilon` method and perform the repeated
measures anova with the corrected *p*-values using
the `ranova` method.