r - 異なる歪度を持つ 2 つの重なり合う分布の理想化されたプロット?

Question

右歪度の増加の影響を示すために、一方の分布が他方より右に歪んでいる 2 つの重なり合うプロットを作成したいと考えています。

だから私は次のようなことをします

library(SuppDists)

tmp1fun <- function(m1=0,m2=1,m3=-0.2,m4=3.7,n=20000) {
    parms <- JohnsonFit(c(m1,m2,m3,m4), 'use')
    rJohnson(n,parms) }

data1.dens <- rbind(data.frame(var1=rep('now',50000),y=tmp1fun(m3=-0.2,n=50000))
    ,data.frame(var1=rep('hist',50000),y=tmp1fun(m3=0.2,n=50000)))

library(ggplot2)
ggplot(data1.dens, aes(y, fill = var1))+
geom_density(alpha = 0.5) + coord_cartesian(xlim=c(-4, 4))

しかし、それは正しく見えません。分布の平均が他の分布に対してシフトしているように見えます。また、密度プロットがもう少し滑らかに見えるようにしたいと思います。

したがって、私の質問には 2 つの部分があります。2. このプロットを作成するのに最適なグラフィック ツールは何ですか? 可能であれば、単なる線ではなく塗りつぶされたグラフが欲しいです。おそらく、よりスムーズな帯域幅が必要です。

ヒントや提案をありがとう。