我目前在qiime2R 中遇到了一些麻烦。我正在关注教程并被困在这里:
unwunifrac <- read_qza("unweighted_unifrac_pcoa_results.qza")
unwunifrac <- read_qza("unweighted_unifrac_pcoa_results.qza")
unwunifrac$data$Vectors %>%
select(SampleID, PC1, PC2) %>%
left_join(metadata) %>%
left_join(shannon)
ggplot(aes(x=PC1, y=PC2, color=`SPT_No`, shape=`Geosmin`, size=shannon_entropy)) +
geom_point(alpha=0.5) +
theme_q2r() +
scale_shape_manual(values=c(16,1), name="Geosmin") +
scale_size_continuous(name="Geosmin") +
scale_color_discrete(name="Location")
这是因为我的 sampleid 在我的元数据和香农向量中称为 sampleid,但在未加权的 unifrac qza 文件中称为 SampleID。
有没有人遇到过这个问题?我该如何解决?最好重命名列吗?
任何帮助将非常感激。
非常感谢