从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值，除非阅读测试所基于的论文并创建自己的函数。抱歉，我帮不上忙。