tick - r plot title




兩個密度圖之間的差異 (2)

看起來你需要花一點時間學習如何使用R(或任何其他語言,就此而言)。 幫助文件是你的朋友。 從?density的輸出:

值[即函數返回的數據]

如果give.Rkern為真,則為數R(K),否則為具有類“密度”的對象,其基礎結構是包含以下組件的列表。

x估計密度點的n個坐標。

y估計的密度值。 這些將是非負的,但可以是零[為簡潔而刪除的“價值”的剩餘部分]

那麼,做:

foo<- density(data1) 
bar<- density(data2)
plot(foo$y-bar$y) 

有沒有簡單的方法來繪製兩個概率密度函數之間的差異?

我可以在同一個圖上繪製我的數據集(都是一維向量,大概有11000個值)的pdf,一起來得到重疊/差異的想法,但是如果我能夠看到區別。

一些沿著下面的線(雖然這顯然不工作):

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

我對R比較陌生,在任何論壇上都找不到我要找的東西。

提前致謝


這應該工作:

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上回答類似問題時沒有採取這一步驟的錯誤。 太糟糕了,我不記得正確的論點。





plot