I want to create the empirical cumulative distribution function for two samples and put the plots in the same plot [R]

Question

I am using this code to generate the empirical cumulative distribution function for the two samples (you can put any numerical values in them). I would like to put them in the same plot but if you run the following commands everything is overlapping really bad [see picture 1]. Is there any way to do it like this [see picture 2] (also I want the symbols to disappear and be a line like the picture 2) .

plot(ecdf(sample[,1]),pch = 1)
par(new=TRUE)        
plot(ecdf(sample[,2]),pch = 2)

picture 1:https://www.dropbox.com/s/sg1fr8jydsch4xp/vanboeren2.png?dl=0

picture 2:https://www.dropbox.com/s/erhgla34y5bxa58/vanboeren1.png?dl=0

Update: I am doing this

  df1 <- data.frame(x = sample[,1]) 
  df2 <- data.frame(x = sample[,2])   
  ggplot(df1, aes(x, colour = "g")) + stat_ecdf()
  +geom_step(data = df2) 
  scale_x_continuous(limits = c(0, 5000)) `

which is very close (in terms of shape) but still can not put them at the same plot.

you have to set the same limits to the x axis for both plots. See if [this](http://stackoverflow.com/questions/42517471/my-parabola-shifts-along-the-x-axis-as-i-change-x-lim/42517775#42517775) helps. — R. Schifini, Apr 19 '17 at 16:33
[Here](http://stackoverflow.com/questions/3606697/how-to-set-limits-for-axes-in-ggplot2-r-plots) is one answer that may help — R. Schifini, Apr 19 '17 at 16:35
still can not do it :) if you can be more specific please. Just copy paste what I have written in the update section and show me please :) — giorgos kapetanidis, Apr 19 '17 at 16:38
If your question is about how to plot multiple datasets in the same plot with `ggplot`, see [here](http://stackoverflow.com/questions/9109156/ggplot-combining-two-plots-from-different-data-frames) — aosmith, Apr 19 '17 at 16:46

score 0 · Accepted Answer · answered Apr 19 '17 at 17:20

0

Try this with basic plot:

df1 <- data.frame(x = runif(200,1,5)) 
df2 <- data.frame(x = runif(200,3,8))

plot(ecdf(df1[,1]),pch = 1, xlim=c(0,10), main=NULL)
par(new=TRUE)        
plot(ecdf(df2[,1]),pch = 2, xlim=c(0,10), main=NULL)

Both graphs have now the same xlim (try removing it to see both superimposed incorrectly). The main=NULL removes the title

Result:

answered Apr 19 '17 at 17:20

R. Schifini

8,782
2
24
31

do you know how to put legends ? I know how to do it with colors, but how can I do it here? – giorgos kapetanidis Apr 19 '17 at 17:30
`legends` is the function that sets the legends on a plot. – R. Schifini Apr 20 '17 at 04:08

I want to create the empirical cumulative distribution function for two samples and put the plots in the same plot [R]

1 Answers1