Hi,
ggpubr is a good solution :
library("ggpubr")
manchots = palmerpenguins::penguins
var_num = names(manchots)[sapply(manchots, is.numeric)]
graphiques = as.list(numeric(6))
k = 0
for (i in seq_along(var_num)) {
for (j in seq_len(i-1)) {
k = k+1
graphiques[[k]] = ggplot(manchots, aes_string(var_num[i], var_num[j], col = "species")) + geom_point()
}
}
gg = ggarrange(plotlist = graphiques, ncol = 2, nrow = 1, common.legend = TRUE)
ggexport(filename = "test.pdf", plotlist = gg, width = 12, height = 6)