Re: [R] overlap cdf plots and add colors and etc
Hi Peter and Jorge, Finally, it works. Thank you for your guidance. ## CDF plots par(mar=c(4,4,2,1.2),oma=c(0,0,0,0),xaxs="i", yaxs="i") plot(ecdf(dt1),verticals=TRUE, horizontals=TRUE,pch=46,do.p=FALSE, col.hor=1, col.vert=1,lwd=2,) lines(ecdf(dt2),verticals=TRUE, horizontals=TRUE,pch=20,do.p=FALSE, col.hor=2, col.vert=2,lwd=2) legend(40, 0.8, legend = c("observed","fitted"), col = c(1,2), pch=c(NA,NA), lty = c(1, 1), lwd=c(3,3),bty="n") What is "plain text"? Thank you. From: Peter Ehlers Cc: Jorge Ivan Velez ; "r-help@r-project.org" Sent: Fri, November 26, 2010 12:52:56 AM Subject: Re: [R] overlap cdf plots and add colors and etc On 2010-11-24 22:24, Roslina Zakaria wrote: > Hi Jorge, > > I tried but still it does not work. Thank you for your time. > Jorge's code works perfectly well. If you prefer lines() over plot(, add = TRUE), then use lines(): plot( ecdf(rnorm(15, sd=3)), verticals=TRUE, col="black") lines(ecdf(rnorm(11, sd=3)), verticals=TRUE, col="red") legend("topright", legend = c("observed","fitted"), col = c("black","red"), pch=c(NA,NA), lty = c(1, 1), lwd=c(3,3),bty="n", pt.cex=2) (Note that "topright" is probably the worst choice you can make. And why set pt.cex when you don't have points??) The key thing is that you can't specify your two colours in the first plot() call. At that point, R has no idea that you want to add to the plot. (I understand that Greg Snow is working on a mind-reading package, but so far it can probably read only his mind, not yours.) Here is a (very brief) reminder of how to post: 1. Use an informative subject line. [[elided Yahoo spam]] 2. *Never* use the phrase "it does not work". That is meaningless. Be specific about your problem. 3. Provide *reproducible* code. We don't have your 'datobs' or 'gam_sum_gen'. 4. Try to make the code *minimal*. It's not likely that anyone cares what labels you use for your axes/title (unless that's the problem). And nobody wants to see reams of data; use rnorm(), etc, or built-in data if possible. 5. Figure out how to set your mail program to send plain text. I know the above (and more) is in the posting guide, but it seems that nobody wants to read that quite brief document. Peter Ehlers > > > From: Jorge Ivan Velez > > Cc: r-help@r-project.org > Sent: Thu, November 25, 2010 4:46:37 PM > Subject: Re: [R] overlap cdf plots and add colors and etc > > Hi Roslina, > > > Try > > par(mar=c(4,4,2,1.2),oma=c(0,0,0,0),xaxs="i", yaxs="i") > plot(ecdf(rnorm(100))) > plot(ecdf(rnorm(100)), add = TRUE, col = 2) > > HTH, > Jorge > > > On Thu, Nov 25, 2010 at 12:18 AM, Roslina Zakaria<> wrote: > > Hi r-users, >> >> I would like to overlap 2 ecdf plots. >> >> I tried this below and it gives me two plots of ecdf but just both just in >> black. >> >> par(mar=c(4,4,2,1.2),oma=c(0,0,0,0),xaxs="i", yaxs="i") >> plot(ecdf(datobs)) >> lines(ecdf(gam_sum_gen)) >> >> Then I try to add colors etc and also the legend but fail. >> >> par(mar=c(4,4,2,1.2),oma=c(0,0,0,0),xaxs="i", yaxs="i") >> plot(ecdf(datobs),main ="CDF of the sum for winter >> season-Hume",cex.axis=1.2,xlab="Rainfall (mm)", pch = >> pch,verticals=TRUE,col=c("black","red"), lty=c(1,1),ylab="Cumulative >> probability", xlim=c(0,800),lwd=1) >> lines(ecdf(gam_sum_gen)) >> legend("topright", legend = c("observed","fitted"), >> col = c("black","red"), pch=c(NA,NA), lty = c(1, 1), >> lwd=c(3,3),bty="n", pt.cex=2) >> box() >> >> I'm not sure why it doesn't show up at all. >> >> Thank you for any help given. >> >> >> >> [[alternative HTML version deleted]] >> >> [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Re: [R] overlap cdf plots and add colors and etc
On 2010-11-24 22:24, Roslina Zakaria wrote: Hi Jorge, I tried but still it does not work. Thank you for your time. Jorge's code works perfectly well. If you prefer lines() over plot(, add = TRUE), then use lines(): plot( ecdf(rnorm(15, sd=3)), verticals=TRUE, col="black") lines(ecdf(rnorm(11, sd=3)), verticals=TRUE, col="red") legend("topright", legend = c("observed","fitted"), col = c("black","red"), pch=c(NA,NA), lty = c(1, 1), lwd=c(3,3),bty="n", pt.cex=2) (Note that "topright" is probably the worst choice you can make. And why set pt.cex when you don't have points??) The key thing is that you can't specify your two colours in the first plot() call. At that point, R has no idea that you want to add to the plot. (I understand that Greg Snow is working on a mind-reading package, but so far it can probably read only his mind, not yours.) Here is a (very brief) reminder of how to post: 1. Use an informative subject line. You've done that. Good! 2. *Never* use the phrase "it does not work". That is meaningless. Be specific about your problem. 3. Provide *reproducible* code. We don't have your 'datobs' or 'gam_sum_gen'. 4. Try to make the code *minimal*. It's not likely that anyone cares what labels you use for your axes/title (unless that's the problem). And nobody wants to see reams of data; use rnorm(), etc, or built-in data if possible. 5. Figure out how to set your mail program to send plain text. I know the above (and more) is in the posting guide, but it seems that nobody wants to read that quite brief document. Peter Ehlers ____________________ From: Jorge Ivan Velez Cc: r-help@r-project.org Sent: Thu, November 25, 2010 4:46:37 PM Subject: Re: [R] overlap cdf plots and add colors and etc Hi Roslina, Try par(mar=c(4,4,2,1.2),oma=c(0,0,0,0),xaxs="i", yaxs="i") plot(ecdf(rnorm(100))) plot(ecdf(rnorm(100)), add = TRUE, col = 2) HTH, Jorge On Thu, Nov 25, 2010 at 12:18 AM, Roslina Zakaria<> wrote: Hi r-users, I would like to overlap 2 ecdf plots. I tried this below and it gives me two plots of ecdf but just both just in black. par(mar=c(4,4,2,1.2),oma=c(0,0,0,0),xaxs="i", yaxs="i") plot(ecdf(datobs)) lines(ecdf(gam_sum_gen)) Then I try to add colors etc and also the legend but fail. par(mar=c(4,4,2,1.2),oma=c(0,0,0,0),xaxs="i", yaxs="i") plot(ecdf(datobs),main ="CDF of the sum for winter season-Hume",cex.axis=1.2,xlab="Rainfall (mm)", pch = pch,verticals=TRUE,col=c("black","red"), lty=c(1,1),ylab="Cumulative probability", xlim=c(0,800),lwd=1) lines(ecdf(gam_sum_gen)) legend("topright", legend = c("observed","fitted"), col = c("black","red"), pch=c(NA,NA), lty = c(1, 1), lwd=c(3,3),bty="n", pt.cex=2) box() I'm not sure why it doesn't show up at all. Thank you for any help given. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Re: [R] overlap cdf plots and add colors and etc
Hi Jorge, I tried but still it does not work. Thank you for your time. From: Jorge Ivan Velez Cc: r-help@r-project.org Sent: Thu, November 25, 2010 4:46:37 PM Subject: Re: [R] overlap cdf plots and add colors and etc Hi Roslina, Try par(mar=c(4,4,2,1.2),oma=c(0,0,0,0),xaxs="i", yaxs="i") plot(ecdf(rnorm(100))) plot(ecdf(rnorm(100)), add = TRUE, col = 2) HTH, Jorge On Thu, Nov 25, 2010 at 12:18 AM, Roslina Zakaria <> wrote: Hi r-users, > >I would like to overlap 2 ecdf plots. > >I tried this below and it gives me two plots of ecdf but just both just in >black. > >par(mar=c(4,4,2,1.2),oma=c(0,0,0,0),xaxs="i", yaxs="i") >plot(ecdf(datobs)) >lines(ecdf(gam_sum_gen)) > >Then I try to add colors etc and also the legend but fail. > >par(mar=c(4,4,2,1.2),oma=c(0,0,0,0),xaxs="i", yaxs="i") >plot(ecdf(datobs),main ="CDF of the sum for winter >season-Hume",cex.axis=1.2,xlab="Rainfall (mm)", pch = >pch,verticals=TRUE,col=c("black","red"), lty=c(1,1),ylab="Cumulative >probability", xlim=c(0,800),lwd=1) >lines(ecdf(gam_sum_gen)) >legend("topright", legend = c("observed","fitted"), > col = c("black","red"), pch=c(NA,NA), lty = c(1, 1), > lwd=c(3,3),bty="n", pt.cex=2) >box() > >I'm not sure why it doesn't show up at all. > >Thank you for any help given. > > > > [[alternative HTML version deleted]] > > >__ >R-help@r-project.org mailing list >https://stat.ethz.ch/mailman/listinfo/r-help >PLEASE do read the posting guide http://www.R-project.org/posting-guide.html >and provide commented, minimal, self-contained, reproducible code. > > [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Re: [R] overlap cdf plots and add colors and etc
Hi Roslina, Try par(mar=c(4,4,2,1.2),oma=c(0,0,0,0),xaxs="i", yaxs="i") plot(ecdf(rnorm(100))) plot(ecdf(rnorm(100)), add = TRUE, col = 2) HTH, Jorge On Thu, Nov 25, 2010 at 12:18 AM, Roslina Zakaria <> wrote: > Hi r-users, > > I would like to overlap 2 ecdf plots. > > I tried this below and it gives me two plots of ecdf but just both just in > black. > > par(mar=c(4,4,2,1.2),oma=c(0,0,0,0),xaxs="i", yaxs="i") > plot(ecdf(datobs)) > lines(ecdf(gam_sum_gen)) > > Then I try to add colors etc and also the legend but fail. > > par(mar=c(4,4,2,1.2),oma=c(0,0,0,0),xaxs="i", yaxs="i") > plot(ecdf(datobs),main ="CDF of the sum for winter > season-Hume",cex.axis=1.2,xlab="Rainfall (mm)", pch = > pch,verticals=TRUE,col=c("black","red"), lty=c(1,1),ylab="Cumulative > probability", xlim=c(0,800),lwd=1) > lines(ecdf(gam_sum_gen)) > legend("topright", legend = c("observed","fitted"), >col = c("black","red"), pch=c(NA,NA), lty = c(1, 1), >lwd=c(3,3),bty="n", pt.cex=2) > box() > > I'm not sure why it doesn't show up at all. > > Thank you for any help given. > > > >[[alternative HTML version deleted]] > > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > > [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
[R] overlap cdf plots and add colors and etc
Hi r-users, I would like to overlap 2 ecdf plots. I tried this below and it gives me two plots of ecdf but just both just in black. par(mar=c(4,4,2,1.2),oma=c(0,0,0,0),xaxs="i", yaxs="i") plot(ecdf(datobs)) lines(ecdf(gam_sum_gen)) Then I try to add colors etc and also the legend but fail. par(mar=c(4,4,2,1.2),oma=c(0,0,0,0),xaxs="i", yaxs="i") plot(ecdf(datobs),main ="CDF of the sum for winter season-Hume",cex.axis=1.2,xlab="Rainfall (mm)", pch = pch,verticals=TRUE,col=c("black","red"), lty=c(1,1),ylab="Cumulative probability", xlim=c(0,800),lwd=1) lines(ecdf(gam_sum_gen)) legend("topright", legend = c("observed","fitted"), col = c("black","red"), pch=c(NA,NA), lty = c(1, 1), lwd=c(3,3),bty="n", pt.cex=2) box() I'm not sure why it doesn't show up at all. Thank you for any help given. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.