beta多样性计算

file=c("salmonRNA.csv","salmonDNA.csv")
for (j in file)
{
  c<-read.csv(j,row.names=1,stringsAsFactors=FALSE)
  vegan::vegdist(t(c),method="bray") %>% ape::pcoa()->pcoa
  tmp<-if (is.null( pcoa$values$Rel_corr_eig))
    pcoa$values$Relative_eig else pcoa$values$Rel_corr_eig
  write.csv(tmp,file=paste(j,"Axexp.csv"))
  write.csv(pcoa$vectors,file=paste(j,"Ax.csv"))
}