0

我目前在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。

有没有人遇到过这个问题?我该如何解决?最好重命名列吗?

任何帮助将非常感激。

非常感谢

4

1 回答 1

0

你遇到的这个问题不清楚。您能否尝试发布一些内容以重现您从 q2R 收到的任何错误消息?

于 2021-05-12T21:44:13.957 回答