从汇总调用中从R中的安德森亲爱测试（ad.test）中检索p值

Question

所以我有以下

ad_test_by_group = aggregate(data, by = group, FUN = "ad.test")

并按组输出测试统计信息。 但是我想改为与每个组关联的p值。

有没有办法做到这一点？

Answer 1

按照@ed_sans的建议，您可以尝试：

library(nortest)
data(iris)
aggregate(. ~ Species, data = iris[,c(1,5)], function(x) ad.test(x)$p.value)
     Species Sepal.Length
1     setosa    0.3352439
2 versicolor    0.4332962
3  virginica    0.1475325

第二列是pvalues。

x<-iris[iris$Species=="setosa", 1]
ad.test(x)$p.value
[1] 0.3352439