Florent Bonneu
2006-Jan-25 10:35 UTC
[R] How to use rfm.test ? (Package MarkedPointProcess)
I would like to compute the MC test (rfm.test) available in the package MarkedPointProcess (for the data BITOEK for example) in order to test the dependence between the marks and their locations. Why the syntax of rfm.test is false here? I have the message : ****************************** ML WARNING! Forbidden values! -- if there are too many warnings try narrower lower and upper bounds for the variables. -9.802347e-17 1364.372 [ 0 0.01808660 , 9.990131 1364.372 ] ****************************** The program : ----- # Packages # library(spatstat) library(RandomFields) library(adapt) library(MarkedPointProcess) # Program # data(BITOEK) win <- ripras(steigerwald$coord) PointProcess <- ppp(x=steigerwald$coord[,1],y=steigerwald$coord[,2],window=win,marks=steigerwald$diam) plot(PointProcess) rfm.test(coord=steigerwald$coord,steigerwald$diam,MCrepetitions=19) ----- regards. Florent Bonneu Laboratoire de Statistique et Probabilit??s bureau 148 b??t. 1R2 Universit?? Toulouse 3 118 route de Narbonne - 31062 Toulouse cedex 9 bonneu at cict.fr
Martin Schlather
2006-Jan-25 10:59 UTC
[R] How to use rfm.test ? (Package MarkedPointProcess)
Dear Florent Bonneu, the optim algorithm with parameter method="L-BFGS-B" (used in rfm.test) does not stay always exactly within the given bounds during the search of the optimum. This happens more frequently when the bounds are too wide. Here. rfm.test notices that -9.802347e-17 is less then the given lower bound of value 0. If the messages appear frequently it appeared to me that the algorithms runs into boundary, local minima more frequently. So, if the messages appear only a very few times, ignore them. If they appear several times, use closer bounds. Cheers, Martin Florent Bonneu wrote:> I would like to compute the MC test (rfm.test) available in the > package MarkedPointProcess (for the data BITOEK for example) in order > to test the > dependence between the marks and their locations. Why the syntax of > rfm.test is false here? I have the message : > > > ****************************** > ML WARNING! Forbidden values! -- if there are too many warnings try > narrower lower and upper bounds for the variables. -9.802347e-17 > 1364.372 [ 0 0.01808660 , 9.990131 1364.372 ] > ****************************** > > > The program : > ----- > # Packages # > library(spatstat) > > library(RandomFields) > library(adapt) > library(MarkedPointProcess) > > > # Program # > data(BITOEK) > win <- ripras(steigerwald$coord) > PointProcess <- > ppp(x=steigerwald$coord[,1],y=steigerwald$coord[,2],window=win,marks=steigerwald$diam) > > plot(PointProcess) > > > rfm.test(coord=steigerwald$coord,steigerwald$diam,MCrepetitions=19) > ----- > > > regards. > > Florent Bonneu > Laboratoire de Statistique et Probabilit??s > bureau 148 b??t. 1R2 > Universit?? Toulouse 3 > 118 route de Narbonne - 31062 Toulouse cedex 9 > bonneu at cict.fr > >-- Prof. Dr. habil. Martin Schlather Helmut-Schmidt-Universit??t Universit??t der Bundeswehr Hamburg Fachbereich WOW Holstenhofweg 85 D -- 22043 Hamburg email: schlather at hsu-hh.de http://www2.unibw-hamburg.de/schlath phone: +49 (0)40 6541 2873 Fax : +49 (0)40 6541 2023