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

Florent Bonneu
Wed, 25 Jan 2006 10:35:31 -0000

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