Вы можете использовать par(mgp = plot_coords)
, как указано в этого сообщения. Например:
par(mfrow = c(2,2))
plot(1:10, col = 1, main = 1)
plot(1:10, col = 2, main = 2)
plot(1:10, col = 3, main = 3)
plot(1:10, col = 4, main = 4)
par(mfg = c(1,2))
points(10:1, col = "hotpink")
![enter image description here](https://i.stack.imgur.com/bpniZ.png)