r - ggplot2 は密度を間違えますか?

Question

コードは次のとおりです。

dat = data.frame(method=gl(3, 100), res=c(rnorm(100), rnorm(100, 1, 1), rnorm(100, 2, 1)))
png('/tmp/a.png')
p = ggplot(dat)
p = p + stat_density(aes(x=res, group=method, color=as.factor(method)), geom='line')
print(p)
dev.off()

png('/tmp/b.png')
res1 = dat[dat$method==1, ]
res2 = dat[dat$method==2, ]
res3 = dat[dat$method==3, ]
plot(density(res1))
lines(density(res2$res), col='green')
lines(density(res3$res), col='red')
dev.off()

結果：

2 番目の図を使用するplot()と、正しいことがわかります。

Answer 1

デフォルトのstat_density()位置は"stack"- なので、これらの 3 行は積み上げられます。plot()使用中と同じ結果を得るにはposition="identity"。

ggplot(dat)+ stat_density(aes(x=res, group=method, color=as.factor(method)), 
         geom='line',position="identity")