[R] Two ecdf with log-scales
Johannes Radinger
JRadinger at gmx.at
Thu May 3 18:32:09 CEST 2012
Thank you Steve,
thats the thing I was looking for....
/Johannes
-------- Original-Nachricht --------
> Datum: Thu, 3 May 2012 08:20:51 -0400
> Von: "Steven Wolf" <wolfste4 at msu.edu>
> An: "\'David Winsemius\'" <dwinsemius at comcast.net>, "\'Johannes Radinger\'" <JRadinger at gmx.at>
> CC: R-help at r-project.org
> Betreff: RE: [R] Two ecdf with log-scales
> I've done it this way before:
>
> eX -> ecdf(distribution 1)
> eY -> ecdf(distribution 2)
> par(mar=c(5,5,2,1),xlog=TRUE)
> plot(eX, do.points=FALSE, verticals=TRUE, col="black", xlab=xlabel,
> xlim=c(1,100000), ylab=ylabel,
> lty=1, cex.lab=1.5, cex.axis=1.5, main="",
> lwd=3,log="x")
> plot(eY, do.points=FALSE, verticals=TRUE, col="blue", add=TRUE,
> xlim=c(1,100000), main="")
>
> Warning: It makes a stair-step that may be difficult to see unless you
> use
> color. I had to change how the ecdf was plotted when I made b/w figures
> for
> my publication so that different dashed lines were distinct.
>
> HTH,
> -Steve
>
> -----Original Message-----
> From: David Winsemius [mailto:dwinsemius at comcast.net]
> Sent: Wednesday, May 02, 2012 10:17 AM
> To: Johannes Radinger
> Cc: R-help at r-project.org
> Subject: Re: [R] Two ecdf with log-scales
>
>
> On May 2, 2012, at 6:14 AM, Johannes Radinger wrote:
>
> > Hi,
> >
> > i want to plot empirical cumulative density functions for two
> > variables in one plot. For better visualizing the differences in the
> > two cumulative curves I'd like to log-scale the axis.
> >
> > So far I found 3 possible functions to plot ecdf:
> >
> > 1) ecdf() from the package 'stats'. I don't know how to successfully
> > set the log.scales? Combining two plots is not a problem:
> >
> > plot(ecdf(x1))
> > lines(ecdf(x2),col.h="red")
> >
> > 2) gx.ecdf() from package 'rgr'. It is easily possible to plot log-
> > scales, but I don't know how to plot two densities?
> >
> > gx.ecdf(x1,log=TRUE,ifqs = TRUE)
> >
> > 3) Ecdf() from package 'Hmisc'. No log-option directly available and
> > here I also don't know how to 'stack' two plots...
> >
> > Ecdf(x1,what="F")
> >
> >
> > Probably there are many more solutions (e.g. ggplot etc.)...
> > ...Has anyone faced a similar task and found a simple solution? Any
> > suggestions are welcome!
>
> Have you searched the Archives? I seem to remember that the log(0) was a
> barrier to persons attempting this in the past. (ISTR a posting in the
> last
> few weeks.) Maybe you could also provide a test data object that has the
> same range as your x1 and x 2 variables.
>
> > and provide commented, minimal, self-contained, reproducible code.
>
> --
>
> David Winsemius, MD
> West Hartford, CT
>
>
>
--
Jetzt informieren: http://mobile.1und1.de/?ac=OM.PW.PW003K20328T7073a
More information about the R-help
mailing list