4

我有一个分布,例如:

d 
#[1]  4 22 15  5  9  5 11 15 21 14 14 23  6  9 17  2  7 10  4

d或者,dput格式中的向量。

d <- c(4, 22, 15, 5, 9, 5, 11, 15, 21, 14, 14, 23, 6, 9, 17, 2, 7, 10, 4)

当我应用 ks.test 时,:

gamma <- ks.test(d, "pgamma", shape = 3.178882, scale = 3.526563) 

这给出了以下警告:

警告消息:在 ks.test(d, "pgamma", shape = 3.178882, scale = 3.526563) 中:Kolmogorov-Smirnov 测试不应存在领带

我试过 put unique(d),但很明显我的数据减少了值,我不希望这种情况发生。
和网上的其他方式和例子一样,这个例子也有,但不同的是测试显示一些带有警告消息的结果,而不仅仅是没有值的消息ks.test

一些帮助?

4

1 回答 1

5

gamma您可以找到您的结果中,警告消息没有被阻止

d <- c(4, 22, 15, 5, 9, 5, 11, 15, 21, 14, 14, 23, 6, 9, 17, 2, 7, 10, 4)
gamma <- ks.test(d, "pgamma", shape = 3.178882, scale = 3.526563)

警告消息:在 ks.test(d, "pgamma", shape = 3.178882, scale = 3.526563) 中:Kolmogorov-Smirnov 测试不应存在领带

gamma

    One-sample Kolmogorov-Smirnov test

data:  d
D = 0.14549, p-value = 0.816
alternative hypothesis: two-sided

您可以在帮助页面中找到对警告的解释??ks.test

关系的存在总是会产生警告,因为连续分布不会产生它们。如果通过四舍五入产生的联系可能近似有效,但即使是适度的四舍五入也会对计算的统计数据产生重大影响。

如您所见,应用了一些舍入,并且测试“大约”有效。

于 2018-08-23T15:38:15.977 回答