I have a point pattern which is clearly inhomogeneous. Furthermore, the inhomogeneity has two components: a large scale effect and a local scale effect. I have constructed a Markov point process model to capture these two first order intensity variations. However, it is clear that the empirical pair correlation function should not be used to determine the second order behavior of the point pattern. So I simply modelled the second order behavior as a Poisson process. Now I fitted the model, obtained simulations from fitted model and try to compare the descrepancies in second order behavior between the simulated data and the real data.
I don't want to use intensity reweighted second order summary statistics as it will very likely to give highly biased results due to potential misfit of the first order intensity from the model.
My question is, in this scenario, does it make sense to compare second order behavior between the simulated data and the real data using any kind of unweighted second order summary statistics? If so, which one would be better? NND distribution or PCF?