inci.matPEtri {pcds} | R Documentation |
Incidence matrix for Proportional Edge Proximity Catch Digraphs (PE-PCDs) - one triangle case
Description
Returns the incidence matrix for the PE-PCD
whose vertices are the given 2D numerical data set, Xp
,
in the triangle tri
.
PE proximity regions are constructed with respect to triangle tri
with expansion parameter
and vertex regions are based on the center
in Cartesian coordinates
or
in barycentric coordinates
in the interior of the triangle
tri
;
default is , i.e., the center of mass of
tri
.
Loops are allowed, so the diagonal entries are all equal to 1.
See also (Ceyhan (2005); Ceyhan et al. (2006); Ceyhan (2011)).
Usage
inci.matPEtri(Xp, tri, r, M = c(1, 1, 1))
Arguments
Xp |
A set of 2D points which constitute the vertices of PE-PCD. |
tri |
A |
r |
A positive real number
which serves as the expansion parameter in PE proximity region;
must be |
M |
A 2D point in Cartesian coordinates
or a 3D point in barycentric coordinates
which serves as a center in the interior of the triangle |
Value
Incidence matrix for the PE-PCD
with vertices being 2D data set, Xp
,
in the triangle tri
with vertex regions based on center M
Author(s)
Elvan Ceyhan
References
Ceyhan E (2005).
An Investigation of Proximity Catch Digraphs in Delaunay Tessellations, also available as technical monograph titled Proximity Catch Digraphs: Auxiliary Tools, Properties, and Applications.
Ph.D. thesis, The Johns Hopkins University, Baltimore, MD, 21218.
Ceyhan E (2011).
“Spatial Clustering Tests Based on Domination Number of a New Random Digraph Family.”
Communications in Statistics - Theory and Methods, 40(8), 1363-1395.
Ceyhan E, Priebe CE, Wierman JC (2006).
“Relative density of the random -factor proximity catch digraphs for testing spatial patterns of segregation and association.”
Computational Statistics & Data Analysis, 50(8), 1925-1964.
See Also
inci.matPE
, inci.matCStri
,
and inci.matAStri
Examples
A<-c(1,1); B<-c(2,0); C<-c(1.5,2);
Tr<-rbind(A,B,C);
n<-10
set.seed(1)
Xp<-runif.tri(n,Tr)$g
M<-as.numeric(runif.tri(1,Tr)$g) #try also M<-c(1.6,1.0)
IM<-inci.matPEtri(Xp,Tr,r=1.25,M)
IM
dom.num.greedy(IM) #try also dom.num.exact(IM)
Idom.num.up.bnd(IM,3)