r - 如何比较两个概率密度函数？

Question

我有两个概率密度函数，我想知道它们的分布是否相似。我知道 R 中的 KS 测试可以做到这一点，但是当我运行代码时，会发生错误。谢谢你的帮助。

set.seed(100)
a=density(sample(x=1:30,size = 30,replace = T))
b=density(sample(x=1:40,size = 35,replace = T))
plot(a)
lines(b)

ks.test(a,b)
Error in ks.test(a, b) : 
 'y' must be numeric or a function or a string naming a valid function

Answer 1

您需要输入两个样本（数据值的数值向量）作为 的参数ks.test。

set.seed(100)
x <- sample(x=1:30,size = 30,replace = T)
y <- sample(x=1:40,size = 35,replace = T)
a=density(x)
b=density(y)
plot(a)
lines(b)

ks.test(x,y)