design <- read.csv("design2.csv",header = T,row.names = 1,sep = ",")
bray_curtis <- read.csv("bray_curtis2.csv",sep = ",",header = T,check.names = F)
row.names(bray_curtis) <- bray_curtis[,1]
bray_curtis <- bray_curtis[,-1]
idx = rownames(design) %in% colnames(bray_curtis)
sub_design = design[idx,]
bray_curtis = bray_curtis[rownames(sub_design), rownames(sub_design)]
pcoa = cmdscale(bray_curtis, k=3, eig=T)
bray <- vegdist(data1, method = "bray")
警告消息:1:在 vegdist(data1, method = "bray"):你有空行:它们的相异性在方法 "bray" 中可能毫无意义 2: 在 vegdist(data1, method = "bray") 中: 结果中缺少值