Этот поток Почему я не могу получить p-значение меньше, чем 2.2e-16? , начинает отвечать на мою проблему, однако моя проблема заключается в использовании F-распределения для генерации p-значения меньше.чем 2.2e-16.
1-pf(50,3,200)
[1] 0
Моя цель состоит в том, чтобы эффективно имитировать levene.test из пакета lawstat, но с использованием функции "lite", чтобы повысить эффективность вычислений на 10-100 000.Я достиг этого, но у меня есть проблема с несколькими значениями p <2.2e-16.Функция <code>levene.test, которая использует функцию anova, способна справиться с этим.
Кроме того, что я возвращаюсь к моей функции и редактирую ее для использования функции anova поверх pf, кто-нибудь можетесть идеи, почему pf ограничен anova нет.Если бы я использовал функцию anova, это удвоило бы время выполнения