Загрузите образцы данных.

load fisheriris

Вектор столбца species состоит из цветков радужки трёх различных видов: сетозы, версиколора и виргиники. Двойная матрица meas состоит из четырёх видов измерений на цветках: длины и ширины чашелистиков и лепестков в сантиметрах соответственно.

Храните данные в табличном массиве.

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'});

Подгоните модель повторных измерений, где измерения являются откликами, а вид является предикторной переменной.

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

Выполните тест Маучли, чтобы оценить предположение о сферичности.

mauchly(rm)

ans=1×4 table
       W       ChiStat    DF      pValue  
    _______    _______    __    __________

    0.55814    84.976     5     7.6149e-17

Небольшое $$$$значение p (в pValue поле) указывает, что сферичность, следовательно, предположение о симметрии соединения, не сохраняется. Для вычисления $$$$значений p для повторяющихся измерений anova следует использовать поправки epsilon. Вы можете вычислить поправки эпсилона, используя epsilon способ и выполнение повторных мероприятий anova с откорректированными $$$$значениями p с использованием ranova способ.