RDA-PLS:多数据集关联分析=as.vector(rda_sum$cont$importance[2,1:2], mode="any")
rdapp=paste(rdap*100, "%",sep="")
#计算物种加权得分并排序 tscore=numeric(length(rownames(sp)))
for (i in 1:length(tscore)) {
tscore[i]=abs(rdap[1]*sp[i,1])+ abs(rdap[2]*sp[i,2])
}
spt=cbind(sp, tscore)
spt=spt[rev(order(spt[,3])),]
sph=spt[1:100,]
根据物种在两个约束排序轴的得分 =paste("RDA1(", rdapp[1], ")", sep="")
rdap2=paste("RDA2(", rdapp[2], ")", sep="")
limx=max(abs(si[,1 , ylab=rdap2, main="RDA analysis", cex=1, las=1)
pointLabel(x=cn[,1]/10, y=cn[,2]/10, labels=rownames 