# Incidence matrix for Proportional Edge Proximity Catch Digraphs (PE-PCDs) - multiple triangle case Returns the incidence matrix of Proportional Edge Proximity Catch Digraph (PE-PCD) whose vertices are the data points in `Xp` in the multiple triangle case. PE proximity regions are defined with respect to the Delaunay triangles based on `Yp` points with expansion parameter $r \ge 1$ and vertex regions in each triangle are based on the center $M=(\alpha,\beta,\gamma)$ in barycentric coordinates in the interior of each Delaunay triangle or based on circumcenter of each Delaunay triangle (default for $M=(1,1,1)$ which is the center of mass of the triangle). Convex hull of `Yp` is partitioned by the Delaunay triangles based on `Yp` points (i.e., multiple triangles are the set of these Delaunay triangles whose union constitutes the convex hull of `Yp` points). For the incidence matrix loops are allowed, so the diagonal entries are all equal to 1. See (if(!exists(".Rdpack.currefs")) .Rdpack.currefs \<-new.env();Rdpack::insert_citeOnly(keys="ceyhan:Phd-thesis,ceyhan:arc-density-PE,ceyhan:dom-num-NPE-Spat2011;textual",package="pcds",cached_env=.Rdpack.currefs) ) for more on the PE-PCDs. Also, see (if(!exists(".Rdpack.currefs")) .Rdpack.currefs \<-new.env();Rdpack::insert_citeOnly(keys="okabe:2000,ceyhan:comp-geo-2010,sinclair:2016;textual",package="pcds",cached_env=.Rdpack.currefs) ) for more on Delaunay triangulation and the corresponding algorithm. ```r inci.matPE(Xp, Yp, r, M = c(1, 1, 1)) ``` ## Arguments - `Xp`: A set of 2D points which constitute the vertices of the PE-PCD. - `Yp`: A set of 2D points which constitute the vertices of the Delaunay triangles. - `r`: A positive real number which serves as the expansion parameter in PE proximity region; must be $\ge 1$. - `M`: A 3D point in barycentric coordinates which serves as a center in the interior of each Delaunay triangle or circumcenter of each Delaunay triangle (for this, argument should be set as `M="CC"`), default for $M=(1,1,1)$ which is the center of mass of each triangle. ## Returns Incidence matrix for the PE-PCD with vertices being 2D data set, `Xp`. PE proximity regions are constructed with respect to the Delaunay triangles and `M`-vertex regions. ## Examples ```r nx<-20; ny<-5; #try also nx<-40; ny<-10 or nx<-1000; ny<-10; set.seed(1) Xp<-cbind(runif(nx,0,1),runif(nx,0,1)) Yp<-cbind(runif(ny,0,.25), runif(ny,0,.25))+cbind(c(0,0,0.5,1,1),c(0,1,.5,0,1)) #try also Yp<-cbind(runif(ny,0,1),runif(ny,0,1)) M<-c(1,1,1) #try also M<-c(1,2,3) r<-1.5 #try also r<-2 IM<-inci.matPE(Xp,Yp,r,M) IM dom.num.greedy(IM) #try also dom.num.exact(IM) #might take a long time in this brute-force fashion ignoring the #disconnected nature of the digraph inherent by the geometric construction of it ``` ## References if(!exists(".Rdpack.currefs")) .Rdpack.currefs \<-new.env();Rdpack::insert_all_ref(.Rdpack.currefs) ## See Also `inci.matPEtri`, `inci.matPEstd.tri`, `inci.matAS`, and `inci.matCS` ## Author(s) Elvan Ceyhan

inci.matPE function