I would like to generate a matrix of boxplots that is the Cartesian product of two parameters. The attached is almost what I want but I would like to declutter it by factoring out the repeated “|V|=xx” into column headings and the “d=yy” into row titles.

I thought that by adding an extra row and column of empty Axis objects I could write the row and column titles into these but it does not appear to work. Could anybody advise please?

I’ll show you how to do it with Label even though you could technically use the upper axes’ titles for the columns:

f = Figure()
for i in 1:4, j in 1:3
Axis(f[i, j])
i == 1 && Label(f[i, j, Top()], "Column $j")
j == 3 && Label(f[i, j, Right()], "Row $i", rotation = pi/2)
end
f

You can add padding like padding = (5, 0, 0, 0) for 5 on the left side.

You could also make whole rows and columns for the labels, but that will lead to problems when you want to align your grid with other stuff next to it.