How to use rfm.test ? (Package MarkedPointProcess)

[Florent Bonneu]

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
[hidden email]

______________________________________________
[hidden email] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

[Martin Schlather-2]

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
> [hidden email]
>
>

--
Prof. Dr. habil. Martin Schlather
Helmut-Schmidt-Universität
Universität der Bundeswehr Hamburg
Fachbereich WOW
Holstenhofweg 85
D -- 22043 Hamburg
email: [hidden email]     http://www2.unibw-hamburg.de/schlath
phone: [phone]     Fax : [phone]

______________________________________________
[hidden email] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

[Florent Bonneu]

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
>>[hidden email]
>>
>>
>>    
>>
>
>
>  
>

______________________________________________
[hidden email] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html