Re: [R-sig-eco] species richness

2012-05-31 Thread Ivailo
On Wed, May 30, 2012 at 10:42 PM, Linda Bürgi patili_bue...@hotmail.com wrote:

 Hi,


 I am interested in testing whether increasing the number of plant
 species sampled increases the number of herbivore species I find, irrespective
 of the number of herbivores collected.


 To do that, I was thinking of fixing the number of
 herbivores collected (e.g. 100) and randomly pulling 99 samples of 100
 herbivores from all possible combinations of 2 plant species, 3 plant species,
 4 plant species, etc. This should then yield a curve with number of plant
 species on the x axis and average number of herbivore species found on the y
 axis, always for a sample of 100 herbivores.


 Does such a function already exist?


 My data (see below) is in matrix form with columns representing
 herbivore species (14) and rows representing plant species (14), with the 
 numbers in the
 cells representing number of specimens collected per herbivore and plant
 species combination.


 I’m not quite sure how to tackle this….


 Thanks!



 Linda

Hi Linda,

your question whether increasing the number of plant species sampled
increases the number of herbivore species seems to require a
contingency table and a corresponding test. I am not sure, however,
how to treat the requirement irrespective of the number of herbivores
collected as usually the number of species increases with the number
of individuals sampled. Perhaps to test additionally if the proportion
of herbivores sampled is independent of the proportion of
hon-herbivores in the sample(s)?

Hope this helps,
Ivailo
-- 
UBUNTU: a person is a person through other persons.

___
R-sig-ecology mailing list
R-sig-ecology@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-ecology


[R-sig-eco] Chi square distribution of Mahalanobis distance -how many degrees of freedom

2012-05-31 Thread Jon Olav Skoien

Dear list,

The squared Mahalanobis distance is defined as:
D2 = (X-m)' C (X-m)
When this is used in ecological modelling, X refers to the environmental 
variables at a location i, m and C are the mean and inverted covariance 
matrix of the same variables at a set of locations where a species has 
been observed. It is then usually assumed that this squared distance is 
Chi square distributed with n-1 degrees of freedom, n being the number 
of variables, as described in the original paper of Clark et al. (1993). 
The same assumption is used in the mahasuhab function of adehabitat. 
However, when reading non-ecological texts about linear models, the 
Mahalanobis distance seems to be a special case of quadratic forms, 
which for an invertible covariance matrix are described to have rank(C) 
= n degrees of freedom, see e.g. Christensen (2011, p9) and different 
hits with the search string:

mahalanobis chi square degrees of freedom
I cannot verify the reliability of all these sites, but it is striking 
that in most ecological links the authors seem to use n-1 degrees of 
freedom whereas in most non-ecological links the authors seem to use n 
degrees of freedom.


Clark et al. (1993) do not give a reference in the description of the 
Chi square distribution, but earlier in the theory section there are 
some references to a book of Morrison (1976). Could this be the source? 
Does anyone know why it appears to be a discrepancy in the  degrees of 
freedom? I have not found any discussions about this, is there anything 
I have missed?


Thanks in advance,
Jon


Christensen, R., 2011. Plane answers to complex questions: The theory of 
linear models, 4 ed. New York: Springer Science.


Clark, J.D., Dunn, J.E.  Smith, K.G., 1993. A multivariate model of 
female black bear habitat use for a geographical information system. 
Journal of Wildlife Management 57, 519-526.


Morrison, D. F., 1976. Multivariate statistical methods. McGraw-Hill 
Book Co., New York, N.Y., 415 pp.


--
Jon Olav Skøien
Joint Research Centre - European Commission
Institute for Environment and Sustainability (IES)
Land Resource Management Unit

Via Fermi 2749, TP 440,  I-21027 Ispra (VA), ITALY

jon.sko...@jrc.ec.europa.eu
Tel:  +39 0332 789206

Disclaimer: Views expressed in this email are those of the individual and do 
not necessarily represent official views of the European Commission.

___
R-sig-ecology mailing list
R-sig-ecology@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-ecology