Skip to content
Prev 665 / 7420 Next

spatstat: help with calculating number of clusters, cluster size

Hi Soumya,

There are several ways. If you are certain that 
your trees are clustered,  you could consider as 
a first option fitting Poisson cluster or 
Thomas  processes by minimum contrast methods. In 
spatstat you can do this, for example, with function thomas.estK .
 From the fitted parameters you will obtain both 
kappa (the intensity of "parent" points , i.e. an 
approximation to the "number" of clusters) and 
sigma (the standard deviation of the distribution 
of "offspring" points around the parents, so you 
could approximate the radius as 2*sigma).

You can also use  the spatstat-based package 
ecespa 
(http://cran.r-project.org/web/packages/ecespa/index.html); 
see the help pages of  the functions pc.estK  and 
ipc.estK  for complete examples of fitting 
homogeneous and inhomogeneous processes.

Cheers,

Marcelino
At 12:49 13/06/2009, Soumya Prasad wrote:
________________________________

Marcelino de la Cruz Rot

Departamento de  Biolog?a Vegetal
E.U.T.I. Agr?cola
Universidad Polit?cnica de Madrid
28040-Madrid
Tel.: 91 336 54 35
Fax: 91 336 56 56
marcelino.delacruz at upm.es
_________________________________