从 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 值。

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

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

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

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

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