Thanks Adrian. I think my confusion is that there are lots of points that
are more than r units away.
Sweden is basically a tall, thin rectangle. It's ~260km wide at its
narrowest point, but it's 1400km tall. And my point density is pretty high,
so for any point, there are lots of other points at more than 160km (just
not in all directions). In other words, the maximum distance from any
point to the boundary is NOT 160km. Thus my confusion.
Also, I've been overriding the `xlim` too -- the x-axis of the plot window
runs all the way to the value I want (300km), and the null poisson model
runs to 300km, it's just the estimates of k(r) that stop.
I don't think I can send pics to the r-geos list, but I'll email the image
to you.
(Ordered your book, btw)
On Wed, Jul 27, 2016 at 3:36 AM Adrian Baddeley <
adrian.baddeley at curtin.edu.au> wrote:
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