老师您好,我根据您提供的物种累积曲线的代码输入后绘制不出结果,以下是我的输出结果而没有得到图:
require(vegan)
> otu <- read.table("C:/Users/CHEN/Desktop/otu_table.txt", header = T,
+ sep = "\t",row.names=1)
> all <- specaccum(otu, method = "random")
> png(file="spe.png",h=2000,w=2000,res=300)
> plot(all, ci.type = "poly", col = "blue", lwd = 2, ci.lty = 0,
+ ci.col = "lightblue", main = "Species Accumulation Curves", xlab = "Number of samples",
+ ylab = "Number of species")
> boxplot(all, col = "yellow", add = TRUE, pch = "+")
> dev.off()
RStudioGD
2