r - 从 R 中的 fitdistrplus 获取 p 值

Question

我不确定如何从 gofstat 函数中获取 p 值。我确实看到了可能有“未计算”、“未拒绝”和“拒绝”的结论的**测试。任何帮助都会很棒。假设使用的数据是连续的。

使用下面的代码

library (fitdistrplus)
fitw <- fitdist(data, "weibull",method='mle')
fw<-gofstat(fitw)

#How to get p-value from this
fw$ad
fw$adtest

Answer 1

试试fw$chisqpvalue。这将为您提供与卡方统计相关的 p 值。

注意：从文档中可以看出，AD 测试和 KS 测试仅在您有连续数据的情况下运行。

编辑：没有注意到fw$ad，等只给出了测试统计数据。

编辑 2：嗯，在搜索了很长一段时间后，我唯一能找到的另一件事是使用dgof的可能替代方案。如果您可以以正确的格式获取数据，那么根据您的分布，这有望为您提供 KS 和 CVM 的正确 p-val，但我找不到任何用于 R 的 AD 测试的东西。

您也可以尝试向作者发送电子邮件fitdistrplus并建议该函数也返回 p 值的更改？

除此之外，我认为没有任何方法可以直接提取/计算 p 值，除非阅读测试所基于的论文并创建自己的函数。抱歉，我无法提供更多帮助。