R Box'M тест на гомоедастичность

Question

Я пытаюсь воспроизвести результаты линейного дискриминантного анализа из SPSS в R, и у меня возникают трудности с поиском способа выполнения теста m-box.

Единственное, что я обнаружил, - это какой-то код, размещенный на форуме, для ручной реализации процесса, но мне было интересно, если для этого нет ничего встроенного в самом языке.

Answer 1

В упаковке биоинструмент вы можете найти функцию boxM (данные, группировка) .Он выполняет M-тест Бокса на однородность ковариационных матриц, полученных из многомерных нормальных данных в соответствии с одним классификационным фактором.Тест основан на приближении хи-квадрат.

Answer 2

Существует код , который можно найти с помощью простого поиска в rseek.Обычно это не делается, потому что очень высокая чувствительность приводит к значительным p-значениям, которые могут не значить много.реализован в пакете биоинструментов с функцией boxM.Это был все еще относительно легкий поиск.И все же верно, что вам, вероятно, не стоит беспокоиться о том, чтобы использовать его как все подобные тесты.Вероятно, вам следует просто внимательно изучить ковариационную матрицу и свои предположения.