rarc
computes rarefaction curves and returns a data frame with bootstrap estimates of species richness, its statistical envelop and the average number of individuals for a given sample size.
1 
matrix 

samplesize 
A vector containing the different sample sizes for which the computations are required. If missing or non vector, 
nrandom 

p1 
The probability value used to compute the upper bound of the statistical envelop using 
p2 
The probability value used to compute the lower bound of the statistical envelop using 
save 
If 
out 
A data frame with 5 columns giving the mean species richness, the lower and upper bound of the richness corresponding to the quantile values for p1 and p2, the mean number of individuals and the sample size. 
bootstrapped.val 
Generated if 
JeanPierre Rossi jeanpierre.rossi@supagro.inra.fr
1 2 3 4 5 6 7 8 9 10 11 12 13 14  ## Not run:
data(ef)
t < rarc(ef,samplesize=c(5,10,15,20,25), nrandom=30, p1=0.975, p2=0.025)
head(t)
# Plot the outputs
plot(t$out[,6],t$out[,1], type="b", ylim=range(c(t$out[,2],t$out[,3])),
xlab="number of sampling units", ylab="richness")
points(t$out[,6] , t$out[,2], type="l", col="red")
points(t$out[,6] , t$out[,3], type="l", col="blue")
## End(Not run)
#See the package vignette for more examples. Type: vignette("rich_introduction").

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.
Please suggest features or report bugs with the GitHub issue tracker.
All documentation is copyright its authors; we didn't write any of that.