我已经做了一个多星期了,它仍然不起作用。我在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)]