[R-sig-Geo] `kest` radius limits?

[Adrian Baddeley]

Nick Eubank writes:

> I'm running kest on points in Sweden (which is about 200km wide at its
> narrowest) with 'border' correction and lots and lots of points.

> However, I can't get estimates of K(r) for radii beyond ~160km. I bumped
> `rmax`, and if I plot my results, the null model line extends  to the
> `rmax` value, but the K(r) estimates just stop at ~160km.

For the border correction, the estimate of K(r) becomes NA (and is not plotted) as soon as there are no data points that lie more than 'r' units away from the boundary. If the maximum distance from any data point to the boundary of Sweden is 160 km, then this is the behaviour I would expect.

All explained in chapter 7 of the spatstat book.

If you print a summary function object, the output shows the "recommended range" of distance values and the full "available range". By default, the plots of summary functions use the "recommended range" of distance values rather than the full "available range". To override this, when you plot the summary function, use 'xlim' to control the range of 'r' values. (xlim should be an argument in the 'plot' call).

Adrian Baddeley



	[[alternative HTML version deleted]]