我在R中运行如下代码:
Abu=read.table('family_nor.out',header=T)Abu<-as.matrix(Abu)#Min support filtering for each taxaAbu[Abu<0.01]<-0#Filter taxa by abundance sumAbu<-Abu[which(rowSums(Abu)>10),]#Filter taxa by Occupancytable<-Abutable[table>0]<-1table1<-Abu[wh...