Hello,
I was uclear before, I'm sory about it. I forgot to add that I'm using duleg...
I used mvpart for multivariate regression trees. My input variables are
environmental parameters, output variables are macrophyte species
(presence=1,absence=0 in conecutive cases=lakes). For obtained classes I used
duleg to find indicator species for every class. I checked the article Dufrene,
M. and Legendre, P. 1997. Species assemblages and indicator species: the need
for a flexible asymmetrical approach. Ecol. Monogr. 67(3):345-366. The authors
used the threshold of indval=0.25(25%) and that's the only hint I've found in
the literature. This threshod seems to reasonable, but still I have impression
that's too low...
best regards
Agnieszka
As Jari indicated, it depends on which function you meant in you
inquiry. The duleg() function implements the Dufrene-Legendre
algorithm, where "indicator" species are indicative of a priori
communities. It this requires a classification, and is biased to find
species which occur in the dataset approximately as often as the mean
cluster size.
The indpsc() function calculates the mean similarity of all samples
a species occurs in. This is slightly biased because we know that the
samples being used to calculate the mean share at least the species that
defines them, but it is still possible to compare those values to the
mean similarity of the whole matrix, or to an expectation of maximum
similarity. Obviously, as species occur more frequently, the harder it
is to have a really high similarity (indicator value), with the extreme
case that a species that occurs in every sample must have the same value
as the mean of the whole matrix.
To tell the truth, I forgot that indspc() was included in the
current version of labdsv. In the new version (due to be released any
day), I have included a permutation test that estimates quantiles of
expected values for different numbers of occurrences. It works, but is
pretty slow. Jari has created a version that uses parametric statistics
to estimate the same envelope, but I haven't had a chance to try it yet.
What research are you doing, and what are you really trying to
determine? Perhaps something altogether different will work better.
On Mon, 2005-09-19 at 09:41 +0200, astrzelczak at ps.pl wrote:
Hi,
I'm trying to find out what threshold of indicator value in labadsv should be
used to accept a specie as an indicator one? So far I assumed that indval=0.5
is high enough to avoid any mistakes but it was based only in my intuition.
I'd be greatful for any advise
best regards
Agnieszka,
R mailing list software appends the following to your message: