The method by Fry for detecting geometric anisotropy in stationary spatial point pattern is investigated. We quantify anisotropy by stretching and compressing the point process about the axis. Using a simulated Strauss point pattern, we first fit ellipsoids to the compressed pattern of pairwise difference vectors to estimate the direction of anisotropy. The strength of compression and the regularity of the point process are varied at different times and the corresponding effect on the estimated ellipsoids is investigated as the level of contour ellipsoid increases. Parameters and contours good for estimation were discussed.
Stationarity and isotropy (translation and rotational invariance of the distribution) have been the assumptions of spatial point pattern that catch the attention of many readers. However, anisotropy has not gotten much attention. The Fry plot gives more idea of anisotropy by creating an empty space of ellipse at the center if the point process is regular, for clustered process, the empty space turns out to be black depending on the shape of the clusters [
The main task of this paper is to evaluate geometric anisotropy estimators for spatial point patterns. The method used is the Fry plot method, thus the plot of pairwise difference vectors of all points of the process. We will create a series of sectors (from point process) beaming out of the origin in a specified direction and count the number of points falling into each of these sectors [
X is assumed to be a stationary and isotropic point process; for a Borel set; B ⊆ R d , let μ be the Lebesgue measure on the Borel set B. See ( [
N X ( B ) : = ∑ x ∈ X 1 B ( x ) (1)
be the counting measure defined by X and v ( B ) = μ ( B ) be the volume of B.
Additionally, we assume that the point pattern observed is of the form;
y : = { y i ; y i ∈ W , i = 1 , ⋯ , n } (2)
W ⊂ R d (where usually d = 2 or d = 3 in applications) is the observation window and comes from the transformed process;
Y : = T X = { T x ; x ∈ X } (3)
where T : R d → R d is an invertible linear mapping. Due to isotropy of X, we can split the mapping into two matrices (see [
T = R C (4)
where R is the rotation matrix that rotates the point process about the origin and C is the diagonal scaling matrix for compression and stretching of dimensions.
Let;
F ( Y ) : = { f i j = y i − y j ; y i , y j ∈ Y } (5)
be the Fry points in R d (see [
l l : = { ( r , u ) ∈ R + × S d − 1 ; r = arg min t { ρ ( 2 ) ( t u ) = l } } (6)
are spheres where l ′ s are the range of ρ ( 2 ) and u is the direction vector [
N Y ( u , r ) = ∑ f i j ∈ F ( Y ) 1 ( f i j ∈ S ε ( u , r ) ) (7)
which counts the Fry points within this r-sector. Let the distance
r l ( u ) : = arg min t { ( N Y ( u , t ) ) = l } (8)
for l = 1 , 2 , ⋯ defines the minimum distance in the sector S ε ( u ) having l points.
Thomas A. et al. (2016) [
the spheroidal transformation is preserved, C = d i a g ( c − 1 2 , c − 1 2 , c ) . i.e. only one free parameter; 0 < c < 1 which compresses the pattern top to down.
In our simulation study, we considered a stationary Strauss process with fixed number of points at n = 300 and the interaction radius were fixed at r = 0.04. We then varied the regularity parameter between γ = 0.0 (most regular case), γ = 0.1 and γ = 0.3 (least regular case). Concerning the parameters for the Fry plot for estimating anisotropy, we varied the strength of compression from c =
0.5, 0.6, 0.7, 0.8, 0.9, 1.0 with the compression model ( 1 c , c ) for compression
in the y-direction and dilation in the x-direction. The number of sectors were fixed at angles = 50 and for the observation window, we considered a unit cube [ 0 , 1 ] × [ 0 , 1 ] . For the level of contours, we considered from 1 up to 20. However, we extented the contours up to 30 in our further simulations to confirm our observation for higher level of contours. We simulated 100 times.
The rotangle is the angle which checks the direction of compression. Since the compression was done vertically, then we expect to get rotation angle around π 2 . If the compression is diagonal (anticlockwise), then you should get rotangle = π 4 .
We started by fixing the regularity parameter, γ and the strength of compression was being varied for only one simulation. We plotted the rotangle against the contours. The marked rotangle was indicated with black horizontal line at
π 2 and different colours were chosen for different strength of compressions.
From our plot, when gamma was fixed at 0.0, only c = 0.5 and c = 0.6 were the parameters for which the rotangle were close to the marked rotangle at π 2 .
From c = 0.7 up to 1.0 the rotangle deviates from the marked rotangle and hence the contours are not informative. The higher levels of contour are more informative than the lower contours.
Now we fixed the strength of compression and the regularity parameter was being varied for only one simulation.
We observed some breakage when compression was fixed at c = 0.5 ( γ = 0.0 ) because some values were missing. From the plot, we observed that the contours were informative at the instances where the compression was higher, thus when compression was fixed at c = 0.5, 0.6, 0.7. At low strength of compression, the contours look more spherical like expected. The worst case was when there was no compression at c = 1.0. As the strength of compression increases, then the contours become more elliptical.
Now we fixed gamma and compression was being varied. We chose the contour levels l = 5, 10, 15, 20, 25, 30 due to higher number of contours and the simulation was done 100 times. The box plot in
From our plot, we observed that as the strength of compression increases from c = 1.0 up to c = 0.5 for each contour level l, the average of the rotangle is around
the marked rotangle at π 2 . Also down the plot as the level of contour increases
from l = 5 up to l = 30, our box plots become shorter and hence the estimated contours look more elliptical. This means that at high level of contours and at high strength of compression, the contours are informative. The lower level contours look spherical in all parameters.
Now we proceed by fixing the strength of compression c and vary the regularity parameter γ . We simulated 100 times and again we chose the contour level l = 5, 10, 15, 20, 25, 30. The box plot in
From our plot, we observed that for each contour level l, as the regularity parameter (γ) of the point process decreases from 0.0 to 0.3, our box plots become
taller and thus the overall rotangle deviates from the marked rotangle at π 2 .
Again, at low strength of compression say at c = 0.9 and irrespective of the contour level l, the box plot becomes taller meaning the overall rotangle deviates from the marked rotangle. This means that in most regular point process and at high strength of compression, the contours are more informative.
The p-value test of equality of the semi-axes is simply a statistical test that tests if the semi axes of the ellipse in the Fry plot are equal. The contrast of interest (null hypothesis) is that the semi axes are equal ( a 1 − a 2 = 0 ) where a i is the length of ith semi-axes (in 2D). Let’s assume that the null hypothesis is true and the value of the test statistic has being calculated. If p-value is higher than α then it means we can accept equality which means that the ellipse is in reality a ball meaning there was no compression. If there is compression and equality is accepted, then the chosen contour is not informative. If the p-value is less than (or equal to) α , we reject the null hypothesis in favour of the alternative hypothesis. ( a 1 − a 2 ≠ 0 ). We tested at 5% level of significance in all simulations ( α = 0.05 ).
We fixed the regularity parameter of the point process and the strength of compression was being varied for one simulation. We plotted the p-value against the levels of contour. The significance level α = 0.05 was indicated with a black horizontal line and different colours were chosen for different strength of compression.
From our plot and for each fixed value of gamma, we observed that at high strength of compression (say c = 0.5 and c = 0.6) and as the level of contour increases, the smaller the p-value (below α = 0.05 ). This means that at high strength of compression and at high level of contour, we reject the null hypothesis that the semi axes of the estimated ellipsoids are equal and hence the contours are informative in these parameters. Also, the p-value increases as we vary the regularity parameter from 0.0 to 0.3 meaning the contours in the least regular patterns are not informative. Low level contours are not as informative as the high contour levels in all parameters.
Now we fixed compression and vary the regularity parameter of the point process and with the same level of significance at α = 0.05 and the same number of contours (nvec = 1:20). The plot in
From our plot and for each fix value of compression, we observed that at high level of contour and in more regular case ( γ = 0.0 and γ = 0.1 ), the p-value of the test decreases (below α = 0.05 ). This means that we can reject the null hypothesis that the semi axes of the estimated ellipsoids are equal in these parameters and hence the contours are informative. Again, at low strength of compression (c = 0.8, 0.9, 1.0), the p-values increase and therefore accepting equality of the semi axes of the estimated ellipsoids in these parameters. For c = 1 (no compression), the test confirms the correct shape of balls/spheres as expected (isotropic case).
Now we fixed gamma and compression was being varied. We chose the contour levels l = 5, 10, 15, 20, 25, 30 due to high number of contours and the simulation was done 100 times. We tested at 5% ( α = 0.05 ). The box plot in
From our plot, we observed that for each fix level of contour, as the strength of compression increases from c = 1.0 to c = 0.5, the overall p-values of the test decrease. Further, as we go down the plot and as the level of contour increases from l = 5 to l = 30, the overall p-value also decreases. This means that at high strength of compression and at high level of contours, we reject the null hypothesis that the semi axes of the estimated ellipsoids are equal and hence the contours are informative in these parameters. Further, as we move across the plot for each fixed level of contour and as we decrease the regularity of the point process, the overall p-value increases and hence the contours are not informative in these parameters.
Now we fixed compression and gamma was being varied. We chose the contour levels l = 5, 10, 15, 20, 25, 30 due to high number of contours and the simulation was done 100 times. We tested at 5% ( α = 0.05 ). The box plot in
From our plot, we observed that for each fix level of contour, the p-value of our test increases as the regularity of the point process decreases. Also as we go down the plot and as the level of contour increases for fix strength of compression, the overall p-value also decreases. This means that in more regular point pattern and at high level of contour, we reject equality of the semi axes of the estimated ellipsoids and hence the contours are informative in these parameters. As we move across the plot and as the strength of compression decreases for fix level of contour, the overall p-value increases and hence the contours are not informative in these parameters.
In theory, the estimated ellipsoids should be more spherical as the level of contour increases [
Again, in all our results, we found out that at low strength of compression (c = 0.8, 0.9, 1.0), the estimated ellipsoids get more spherical as expected. At high strength of compression (c = 0.5, 0.6, 0.7), we found out that our estimated ellipsoids are informative rendering our contours meaningful in these set of parameters. c = 0.5 is the best case.
Furthermore, in all our results, we found out that as the regularity of the point process increases (γ = 0.0, 0.1), the estimated ellipsoids get more elliptical and hence the chosen level of contour is informative in these set of parameters. γ = 0.0 is the best case. In the least regular case (γ = 0.3), the estimated ellipsoids look more spherical and hence the contours are not informative in this parameter. In summary, we found out that high levels of contour should be considered in all parameters.
In conclusion and due to the model we considered, at high level of contour coupled with high strength of compression (c = 0.5, 0.6, 0.7) and in more regular point process (γ = 0.0, 0.1), our estimated ellipsoids are more elliptical (not like sphere) and hence these parameters are good for estimating anisotropy geometrically by the method of Fry. At low level of contour coupled with low strength of compression (c = 0.8, 0.9, 1.0) and in low regular point process (γ = 0.3), the estimated ellipsoids look more spherical like expected.
The author declares no conflicts of interest.
Kumi, W. (2021) Simulation Study on Geometric Anisotropic Estimators for Spatial Point Process. Open Access Library Journal, 8: e8095. https://doi.org/10.4236/oalib.1108095