此线程 Why can't I get a p-value smaller than 2.2e-16?开始回答我的问题，但是我的问题是使用F分布生成小于2.2e-16的p值。
1-pf(50,3,200)
[1] 0
我的目标是有效地模仿levene.test包中的lawstat，但在应用程序中包装“ lite”功能可以更有效地进行10-100,000的计算。
我已经实现了，但是有一些p值<2.2e-16的问题。使用levene.test函数的anova函数可以处理这些问题。
除了返回我的功能并在anova上使用pf函数对其进行编辑之外，还有谁知道pf受anova限制的原因吗？没错
如果我使用anova函数，它将使运行时间加倍