从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值，但是我找不到任何用于R的A-D测试。

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

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