r - 具有 Kolmogorov-Smirnov 双样本、Cramer-von Mises 和 Kuiper 检验的 ECDF

Question

在创建等级检验的经验分布时，我有兴趣重新创建 SAS NPAR1WAY 结果。

我找到了以下代码，了解如何使用Easier way 中的 iris 数据集创建等级图的经验分布，以在 ggplot 中绘制累积频率分布？

library(plyr)
data(iris)
#Ecdf within species
iris.species <- ddply(iris, .(Species), summarize,
Sepal.Length = unique(Sepal.Length),
ecdf = ecdf(Sepal.Length)(unique(Sepal.Length)))
ggplot(iris.species, aes(Sepal.Length, ecdf, color = Species)) + geom_step()

但我不确定如何对按物种分类的 Sepal.Length 的 Kolmogorov-Smirnov 检验、Cramer-von Mises 检验或 Kuiper 检验进行编码。

任何帮助将不胜感激！

标签： r