mauchly

Класс: RepeatedMeasuresModel

Тест Мочли для шарообразности

расширьте все на странице

Синтаксис

tbl = mauchly(rm)

tbl = mauchly(rm,C)

Описание

пример

tbl = mauchly(rm) возвращает результат теста Мочли для шарообразности для модели rm повторных измерений.

Это тестирует нулевую гипотезу, что предположение шарообразности верно для переменных отклика в rm.

Для получения дополнительной информации смотрите Тест Мочли Шарообразности.

tbl = mauchly(rm,C) возвращает результат теста Мочли на основе контрастного матричного C.

Входные параметры

развернуть все

`rm` — Модель повторных измерений
`RepeatedMeasuresModel` объект

Модель повторных измерений, возвращенная как RepeatedMeasuresModel объект.

Для свойств и методов этого объекта, смотрите RepeatedMeasuresModel.

`C` — Контрасты
матрица

Контрасты в виде матрицы. Значение по умолчанию C Q, включают разложение QR матрицы M, где M задан так, чтобы Y*M был различием между всеми последовательными парами столбцов матрицы Y повторных измерений.

Типы данных: single | double

Выходные аргументы

развернуть все

`tbl` — Результаты теста Мочли шарообразности
таблица

Результаты теста Мочли для шарообразности для модели rm повторных измерений, возвращенный как table.

tbl содержит следующие столбцы.

ColumnName	Определение
`W`	Значение статистической величины W Мочли
`ChiStat`	Значение статистической величины хи-квадрата
`DF`	Степени свободы статистической величины Хи-квадрата
`pValue`	p-, соответствующее статистической величине Хи-квадрата

Типы данных: table

Примеры

развернуть все

Выполните тест Мочли

Попробовать в MATLAB

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

load fisheriris

Вектор-столбец species состоит из ирисовых цветов трех различных разновидностей: setosa, versicolor, и virginica. Двойной матричный 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 метод и выполняет повторные измерения anova с откорректированным $p$ - значения с помощью ranova метод.