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

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)