求問各位前輩大佬,以下是我想通過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
轉載請註明出處,本文鏈接:https://www.uj5u.com/qita/186522.html
標籤:其他開發語言
