R 如何向物种累积图中添加估计器（Chao1、Chao2、ACE和ICE）

我有下面的脚本，显示了森林环境的累积曲线图，我很想在同一累积曲线中集成Chao1、Chao2、ACE和ICE指标图（使用“化石”packe函数spp.est（x）生成的矩阵），我需要做什么来包括它们

    x = structure(list(
        Achromoporus.occultus = c(0L, 0L, 0L, 0L, 0L), 
        Achromoporus.platyurus = c(0L, 0L, 0L, 0L, 0L), 
        Achromoporus.andujari = c(0L, 0L, 0L, 0L, 0L), 
        Achromoporus.occultus.1 = c(0L, 0L, 0L, 0L, 0L),
        Achromoporus.platyurus.1 = c(0L, 0L, 0L, 0L, 0L), 
        Achromoporus.vallenuevo = c(0L, 0L, 0L, 0L, 0L),
        Acromoporus.Andujari = c(4L, 3L, 1L, 5L, 0L),
        Acromoporus.concolor = c(0L, 0L, 0L, 0L, 1L),
        Acromoporus.occultus = c(2L, 3L, 10L, 16L, 1L),
        Acromoporus.Platyurus = c(1L, 6L, 3L, 15L, 3L),
        Acromoporus.Vallenuevo = c(4L, 0L, 0L, 2L, 0L),
        Docodesmus.angustus = c(0L, 0L, 0L, 1L, 0L),
        Docodesmus.Sp. = c(0L, 1L, 0L, 0L, 0L),
        Henicomus.Septiporus = c(0L, 0L, 0L, 0L, 0L),
        Morfo.A = c(0L, 1L, 1L, 0L, 0L),
        Morfo.B = c(0L, 0L, 0L, 0L, 0L),
        Prostemmiulus.gracilipes = c(0L, 0L, 0L, 0L, 0L),
        Prostemmiulus.Scaurus = c(0L, 0L, 4L, 3L, 5L),
        Prostemmiulus.setosus = c(1L, 5L, 3L, 6L, 6L),
        Siphonophora.platops = c(0L, 0L, 0L, 0L, 0L),
        Spirobolus.sp = c(0L, 0L, 0L, 0L, 0L),
        Spirobolus.sp.Morfo.A = c(0L, 0L, 0L, 0L, 0L),
        Syphonophora.Platops = c(0L, 0L, 1L, 0L, 0L)), 
        .Names = c(
            "Achromoporus.occultus",
            "Achromoporus.platyurus", 
            "Achromoporus.andujari",
            "Achromoporus.occultus.1",
            "Achromoporus.platyurus.1", 
             "Achromoporus.vallenuevo",
            "Acromoporus.Andujari",
            "Acromoporus.concolor", 
            "Acromoporus.occultus",
            "Acromoporus.Platyurus",
            "Acromoporus.Vallenuevo", 
            "Docodesmus.angustus",
            "Docodesmus.Sp.",
            "Henicomus.Septiporus", 
            "Morfo.A",
            "Morfo.B",
            "Prostemmiulus.gracilipes",
            "Prostemmiulus.Scaurus", 
            "Prostemmiulus.setosus",
            "Siphonophora.platops",
            "Spirobolus.sp", 
            "Spirobolus.sp.Morfo.A",
            "Syphonophora.Platops"), 
             row.names = c(NA, 5L), class = "data.frame")


    x = read.csv("Desktop/pine.csv")
    plot(specaccum(x), marker='.')
    legend('bottomright', c('Pine Forest'),
           col=1:3, lty=1, bty='n', inset=0.025)

---

类似于这张照片的东西。

不清楚您想要什么。因为您没有包含任何示例数据，所以我们无法运行代码或查看输出。这个代码产生了什么？你想要的和生产的有什么不同？是否有您试图在某处创建的绘图类型的示例？可能您希望使用基于此处所述程序的

素食者

包装和绘图曲线：（我自己不使用这些方法进行物种丰富度估计。）上述论文中所述的方法也可能在此处可用。不过，我并不是通过提供可能对它们进行评估的软件名称来支持它们的使用。@MrFlick哪种方式最能让您了解我的样本数据？有关建议，请参阅。