fossil.dist: Distance calculation between fossil samples and modern synthetic assemblages

Description

This function calculates the distance between each fossil sample and the synthetic assemblages, according to the selected distance method.

Usage

fossil.dist(modern, fossil, method = "canberra")

Arguments

modern

An object containing the synthetic assemblages produced by the function synthetic.

fossil

A matrix containing the fossil pollen percentages with samples in rows and taxa in columns.

method

Corresponds to methods available in vegdist of vegan. Available methods are "manhattan", "euclidean", "canberra", "bray", "kulczynski", "jaccard", "gower", "morisita", "horn", "mountford", "raup" , "binomial" or "chao". See vegdist (package vegan) for details.

Value

A matrix containing the distance or dissimilarity between fossil samples and each one of the modern synthetic assemblages. Fossil samples are in rows and modern synthetic assemblages are in columns. The first row contains the values of the environmental gradient associated with each synthetic assemblage.

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,fossilq)
# 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))
#calculate distance from synthetic assemblages to each fossil sample
fossil.dist(syas.q,fossilq)