如何将估算器（Chao1，Chao2，ACE和ICE）添加到物种累积图中

Question

我有以下脚本显示森林环境的积累曲线图，我很想整合Chao1，Chao2，ACE和ICE指标图（使用＆＃39;化石＆＃39; packe function 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)

像这张照片。