An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-sig-ecology/attachments/20120216/25166d47/attachment.pl>
taxonomic diversity and distinctness - taxondive() in vegan package
6 messages · Jari Oksanen, Gian Maria Niccolò Benucci
On 16/02/2012, at 18:39 PM, Gian Maria Niccol? Benucci wrote:
Hi everybody, I am here for asking some kind of help about taxonimic diversity and distinctness indices, from taxondive()<http://127.0.0.1:25979/library/vegan/html/taxondive.html>function in vegan. The function produced also expected values for Delta, Dela* and Delta + indices but no metodology about how these expected data are calculed is reported. How are they calculated?
Dear Giain Maria Niccol? Benucci,
The ?taxondive help page reads:
Clarke & Warwick (1998, 2001) suggested several alternative
indices of taxonomic diversity or distinctness. Two basic indices
are called taxonomic diversity (Delta) and distinctness (Delta^*):
Delta = (sum sum_{i<j} omega_{ij} x_i x_j)/(n (n-1) / 2)
Delta^* = (sum sum_{i<j} omega_{ij} x_i x_j)/(sum sum_{i<j} x_i x_j)
The equations give the index value for a single site, and
summation goes over species i and j. Here omega are taxonomic
distances among taxa, and x are species abundances, and n is the
total abundance for a site. With presence/absence data both
indices reduce to the same index Delta^+
The intention of this passage was to give equations for Delta and Delta^*, and explanation of Delta^+ (and the pdf version has nicer formatting of equations). If this is not enough to explain the methodology of calculation, please consult the original publications referred to in the help page and the R code of taxondive().
According Clarke and Warwick (1998) it is also possible to plot a confidence funnel of the result with plot.taxondive, where Delta * values are plotted against corresponding number of species at each site.
It is a long time since I read these papers, but I think C&W suggested plotting Delta^+ values (not Delta^*), and that plot is produced by plot() command for taxondive() results in vegan, like displayed in the examples of the help page of taxondive. However, the labels cannot be displayed with the standard plot() command: you have to do that yourself. Look at the model in vegan:::plot.taxondive to see how simple it is. Cheers, Jari Oksanen
I argued that I only 15 points, instead of 24 that is the overall sample number, becasue there is an overplotting...do is it possible? Moreover, there is the possibility to plot in the funnel diagram the name of the samples? Thank you to everyone, G. -- Gian Maria Niccol? Benucci, Ph.D. Department of Applied Biology University of Perugia [[alternative HTML version deleted]]
_______________________________________________ R-sig-ecology mailing list R-sig-ecology at r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-ecology
Jari Oksanen, Dept Biology, Univ Oulu, 90014 Finland jari.oksanen at oulu.fi, Ph. +358 400 408593, http://cc.oulu.fi/~jarioksa
An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-sig-ecology/attachments/20120217/87e98a40/attachment.pl>
4 days later
An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-sig-ecology/attachments/20120221/8524ba06/attachment.pl>
1 day later
On 21/02/2012, at 13:38 PM, Gian Maria Niccol? Benucci wrote:
Dear All, I just want to ask if someone can explane me how the steps are calculated in taxa2dist(). I would know how steps are calculated because it is foundamental to define the weights [I'm referring to Clarke & Warwick (1999)] given to taxonomic relationships. Why I do have the same step value for genus and class ranks? Below is reported the code.
When I implemented step lengths, I followed Clarke & Warwick (1999): probably they explain the scaling, and I would read that reference if I needed to explain how it works. The docs say that the step length is relative to proportional decrease in the number of classes. Does this apply to your genus and class ranks? Cheers, Jari Oksanen
Jari Oksanen, Dept Biology, Univ Oulu, 90014 Finland jari.oksanen at oulu.fi, Ph. +358 400 408593, http://cc.oulu.fi/~jarioksa
An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-sig-ecology/attachments/20120223/64808f62/attachment.pl>