Hello all: Is there any way to add points to a point pattern, while keeping a given minimum distance amongst this new points and the pre-existing points? Thanks, Frederico [[alternative HTML version deleted]]
Perhaps.
What does this have to do with R?
---------------------------------------------------------------------------
Jeff Newmiller The ..... ..... Go Live...
DCN:<jdnewmil at dcn.davis.ca.us> Basics: ##.#. ##.#. Live
Go...
Live: OO#.. Dead: OO#.. Playing
Research Engineer (Solar/Batteries O.O#. #.O#. with
/Software/Embedded Controllers) .OO#. .OO#. rocks...1k
---------------------------------------------------------------------------
Sent from my phone. Please excuse my brevity.
Frederico Mestre <mestre.frederico at gmail.com> wrote:
>Hello all:
>
>
>
>Is there any way to add points to a point pattern, while keeping a
>given
>minimum distance amongst this new points and the pre-existing points?
>
>
>
>Thanks,
>
>
>
>Frederico
>
>
>
>
> [[alternative HTML version deleted]]
>
>______________________________________________
>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.
On 04/09/12 11:58, Frederico Mestre wrote:> Hello all: > > > > Is there any way to add points to a point pattern, while keeping a given > minimum distance amongst this new points and the pre-existing points? >Presumably this is a question about the "spatstat" package. If so, the answer is "yes". At least two ways: (1) Using rmh() --- need to specify "beta" (the "chemical" activity parameter). E.g. # Generate a "test pattern" to which to add points. set.seed(42) X <- rSSI(0.05,100) # Add points. M <- rmhmodel(cif="hardcore",par=list(beta=300,hc=0.05)) Y <- rmh(M,start=list(x.start=X),expand=1,control=list(x.cond=as.data.frame(X))) plot(Y,main="Pattern with points added via rmh") plot(X,add=TRUE,chars=20,cols="red") print(min(nndist(Y))) (2) Using rSSI() --- need to specify the *total number* of points desired. E.g. # Add points to the same "test pattern", in a different way. Z <- rSSI(0.05,npoints(X)+42,x.init=X) plot(Z,main="Pattern with points added via rSSI") plot(X,add=TRUE,chars=20,cols="red") print(min(nndist(Z))) HTH cheers, Rolf Turner