R plot: draw polygon with coordinates based on second y-axis values

Question

I have a plot like the following:

plot(x, y)
par(new=T)
plot(x, z, xlab=NA, ylab=NA, xaxt="n", yaxt="n")
axis(side=4, at = pretty(range(z)))
<