如何使用素食主义者在 R 中的 plot 稀疏曲线？

Question

我有一个数据集，其中包含在不同参考基因组中鉴定的基因。 因此，参考基因组在行中，基因在表的列中。 该表被编码为二进制，其中0表示基因不存在， 1表示基因存在。 我做了基因积累曲线，这表明每个基因组的基因数量正在接近一个平台期。 现在，我正在尝试使用 R-package vegan对稀疏曲线进行 plot 。 我使用了以下代码：

b<-read.csv("data.csv", header = T, check.names = F)
S <- specnumber(b) # observed number of species
(raremax <- min(rowSums(b)))
Srare <- rarefy(b, raremax)
plot(Srare, xlab = "Observed No. of genes", ylab = "Rarefied No. of genes")
abline(0, 1)
rarecurve(b, step = 15, sample = raremax, col = "blue", cex = 0.6)

数据集如下：

          gene1 gene2 gene3
#genome1    0     1     0
#genome2    1     0     1
#genome3    1     0     1

但是，使用此代码我没有得到任何令人满意的 output。 我只得到一条通过对角线的直线。 我在下面附上了 output。

有人可以建议我如何纠正 output？

谢谢你。

Answer 1

稀有rarefy稀有数据的各个行：它对每行中出现的事件（“个人”）进行子样本。 如果所有这些抽样个体的值都是 1，那么您将有一个子样本，而一个子样本的总和就是样本量：这就是您得到的。 稀疏向量没有有意义的方法：您需要一些计数 > 1 的计数数据。

在对矩阵的行进行二次采样时，您可能正在寻找整个数据集中的基因积累。 这是在素食主义者function specaccum （参数method = "exact" ）中完成的，它有自己的plot等方法。