密度プロットをy軸にプロットする方法は？

Question

が、私はこのように、異なるサンプルサイズを持つ複数の変数の密度プロットを描画するためにRにSMパッケージを使用しています：

var1 <- density(vars1[,1]) 
var2 <- density(vars2[,1]) 
var3 <- density(vars3[,1]) 

pdf(file="density.pdf",width=8.5,height=8) 
plot(var1,col="BLUE") 
par(new=T) 
plot(var2,axes=FALSE,col="RED") 
par(new=T) 
plot(var3,axes=FALSE,col="GREEN") 
dev.off() 

私がいる問題を、私はy軸を表示したいということです私は、異なる変数をより意味のある方法で互いに比較することができるように、比率を指定します。 3つの密度プロットの最大値はまったく同じになりました.Y軸に比率が表示されていても、そうでないと確信しています。助言がありますか？どうもありがとう！

編集：

私はちょうど私が今のコードのプロットの部分は次のようになり、既存のプロットの上にプロットべきではないということを学びました。しかし

pdf(file="density.pdf",width=8.5,height=8) 
plot(var1,col="BLUE") 
lines(var2,col="RED") 
lines(var3,col="GREEN") 
dev.off() 

これらの行の最大値サンプルサイズの違いに非常に似ています。 3つの変数すべてに対してy軸の比率を設定する方法はありますか？そのため、曲線下の面積は3つの変数すべてで等しくなりますか？どうもありがとう！

Answer 1

軸が異なる可能性があるため、既存のプロットの上にプロットしないでください。代わりに、lines()を使用して、最初のプロット後に2番目と3番目の濃度をプロットします。必要に応じて、ylimのパラメータをplot()に調整してください。

どのようにサンプルサイズがない問題べきための例：

set.seed(1) 
    D1 <- density(rnorm(1000)) 
    D2 <- density(rnorm(10000)) 
    D3 <- density(rnorm(100000)) 
    plot(D1$x,D1$y,type='l',col="red",ylim=c(0,.45)) 
    lines(D2$x,D2$y,lty=2,col="blue") 
    lines(D3$x,D3$y,lty=3,col="green") 

Answer 2

あなたは限界にハードコーディングしないことで、もう少し柔軟なティムのソリューションを作ることができます。

plot(D1$x,D1$y,type='l',col="red",ylim=c(0, max(sapply(list(D1, D2, D3), 
    function(x) {max(x$y)})))) 

これは、密度関数が必ずしもその範囲内に制約されていないことをVincentが指摘している。