simulat.t: Model desaturation by samples simulation
In paleoMAS: Paleoecological Analysis

Description Usage Arguments Details Value Author(s) References See Also Examples

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

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

`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.

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.

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.

Alexander Correa-Metrio, Kenneth R. Cabrera.

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.

This function is an iteration of simulat

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