synthetic: Calculation of synthetic assemblages

Description

This fucntion calculates the ideal pollen assemblages associated with a series of values of a given environmental factor.

Usage

synthetic(x, y, z, delta, plot = FALSE, percen, upmargin = 0.1, widths = c(1, 4), steps = 0.05)

Arguments

A vector containing the environmental gradient value for each sample or location.

A matrix containing the observed taxa abundances along the environmental gradient. Taxa are in columns and samples in rows. Taxa included in this analysis should be selected by evaluating dispersion and residuals distribution through the function dispersion.all.

A matrix containing the parameters for the LOESS regression of each species ($\alpha$ and degree in that order). Objects produced by bcrossv.all and akaike.all could be used as z.

delta

Increments for the LOESS regression prediction. See predict.loess for details.

plot

Logical indicating whether or not a plot of the synthetic assemblages is desired.

percen

A vector containing total percentage of pollen in each modern sample. As some taxa might have been excluded from the pollen sum, this value could be larger than 100. See details below.

upmargin

A number indicating the margin size for the upper part of each panel. par function, mai parameter.

widths

A vector of values for the widths of columns in the plotting device. layout parameter.

steps

The proportion of the range to draw lines separating each synthetic assemblage.

Value

loess.f: A list with the LOESS regression for each taxon.
predicted: A matrix containing the predicted percentages of each taxon along the evaluated environmetal gradient. The first column contains the values of the environmental variable, and taxa are contained in the subsequent columns.
avg.info: if plot=TRUE, this value is the average pollen percentage used from each sample in the construction of the synthetic assemblages.

Details

When applying the synthetic assemblages technique, it is desirable to know the proportion of data from each sample that was used in the the construction of the modern ideal assemblages. The avg.info returns the proportion of pollen data that was used from each sample. If no taxa are excluded from the pollen sum, such proportion will coincide with the sum of the percentages of the taxa included. In this case, the argument percent can be filled by repeating 100 as many times as number of of modern samples (percent = rep(100, nrow(y))). When one or more taxa are excluded from the pollen sum, samples might have total percentages that are larger than 100. Therefore, the percent has to be calculated by adding all the percentages in each sample (percent = apply (Y, 1, sum), where Y is the original percentage matrix including all taxa).

References

Correa-Metrio, A., M.B. Bush, L.Perez, A. Schwalb, and K. Cabrera. Accepted. Pollen distribution along climatic and biogeographic gradients in northern Central America. The Holocene

Examples

Run this code

data(modernq)
# Calculate percentages
perq<-percenta(modernq,first=2,last=39)[,2:55]
# filter data set to include only samples with at least 0.5
# percent in 20 percent of the samples
perq1<-filter.p(perq,presen=0.5,persist=0.2)$filtered
# calculate alpha and degree for each taxon through AIC
a.d<-akaike.all(modernq[,1],perq1)
# select taxa that have acceptable dispersion and normally 
# distributed residuals in percentages and AIC matrices
perq2<-perq1[,-c(3:5,8,17,19)]
a.d1<-a.d[-c(3:5,8,17,19),]
#build synthetic assemblages
syas.q<-synthetic(modernq[,1],perq2,a.d1,delta=25,plot=TRUE,
percen=rep(100,53))
# predic percentages of Cyperacea for an elevation sequence from
# 100 to 550 in 25-m increments
predict(syas.q[[1]]$Cyperaceae,seq(100,550,25))

Run the code above in your browser using DataLab

Get 50% off unlimited learning