Skip to content
Prev 8347 / 29559 Next

existance of a specific non-overlapping marked spatial model in spatstat or other R package?

On 28/05/2010, at 8:53 AM, Quets Jan wrote:

            
The only way I can think of doing this is sequentially.  I.e. select a radius,
then repetitively select a point to add, rejecting the selected point if its
circle (with the chosen radius) intersects any of the existing circles.  If
after ``giveup'' attempts the function has failed to find an acceptable point,
it gives up (and throws an error).

I cobbled together some code to effect this; it took one minute of user time
to do 100 replications of choosing 100 points with non-overlapping circles
in the unit square, with the radii selected uniform-randomly from (1:5)/100
and giveup=1000.

I believe that it would require a lot of work and great cleverness to get
anything substantially faster.  Of course, I could be wrong --- I was once,
back in 1968 when I thought I'd made a mistake and I hadn't. :-)

	cheers,

		Rolf Turner

P. S.  Let me know if you want my code.  I suspect it's not much different
from what you have built yourself.

		R.
######################################################################
Attention: 
This e-mail message is privileged and confidential. If you are not the 
intended recipient please delete the message and notify the sender. 
Any views or opinions presented are solely those of the author.

This e-mail has been scanned and cleared by MailMarshal 
www.marshalsoftware.com
######################################################################