I don't find where it's possible to configure "closer bounds" in the algorithm. I have obtained some results with
rfm.test(coord=steigerwald$coord,steigerwald$diam,MCrepetitions=19,n.hypo=10) but I don't know how to interpret them. Where is the result of the MC test ? Could you give me the syntax of rfm.test for this example (Steigerwald) that you have used in your paper "Detecting dependence between Marks and Locations of Marked Point Processes ? regards. Martin Schlather wrote: >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] >> >> >> >> > > > > ______________________________________________ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
