preseqR.interpolate.rSAC: Interpolation

Description Usage Arguments Details Value Author(s) References Examples

View source: R/rSAC.R

Description

Interpolating the number of species represented at least r times in a subsample given an initial sample

Usage

1

Arguments

n

A two-column matrix. The first column is the frequency j = 1,2,…; and the second column is N_j, the number of species with each species represented exactly j times in the initial sample. The first column must be sorted in an ascending order.

ss

A positive double equal to the step size between subsamples.

r

A positive integer. Default is 1

Details

The expected number of species represented at least r times in the subsample is estimated based on an expended version of the formula by Heck Jr, KL. et al. (1975).

Value

A two-column matrix for the number of species represented at least r times in a random sample. The first column is the size of the random sample; the second column is the expected number of species represented at least r times in the sample.

NULL if failed.

Author(s)

Chao Deng

References

Heck Jr, K. L., van Belle, G., & Simberloff, D. (1975). Explicit calculation of the rarefaction diversity measurement and the determination of sufficient sample size. Ecology, 1459-1461.

Examples

1
2
3
4
5
6
7
8
9
## load library
library(preseqR)

## import data
data(Shakespeare)

## The expected number of distinct words represented twice or more in the
## subsample 
preseqR.interpolate.rSAC(n=Shakespeare, ss=1e5, r=2)

Example output

     sample.size interpolation
[1,]       1e+05      4923.340
[2,]       2e+05      7595.853
[3,]       3e+05      9620.429
[4,]       4e+05     11298.566
[5,]       5e+05     12752.162
[6,]       6e+05     14046.429
[7,]       7e+05     15221.208
[8,]       8e+05     16302.507
There were 50 or more warnings (use warnings() to see the first 50)

preseqR documentation built on May 2, 2019, 6:39 a.m.