Оценка ядра ggplot2 logdensity - PullRequest
       11

Оценка ядра ggplot2 logdensity

2 голосов
/ 25 ноября 2011

Как лучше всего добиться того же результата, что и: Оценка ядра log-pdf

используя ggplot2 вместо базовой графики? Я хотел бы сделать это с чем-то похожим на stat_density, так как в противном случае все быстро усложняется в сочетании с огранкой и т. Д.

Ответы [ 2 ]

3 голосов
/ 25 ноября 2011
library(ggplot2)

x = rnorm(100)

dev.new(width=4, height=4)
qplot(x, stat='density', geom='line') + scale_y_log10()

enter image description here

0 голосов
/ 25 ноября 2011

Не могли бы вы просто присвоить dlog_y= log(density(x)$y) и вызвать gglplot с этим в качестве вашей переменной y?

x = rnorm(100); dlog_y= log()
dev.new(width=4, height=4)
qplot(x = density(x)$x, y=dlog_y)
...