Создание граненого графика с помощью функции phyloseq plot_composition

Question

У меня есть объект phyloseq с 2 факторами (время и лечение). Я сделал график относительной численности с этим кодом:

p1 <- plot_composition(pseqTSS$Rumen$phyTR.core, average_by = "Treat") +
      theme(legend.text = element_text(colour = "black", size = 7)) +
      theme(legend.key.size = unit(.3, "cm")) + 
      guides(fill = guide_legend(ncol = 3))

С этим выводом:

Когда я добавляю facet_wrap с другим фактором, как это,

p1 <- p1 + facet_wrap(eval(parse(text = paste0("~","Time"))))

ошибки вывода кода:

Error: At least one layer must contain all faceting variables: `Time`.
* Plot is missing `Time`
* Layer 1 is missing `Time`

Можно ли получить сюжет как первый, добавив facet_wrap с другим фактором?

Спасибо!