pcds (version 0.1.6)

PEarcdens.tetra: Arc density of Proportional Edge Proximity Catch Digraphs (PE-PCDs) - one tetrahedron case

Description

Returns the arc density of PE-PCD whose vertex set is the given 2D numerical data set, Xp, (some of its members are) in the tetrahedron th.

PE proximity region is constructed with respect to the tetrahedron th and vertex regions are based on the center M which is circumcenter ("CC") or center of mass ("CM") of th with default="CM". For the number of arcs, loops are not allowed so arcs are only possible for points inside the tetrahedron th for this function.

th.cor is a logical argument for tetrahedron correction (default is TRUE), if TRUE, only the points inside the tetrahedron are considered (i.e., digraph induced by these vertices are considered) in computing the arc density, otherwise all points are considered (for the number of vertices in the denominator of arc density).

See also (ceyhan:Phd-thesis,ceyhan:comp-geo-2010;textualpcds).

Usage

PEarcdens.tetra(Xp, th, r, M = "CM", th.cor = FALSE)

Value

Arc density of PE-PCD whose vertices are the 2D numerical data set, Xp; PE proximity regions are defined with respect to the tetrahedron th and M-vertex regions

Arguments

Xp

A set of 2D points which constitute the vertices of the PE-PCD.

th

A \(4 \times 3\) matrix with each row representing a vertex of the tetrahedron.

r

A positive real number which serves as the expansion parameter in PE proximity region; must be \(\ge 1\).

M

The center to be used in the construction of the vertex regions in the tetrahedron, th. Currently it only takes "CC" for circumcenter and "CM" for center of mass; default="CM".

th.cor

A logical argument for computing the arc density for only the points inside the tetrahedron, th. (default is th.cor=FALSE), i.e., if th.cor=TRUE only the induced digraph with the vertices inside th are considered in the computation of arc density.

Author

Elvan Ceyhan

References

See Also

PEarcdens.tri and NumArcsPEtetra

Examples

if (FALSE) {
A<-c(0,0,0); B<-c(1,0,0); C<-c(1/2,sqrt(3)/2,0); D<-c(1/2,sqrt(3)/6,sqrt(6)/3)
tetra<-rbind(A,B,C,D)
n<-10  #try also n<-20

set.seed(1)
Xp<-runif.tetra(n,tetra)$g

M<-"CM"  #try also M<-"CC"
r<-1.5

NumArcsPEtetra(Xp,tetra,r,M)
PEarcdens.tetra(Xp,tetra,r,M)
PEarcdens.tetra(Xp,tetra,r,M,th.cor = FALSE)
}