Включить вывод лиммы plotMDS() в grid.arrange

Question

Я пытаюсь включить два отдельных вывода из plotMDS() в объект grid.arrange, но получаю ошибки и предупреждения следующим образом:

Error in gList(list(dim.plot = 1:2, distance.matrix = c(0, 3.88481791045152,  : 
  only 'grobs' allowed in "gList"
In addition: Warning messages:
1: In class(grob) <- c("gTableChild", class(grob)) :
  Setting class(x) to multiple strings ("gTableChild", "MDS", ...); result will no longer be an S4 object
2: In class(grob) <- c("gTableChild", class(grob)) :
  Setting class(x) to multiple strings ("gTableChild", "MDS", ...); result will no longer be an S4 object

MWE показан ниже:

library(gridExtra)
library(limma)

colors <- c(rep(c("purple"), 2), rep(c("blue"), 2))
data = data.frame(A = runif(100, 1, 10), B = runif(100, 1, 10), C = runif(100, 1, 10), D = runif(100, 1, 10))
data2 = data.frame(E = runif(100, 1, 10), F = runif(100, 1, 10), G = runif(100, 1, 10), H = runif(100, 1, 10))

plotMDS(data, col=colors, xlab="Dim 1", ylab="Dim 2", labels=colnames(data))
plotMDS(data2, col=colors, xlab="Dim 1", ylab="Dim 2", labels=colnames(data2))

grid.arrange(plotMDS(data, col=colors, xlab="Dim 1", ylab="Dim 2", labels=colnames(data)), plotMDS(data2, col=colors, xlab="Dim 1", ylab="Dim 2", labels=colnames(data2)), ncol = 2)

Я видел похожие посты без ответа ( https://www.biostars.org/p/163572/), и мне было любопытно узнать, есть ли у кого-нибудь совет по этому вопросу. Спасибо!