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]]
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
