用标准误差绘制的rarecurve()

时间:2015-06-04 23:08:38

标签: r vegan

rarecurve()(menuButton.addTarget(self.revealViewController(), action: "revealToggle:", forControlEvents: UIControlEvents.TouchUpInside) )是否接受绘图的标准错误?

如果是这样,我该如何绘制这样一条曲线?

我正在遵循一个经典的脚本,使用BCI数据集:

vegan

从统计学角度讲,促进这一功能对大多数S <- specnumber(BCI) (raremax <- min(rowSums(BCI))) Srare <- rarefy(BCI, raremax) plot(S, Srare, xlab = "Observed No. of Species", ylab = "Rarefied No. of Species") abline(0, 1) rarecurve(BCI, step = 20, sample = raremax, col = "blue", cex = 0.6) 用户都有帮助。

谢谢! 安德烈

1 个答案:

答案 0 :(得分:1)

rarecurve没有给你SE。原因很明显,已经给你了:没有额外的曲线就有足够的杂乱。如果你真的想这样做,你必须手动完成。这并不太复杂,因为rarefy函数接受矢量样本大小并为您提供所需的所有数字。下面使用Barro Colorado数据集的一个站点绘制基本图:

library(vegan)
data(BCI)
sum(BCI[1,]) # site 1, 448 tree stems
N <- seq(2, 448, by=8)
S <- rarefy(BCI[1,], N, se = TRUE)
plot(N, S[1,], type="l", lwd=3)
lines(N, S[1,] + 2*S[2,]) ## 2*SE is good enough for 95% CI
lines(N, S[1,] - 2*S[2,])

从统计学上讲,这只给出了由子采样过程引起的误差,假设观察到的数据没有随机变化。对我来说这没有多大意义,我发现稀有的SE是误导性和无意义的。这并不能阻止我在素食主义者中提供它们。