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

[Martin Schlather]

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