paleoMAS (version 2.0-1)

simulat.t: Model desaturation by samples simulation

Description

This function simulates new samples for a sequence of fossil data to desaturate models.

Usage

simulat.t(x, pop = 1000, nsamples = 10, ssample = 300, percenta = TRUE, last)

Arguments

x

A percentage matrix containing samples in rows and taxa in columns.

pop

Size of the population to be simulated.

nsamples

Number of samples to be drawn.

ssample

Sample size.

percenta

Logical indicating whether or not the outcome matrix should be one of percentages.

last

If percenta=TRUE, last refers to the column where the last species to be taken into account is located. Only relevant if percenta=TRUE. See percenta for details.

Value

When percenta=FALSE, samples.t is returned. When percenta=TRUE an additional element percentages is returned.

samples.t

An absolute abundance matrix with samples in rows and species in columns. The first column contains sample number.

percentages

A relative abundance matrix with samples in rows and species in columns. Each species is repeated nsamples times.

Details

When all species are taken into the pollen count, last corresponds to the number of columns. If only some species are to be taken into the pollen count, they must be contiguous and start in the first column.

References

Correa-Metrio, A., K.R. Cabrera, and M.B. Bush. 2010. Quantifying ecological change through discriminant analysis: a paleoecological example from the Peruvian Amazon. Journal of Vegetation Science 21: 695-704.

See Also

This function is an iteration of simulat

Examples

data(quexilper)
colnames(quexilper)
# Cyperaceae and spores are left out of the pollen sum.
# For absolute abundance data
simulat.t(quexilper,percenta=FALSE,last=73)
# or
simulat.t(quexilper,last=73)$samples.t
# For percentages data
simulat.t(quexilper,last=73)$percentages