[R-sig-Geo] determining used default bandwidth in pcf function (spatstat)

Rolf Turner r.turner at auckland.ac.nz
Mon Jun 20 23:07:46 CEST 2011


On 21/06/11 04:52, Quets Jan wrote:
> Dear readers,
>
> I would like to know the default value of the bandwidth of the Epanechnikov kernel which is used in calculating a pcf.
>
> For example
>
> pp = rSSI(0.1,50)     #a point pattern
> plot(pcf(pp))              #the pcf of this point pattern
>
> There is a kernel bandwidth calculated here according to Stoyan's rule of thumb. But what is its value, is it possible to retrieve this in a simple way?
>
> Stoyan's rule of thumb = bw = c/sqrt(lambda).
>
> lambda is the point density and is known,
> c is a value between 0.1 and 0.2 according to the documentation, but I don't know how c is exactly calculated.
>
> I could search for the book of Stoyan but perhaps there is a simpler way?

The documentation for pcf.ppp() states clearly that the default value of "c"
(given as argument "stoyan") is 0.15.

Does this tell you what you want to know?

     cheers,

         Rolf Turner



More information about the R-sig-Geo mailing list