Re: [R] How to use rfm.test ? (Package MarkedPointProcess)

From: Florent Bonneu <bonneu_at_cict.fr>
Date: Thu 26 Jan 2006 - 01:10:25 EST

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
>>bonneu@cict.fr
>>
>>
>>
>>
>
>
>
>



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 Received on Thu Jan 26 01:32:33 2006

This archive was generated by hypermail 2.1.8 : Thu 26 Jan 2006 - 08:07:08 EST