将ecdf和密度绘制到具有不同y轴的同一图形中

Question

Consider the two following plots:

x <- rnorm(1000, sd=10)
plot(density(x))
plot(ecdf(x))

Please note the different y-axis ranges:

Now I would like to have these two in one graph and two y-axis.

I know the basics ( add=TRUE , or par(new=TRUE) ), but add=TRUE does not work due to the different scales, par(new=TRUE) requires specifying a x-axis to avoid over plotting.

Is there an easy way of getting these two into one graph, left y-axis density, right y-axis ecdf?

Answer 1

Create plots with two different Y scales is not generally recommended, but the easiest way in this case would be

    x <- rnorm(1000, sd=10)
    xlim=c(-50, 40)
    par(mar = c(5,4,4,5))
    plot(density(x), xlim=xlim)
    par(new=TRUE)
    plot(ecdf(x), axes=F, xlab="", ylab="", main="", xlim=xlim)
    axis(4)
    mtext("ECDF",4, 3)

If you have a comelling reason not to want to set identical xlim values, you can borrow the answer from this question to extract the xlim from the previous plot

getplotlim<-function() {
    usr <- par('usr')
     xr <- (usr[2] - usr[1]) / 27 # 27 = (100 + 2*4) / 4
     yr <- (usr[4] - usr[3]) / 27
     list(
         xlim = c(usr[1] + xr, usr[2] - xr),
         ylim = c(usr[3] + yr, usr[4] - yr)
     )
}


x <- rnorm(1000, sd=10)
par(mar = c(5,4,4,5))
plot(density(x))
par(new=TRUE)
plot(ecdf(x), axes=F, xlab="", ylab="", main="", xlim=getplotlim()$xlim)
axis(4)
mtext("ECDF",4, 3)