In [1]:

from plotnine import *
from plotnine.composition import Stack
from plotnine.data import huron

In [2]:

p1 = (
    ggplot(huron, aes("year", "level"))
    + geom_line()
)

p2 = (
    ggplot(huron, aes("factor(decade)", "level"))
    + geom_boxplot()
    + labs(x="decade")
)

p3 = (
    ggplot(huron, aes("level"))
    + geom_histogram(bins=12)
    + scale_y_continuous(expand=(0.015, 0))
    + coord_flip()
)

p1 / p2 / p3

Which is equivalent to

In [3]:

Stack([p1, p2, p3])