[R-sig-Geo] problem with enveloped test in spatstat

[David Unwin]

Can any *spatstat* user explain to me why the two p-values obtained below
for an envelope test against CSR are so different?



> data(swedishpines)

> d<-swedishpines

> plot(d)

> mad.test(d,Kest,nsim=999,verbose=F)



        Maximum absolute deviation test of CSR

        Monte Carlo test based on 999 simulations

        Summary function: K(r)

        Reference function: theoretical

        Alternative: two.sided

        Interval of distance values: [0, 24] units (one unit = 0.1 metres)

        Test statistic: Maximum absolute deviation

        Deviation = observed minus theoretical



data:  d

mad = 150.69, rank = 216, p-value = *0.216*



> mad.test(d,*Lest*,nsim=999,verbose=F)



        Maximum absolute deviation test of CSR

        Monte Carlo test based on 999 simulations

        Summary function: L(r)

        Reference function: theoretical

        Alternative: two.sided

        Interval of distance values: [0, 24] units (one unit = 0.1 metres)

        Test statistic: Maximum absolute deviation

        Deviation = observed minus theoretical



data:  d

mad = 2.9921, rank = 9, p-value = *0.009*



*!!!*



These data are dispersed relative to CSR:

>kl<-envelope(d,Kest,nsim=999,correction="border")

>plot(kl)


Dave Unwin

	[[alternative HTML version deleted]]