Increase margin of combined plot made out of subplots with Plots.jl

Specific Domains Visualization
1

I am plotting multiple plots in a combined plot like this:

using Plots
using Plots.Measures: mm

plots = Vector{Any}()
for (i,j) in Iterators.product(1:4, 1:4)
    p = scatter(rand(100), xlabel="X-axis $i,$j", ylabel="Y-axis $i,$j", title="Plot $i,$j")
    push!(plots, p)
end
plot(plots..., layout=(4,4), size = (2480, 3508) .÷ 4)

image
image867×1227 334 KB

However as can be seen the labels for the axes are cut off.

Therefore I would like to increase the margin on the left.

However when I use the margins keyword (or similar) like this:

plot(plots..., layout=(4,4),size = (2480, 3508) .÷ 4, left_margin=15mm, bottom_margin=15mm)

it will increase the margin of all subplots:

image
image867×1227 192 KB

I actually intend to change the margin of the total, outside plot. So how could I fix this?