Skip to content
Prev 17389 / 29559 Next

spatstat kstest

Dear Rolf,

Unfortunately, the data I have is under nondisclosure agreement.
Therefore, I hesitate to offer.

Let X={(x1,y1),(x2,y2),..} be the data set, and X_x={x1,x2,...} be its
projection to the x-axis.
Try, for example, 
fit <- ppm(X, ~x+y+I(x^2)+I(y^2), Poisson())
(that is, the inhomogeneous Poisson process with the intensity
\lambda(x,y)=exp(a+bx+cy+dx^2+ey^2))
and test
kstest(fit, covariate="x")

In my understanding, this is the KS-test whether X_x and the
inhomogeneous Poisson process with the intensity \lambda(x)=\int
\lambda(x,y)dy.
Am I wrong ?

What I would like to say is that kstest is not a function of the test of
the agreement between the model defined over (x,y) and the data set on
the plane, but the function of the test of the agreement between the
projected model and the projected model on a one-dimensional space.
Am I wrong?

Even so, you may have a good way to use kstest to test the two
dimensional data by, for example, appropriately setting the covariate in
kstest.
If so, please give me some pieces of advice.

Regards,
Hiroshi Saito

On Fri, 08 Feb 2013 13:51:55 +1300
Rolf Turner <r.turner at auckland.ac.nz> wrote:

            
*************************
Hiroshi Saito:
E-MAIL: saito.hiroshi at lab.ntt.co.jp
http://www9.plala.or.jp/hslab/
PHONE: +81 422 59 4300
FAX: +81 422 59 6364