通过R语言picante包求NTI。

Question

求问各位前辈大佬，以下是我想通过R语言的picante包的ses.mntd求NTI的代码，出现了如下所示的问题，请问该怎么解决呢？
谢谢各位。
library(picante)
mydata<-read.csv("C:\Users\PRCLUSHAN\Desktop\OTUas.csv", header=TRUE, row.names=1)
comm<-t(mydata)
tree<-read.tree("C:\Users\PRCLUSHAN\Desktop\Newick Export1.NWK")
phydist <- cophenetic(tree)
ses.mntd(comm, phydist, null.model = c("taxa.labels", "richness", "frequency", "sample.pool", "phylogeny.pool", "independentswap", "trialswap"), abundance.weighted=FALSE, runs = 999, iterations = 1000)

Error in data.frame(ntaxa = specnumber(samp), mntd.obs, mntd.rand.mean, :
'row.names' should specify one of the variables

Answer 1

可能是树里面的OTU名称和OTU表中的名称不完全一样，毕竟树是在序列分析中产生的，而OTU表会有抽平和出掉低丰度OTU的过程