[R] How to use rfm.test ? (Package MarkedPointProcess)
Florent Bonneu
bonneu at cict.fr
Wed Jan 25 11:35:31 CET 2006
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
More information about the R-help
mailing list