通过R语言picante包求NTI

[Arcgis小白学习经验分享]

求问各位前辈大佬，以下是我想通过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

[Janet8]

@heroy26 请问植物群落数据一般选择哪种零模型更好，样方内植物生物量（干重）

[heroy26]

null.model = c("taxa.labels", "richness", "frequency", "sample.pool", "phylogeny.pool", "independentswap", "trialswap")
这里面是7种0模型，选择一种就可以，比如null.model = "independentswap"

[weixin_53592269]

@Janet8 应该是你的树上的物种名与表格里的物种名没有对应上

[Janet8]

@RickyKing666 好像没用

[RickyKing666]

@Janet8 丰度表，转置一下样品名和OTU名的行列

[Janet8]

请问下标出界是怎么回事

[m0_61633985]

您好，想请教一下怎么计算NTI指数，可以加下qq吗 996676380