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 > [EMAIL PROTECTED] > > -- Prof. Dr. habil. Martin Schlather Helmut-Schmidt-Universität Universität der Bundeswehr Hamburg Fachbereich WOW Holstenhofweg 85 D -- 22043 Hamburg email: [EMAIL PROTECTED] http://www2.unibw-hamburg.de/schlath phone: +49 (0)40 6541 2873 Fax : +49 (0)40 6541 2023 ______________________________________________ [email protected] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
