根据R, OTU表不匹配元数据表(头)



我已经做了一个多星期了,它仍然不起作用。我在r中使用包"RAM",更具体地说是OTU。多样性功能:

*总结OTU表的多样性指数

描述

这些函数计算所有样本的多样性指数,并将输出作为新列追加到元数据表中。

使用

OTU。多样性(元数据)参数一个OTU表的列表。Meta元数据以附加输出。

详细信息

该函数总结了给定otu表的以下多样性指数:specnumber、shannon、simpson、invsimpson、truediversity、偶数、chao和ACE指数。

价值

该函数返回每个样本的多样性指数向量,这些向量被附加到给定的元数据表中。*

我的命令行是:

tabletemp <- read.table(paste(pathDataAn, "ITS.table.tsv", sep = ""),
                  sep = "t", header = TRUE, dec = ".", comment.char = "", quote = "", stringsAsFactors = TRUE,
                  as.is = TRUE, check.names = FALSE, colClasses=c("row.names"="character", "taxonomy"="character"))
tabletemp$row.names <- NULL
metatemp <-  read.table(paste(pathDataAn, "ITS.meta.tsv", sep = ""),
                  sep = "t", header = TRUE, dec = ".", comment.char = "", quote = "", stringsAsFactors = TRUE,
                  as.is = TRUE)
temp2 <- OTU.diversity(list(data=tabletemp), metatemp)

,输出是这样的:

.valid错误。Meta (otu, Meta = Meta):otu1和meta的样例名称不同。它们可能出了问题。在otu1,而不是在meta:Em13s01rv em13s02rv em13s03fw em13s03rv…(所有样品的清单)在meta中,不在otu1中:

1,2,3,…

我已经验证了每个标头是相同的,并且顺序完全相同。它们都是正确排列的。请帮助我,我不明白为什么它似乎不能正确读取这些表

在周末,我认为我需要告诉软件列标头(OTU表)匹配行名(元数据表)

- colnames(tabletemp)[-ncol(tabletemp)]