[R] Ripley's K function within envelope--how to change the maximum distance?

From: Steve Markofsky <s.markofsky_at_ucl.ac.uk>
Date: Tue, 04 Mar 2008 10:37:59 +0000


Hi--I'm trying to run a Ripley's K analysis on a point pattern within a window of 1 square km. The maximum distance I want to use is the diagonal of the window, around 1400m, run in 50m increments.

The function defaults to a 250m maximum--I thought this was simply the way I had plotted it, but when I exported the data there was nothing over 250m.

I've looked into seq and ginterval, and I think that seq may be my answer, but I'm not exactly sure how to use it properly within the envelope function. Here's two versions of the code I'm working with:

ERip<-envelope(EforRip,fun=Kest,correction=c("Ripley"),nsim=99,rank=1,transform=expression(sqrt(./pi)),global=TRUE)

ERip<-envelope(EforRip,fun=Kest,r=seq(0,1400,50),correction=c("Ripley"),nsim=99,rank=1,transform=expression(sqrt(./pi)),global=TRUE)

I think I'm missing something obvious, but couldn't find the solution.

Thanks very much--

Steve

Steve Markofsky
Institute of Archaeology
University College London

Mobile: 07779-789-720
Email: s.markofsky_at_ucl.ac.uk



R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Tue 04 Mar 2008 - 10:46:07 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Thu 06 Mar 2008 - 05:30:19 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive