subplot() doesn't account for dual axes

I think there is something wrong with the subplot function when there is yaxis2

I tried to use subplot with 2 plots that have yaxis2. Each plot looks like this:

foo_plot <- function(n) {
  x <- rnorm(n)
  den <- density(x)
  p <- plot_ly()
  p <- add_trace(p, x=x, type="histogram", name="Histogram",
    marker=list(line=list(color="#000", width=1), color="#555"))
  p <- add_trace(p, x=den$x, y=den$y, yaxis="y2", name="Density",
    type="scatter", mode="lines", line=list(color="#222"))
  p <- layout(p,
    yaxis2=list(overlaying="y", side="right", rangemode="nonnegative",
      showgrid=F, showline=T, ticks="outside"),
    yaxis=list(showline=T, ticks="outside"))
  return(p)
}

p <- subplot(foo_plot(1000), foo_plot(1000), nrows=2, shareX=F, shareY=F)

The density curves were missing when I merge the two plots.

This is the result of sessionInfo()

Start by running the issue's foo_plot and subplot example, comparing each plot with the combined two-row result. Trace subplot's handling of yaxis2 and verify that both density curves remain visible after merging; done means dual-axis traces render correctly in the combined plot.