2011-06-14 4 views
8

私は2つの曲線の違いを説明するために、お互いの上にプロットしたい累積経験密度関数をいくつか持っています。 in a previous questionを指摘したとおり、関数ecdfを描画する機能は、単にplot(Ecdf())と私は細かいマニュアルページを読んで、私は次のようなものを使用して、互いの上に複数のECDFsをプロットすることができていると判断:R:異なる色の別のものの上に1つのECDFをプロットする

require(Hmisc) 
set.seed(3) 
g <- c(rep(1, 20), rep(2, 20)) 
Ecdf(c(rnorm(20), rnorm(20)), group=g) 

しかし、私の曲線は、時には少し重なり、ちょうどこのグラフを生成する上記の例のように、これはどの伝えるのは難しいことができます:

enter image description here

私は本当に異なるこれら二つのCDFの色を作りたいと思います。しかし、私はそれを行う方法を理解することはできません。任意のヒント?

+1

Guinessからコーヒー、JDに切り替える時間。 'col ='引数があります。私の編集された答えを見てください。 –

答えて

12

メモリがあれば、私は過去にこれを行っています。私が思い出したように、あなたはEcdf()がparamterisedされているので、それをトリックする必要がありました。私はhelp(ecdf)で、それはステップ関数のプロットであることを示唆しているので、2つ以上のecdfsを見積もり、プロット1を作成してlines()で注釈を付けることができます。

編集は、ヘルプページを明確にcol=引数を述べたとして、それは

R> Ecdf(c(rnorm(20), rnorm(20)), group=g, col=c('blue', 'orange')) 

のと同じくらい簡単ですが判明します。しかし私はplot.stepfun()を明示的に使ったスクリプトレットもいくつか見つけました。

+2

マニュアルにはっきりと記載されていますが、私はそれを読んだと思います。明らかに、私のギネスのガラスは画面のその部分をブロックしていました。 –

+0

うわー、彼らは今では64オンスの「麻薬」眼鏡をしていますか? –

+0

私がDUBに着陸したとき、彼らは私をギネスのIV袋に引っ掛けました。 IV袋が私の画面を塞いだのです。 –

7

各曲線を一度に1つずつ(それぞれ独自のスタイルで)追加することができます。

Ecdf(rnorm(20), lwd = 2) 
Ecdf(rnorm(20),add = TRUE, col = 'red', lty = 1) 
+1

'col'、' lty'、...がすべてベクター化されているので、1つの操作もできます。 –