Re: [R] spatstat rmh problem

2014-09-16 Thread Sebastian Schutte
You should use: plot(density(rmh(mod,new.coef=c(1,200 Sorry, my bad, typo in the example code. (2) However, even when the correct call is given you still wind up with identical densities!!! Hm. I think this may be a bug; I'll will check with the other authors of spatstat and repor

Re: [R] spatstat rmh problem

2014-09-16 Thread Rolf Turner
There was indeed a bug in rmh() w.r.t. the "new.coeff" argument. The bug has been fixed and will not be present in the next release of spatstat. cheers, Rolf Turner On 16/09/14 16:30, Sebastian Schutte wrote: Thanks so much for your comments. Sorry for not having sent a running example from

Re: [R] spatstat rmh problem

2014-09-15 Thread Rolf Turner
OK. Two things are going wrong. (1) There is an error in your code. You are passing the new.coef argument to density() and not to rmh(). The function density() has no such argument, but has a "..." argument, so "new.coef" simply gets ignored. You should use: plot(density(rmh(mod,new.coef=

Re: [R] spatstat rmh problem

2014-09-15 Thread Sebastian Schutte
Thanks so much for your comments. Sorry for not having sent a running example from the start. Here it is: library (spatstat) #Load example data data(demopat) #Generate a random point pattern within the polygon set.seed(12345) pdat <- rpoint(200,win=demopat$window) #Generate a distmap, which wi

Re: [R] spatstat rmh problem

2014-09-15 Thread Rolf Turner
Your example is not reproducible. We don't have "cshape" or "im.pop" (and are possibly lacking other bits and pieces; I didn't check the details since the example fails to run from the get-go). Please provide a *reproducible* example. Also I am puzzled by the line mod <- ppm (ppp, ~ pop

[R] spatstat rmh problem

2014-09-15 Thread Sebastian Schutte
Dear R and spatstat developers, Thanks so much for the time and effort that you invest into this awesome software. I have a problem simulating from a Point Process Model in spatstat. In summary, the option "new.coef" should allow me to use a fitted model and change its beta coefficients before

[R] spatstat rmh problem

2014-09-15 Thread Sebastian Schutte
Dear R and spatstat developers, Thanks so much for the time and effort that you invest into this awesome software. I have a problem simulating from a Point Process Model in spatstat. In summary, the option "new.coef" should allow me to use a fitted model and change its beta coefficients before