У меня есть объект 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
с другим фактором?
Спасибо!