如何解决Phyloseq sample_data和sample_names无法将样本合并到phyloseq对象中
我的问题与使用phyloseq()建立一个phyloseq对象有关。我遵循了https://joey711.github.io/phyloseq/import-data.html和https://vaulot.github.io/tutorials/Phyloseq_tutorial.html ,并且该代码以前已经起作用。
具体来说,当我的样本在otu_mat和samples_df中都为SAME时,为什么sample_names()返回一个NULL值向量
以下是开始文件:
> head(otu_mat)
Species 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
1 Absiella dolichum 44 14 3 152 50 106 4 4 14 41 2 0 1 21 1 2 4 2 2 4 0 4 7 0
2 Acetivibrio ethanolgignens 50 4 10 65 30 96 71 88 35 65 16 21 108 44 12 40 8 64 42 77 29 24 21 31
3 Acetoanaerobium sticklandii 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
4 Acetobacterium bakii 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
5 Acholeplasma axanthum 53 26 2 45 9 19 8 7 4 36 3 4 75 16 2 18 6 22 4 2 25 13 3 13
6 Acholeplasma brassicae 0 0 0 0 1 0 0 0 0 1 0 0 2 0 0 2 1 1 1 1 0 3 0 3
> head(tax_mat)
Kingdom Phylum Class Order Family Genus Species
<chr> <chr> <chr> <chr> <chr> <chr> <chr>
1 Bacteria Firmicutes Clostrid… Clostridi… Lachnospiraceae unclassified Lachnospir… Lachnospiraceae ba…
2 Bacteria Bacteroid… Bacteroi… Bacteroid… Tannerellaceae Parabacteroides unclassified Parab…
3 Bacteria Firmicutes Clostrid… Clostridi… Lachnospiraceae Blautia Ruminococcus gnavus
4 Bacteria Firmicutes Clostrid… Clostridi… Ruminococcaceae Ruminococcus Ruminococcus bromii
5 Bacteria Firmicutes Clostrid… Clostridi… Clostridiales inc… Clostridiales Family XI… Emergencia
6 Bacteria Firmicutes Clostrid… Clostridi… Ruminococcaceae Anaerotruncus unclassified Anaer…
> head(samples_df)
OG_Sample Sample Activity Diet Sex Time Weight Treatment
1 101 1 Sedentary Control Female Week_3 21 Sedentary_Control_Female
2 102 2 Sedentary Control Female Week_3 20 Sedentary_Control_Female
3 103 3 Sedentary Control Female Week_3 21 Sedentary_Control_Female
4 104 4 Sedentary Control Female Week_3 22 Sedentary_Control_Female
5 105 5 Sedentary Control Female Week_3 19 Sedentary_Control_Female
6 106 6 Sedentary Control Female Week_3 19 Sedentary_Control_Female
# Here is the code:
otu_mat <- w3.long %>% rownames_to_column("Species")
# make tax_otu file (OTUs as rows,taxnomoy as columns)
tax_mat1 <- l %>%
select(Kingdom,Phylum,Class,Order,Family,Genus,Species) %>%
group_by(Species)
# This filters tax_mat1 by Species in otu_mat
tax_mat <- tax_mat1 %>%
filter(Species %in% otu_mat$Species) %>%
distinct(.keep_all = FALSE)
# Subset of Metadata file for conformable arrays
w3.env <- as.data.frame(GS2.env[1:56,])
samples_df <- w3.env %>% select(Sample,Activity,Diet,Sex,Time,Weight,Treatment)
rownames(samples_df) <- samples_df$Sample
# --------------------- Preprocessing for Phyloseq ------------------------------------------
row.names(otu_mat) <- otu_mat$Species
otu_mat <- otu_mat %>% select(-Species)
row.names(tax_mat) <- tax_mat$Species
otu_mat <- as.matrix(otu_mat)
tax_mat <- as.matrix(tax_mat)
OTU = otu_table(otu_mat,taxa_are_rows = TRUE)
TAX = tax_table(tax_mat)
samples = sample_data(samples_df)
# Phyloseq uses the phyloseq() command to rapidly combine the otu_table,# tax_table and samples
carbom <- phyloseq(OTU,TAX,samples)
carbom
这是错误:
validobject(.Object)中的错误:无效的类“ phyloseq”对象: 组分样品名称不匹配。 尝试sample_names()
我使用sample_names()重试,它返回一个NULL值的向量
samples = sample_names(samples_df)
在这里验证样品名称是否匹配
> # find the overlapping samples
> common.ids <- intersect(rownames(samples),colnames(otu_mat))
> dim(common.ids)
NULL
非常感谢您的帮助!
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。