两个密度图之间的差异

Question

是否有一种简单的方法来绘制两个概率密度函数之间的差异？

我可以在同一个绘图上绘制我的数据集的pdf(两者都是一维向量,大约有11000个值),以便了解重叠/差异,但如果我能看到一个情节,对我来说会更有用.区别.

下面的内容(虽然这显然不起作用):

> plot(density(data1)-density(data2))

我对R比较陌生,并且无法在任何论坛上找到我想要的东西.

提前致谢

Answer 1

这应该工作:

plot(x =density(data1, from= range(c(data1, data2))[1], 
                       to=range(c(data1, data2))[2] )$x, 
  y=  density(data1, from= range(c(data1, data2))[1], 
                     to=range(c(data1, data2))[2] )$y-
       density(data2,  from= range(c(data1, data2))[1], 
                      to=range(c(data1, data2))[2] )$y )

诀窍是确保密度具有相同的限制.然后你可以在相同的位置绘制他们的差异.我对相同限制的需要的理解来自几年前在Rhelp回答类似问题时没有采取这一步骤的错误.太糟糕了,我记不起正确的论点了.