核估計

Question

很容易積內核估計給出通過觀察向量的隨機變量的概率密度函數（pdf）：

plot(density(x)) 

我願做相同的，但對於log-pdf（概率密度函數的對數）。

Answer 1

你可以做

d <- density(x) 
plot(d$x, log(d$y), type="l") 

也許你寧願

plot(d, log="y")