Как лучше всего добиться того же результата, что и: Оценка ядра log-pdf
используя ggplot2 вместо базовой графики? Я хотел бы сделать это с чем-то похожим на stat_density, так как в противном случае все быстро усложняется в сочетании с огранкой и т. Д.
stat_density
library(ggplot2) x = rnorm(100) dev.new(width=4, height=4) qplot(x, stat='density', geom='line') + scale_y_log10()
Не могли бы вы просто присвоить dlog_y= log(density(x)$y) и вызвать gglplot с этим в качестве вашей переменной y?
dlog_y= log(density(x)$y)
x = rnorm(100); dlog_y= log() dev.new(width=4, height=4) qplot(x = density(x)$x, y=dlog_y)