Загрузите выборочные данные.

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 field) указывает, что сферичность, следовательно, предположение составной симметрии, не выполняется. Вы должны использовать коррекции epsilon, чтобы вычислить $$p$$-значения для повторных измерений anova. Вычислить коррекции эпсилона можно с помощью epsilon метод и выполните повторные измерения anova с исправленным $$p$$-значения с использованием ranova способ.