generalized_Tradidiss: Plot-to-plot dissimilarity taking account of functional...

In adiv: Analysis of Diversity

Plot-to-plot dissimilarity taking account of functional dissimilarities between species

Description

Given a matrix of S species' relative or absolute abundance values in N plots, together with an S x S (functional) dissimilarity matrix, the function generalized_Tradidiss calculates a semimatrix with the values of a plot-to-plot dissimilarity index, as proposed in Pavoine and Ricotta (2019).

Usage

generalized_Tradidiss(comm, dis, method = c("GC", "MS", "PE"),
    abundance = c("relative", "absolute", "none"), 
    weights = c("uneven", "even"), tol = 1e-8)

Arguments

comm	a data frame typically with communities as rows, species as columns and an index of abundance as entries. Species must be labeled as in the object dis.
dis	an object of class dist or a matrix containing pair-wise (functional or phylogenetic) dissimilarities between species [obtained by functions like vegdist in package vegan (Oksanen et al. 2013), gowdis in package FD (Laliberte and Shipley 2011), or dist.ktab in package ade4 for functional dissimilarities (Dray et al. 2007), or functions like cophenetic.phylo in package ape (Paradis et al. 2004) or distTips in package adephylo (Jombart and Dray 2010) for phylogenetic dissimilarities]. If the dissimilarities are outside the range 0-1, a warning message is displayed and each dissimilarity is divided by the maximum over all pairwise dissimilarities.
method	one of the following strings: "GC", "MS", "PE". See Details.
abundance	a string with three possible values: "relative" for the use of relative species abundance, "absolute" for the use of absolute species abundance, and "none" for the use of presence/absence data (1/0).
weights	a string. Two types of weights are available in the function: "uneven" (Eq. 5 in Pavoine and Ricotta (2009)) or "even" (1/S, where S is the number of species in the two compared plots).
tol	numeric tolerance threshold: values between -tol and tol are considered equal to zero.

Details

The plot-to-plot dissimilarity coefficients used in this function are as follows: "GC": Equation 6 in Pavoine and Ricotta (2019) "MS": Equation 8 in Pavoine and Ricotta (2019) "PE": Equations 9 and 10 in Pavoine and Ricotta (2019)

Value

The function returns an object of class "dist" with the values of the proposed dissimilarities for each pair of plots.

Author(s)

Sandrine Pavoine [e-mail]

References

Pavoine, S. and Ricotta, C. (2019) Measuring functional dissimilarity among plots: adapting old methods to new questions. Ecological Indicators, 97, 67–72.

Examples

## Not run: 
if(require(ade4) && require(adephylo) && require(ape)){

data(birdData)
phy <- read.tree(text=birdData$tre)
phydis <- sqrt(distTips(phy, method="nNodes")+1)
fau <- birdData$fau[1:6, phy$tip.label]
disGC <- generalized_Tradidiss(fau, phydis, method="GC")
disGC

### The second example is a bit TIME CONSUMING
data(mafragh)
namspe <- rownames(mafragh$traits[[1]])
M <- mafragh$flo
colnames(M) <- namspe
Bin <- prep.binary(mafragh$traits$tabBinary, c(3, 4))
distraits <- dist.ktab(ktab.list.df(list(mafragh$traits$tabOrdinal[,2:3], Bin)), 
    c("O","B"), scan=FALSE)
disGC <- generalized_Tradidiss(M, distraits, method="GC")
pcoGC <- dudi.pco(as.dist(cailliez(disGC)), full=TRUE)
s.value(mafragh$xy, pcoGC$li[,1])
}

## End(Not run)

adiv documentation built on May 29, 2024, 7:11 a.m.