pcds (version 0.1.2)

IndCSTe: The indicator for the presence of an arc from a point to another for Central Similarity Proximity Catch Digraphs (CS-PCDs) - standard equilateral triangle case

Description

Returns \(I(\)pt2 is in \(N_{CS}(pt1,t))\) for points pt1 and pt2, that is, returns 1 if pt2 is in \(N_{CS}(pt1,t)\), returns 0 otherwise, where \(N_{CS}(x,t)\) is the CS proximity region for point \(x\) with expansion parameter \(t>0\).

CS proximity region is defined with respect to the standard equilateral triangle \(T_e=T(v=1,v=2,v=3)=T((0,0),(1,0),(1/2,\sqrt{3}/2))\) and edge regions are based on the center \(M=(m_1,m_2)\) in Cartesian coordinates or \(M=(\alpha,\beta,\gamma)\) in barycentric coordinates in the interior of \(T_e\); default is \(M=(1,1,1)\) i.e., the center of mass of \(T_e\). re is the index of the edge region pt1 resides, with default=NULL.

If pt1 and pt2 are distinct and either of them are outside \(T_e\), it returns 0, but if they are identical, then it returns 1 regardless of their locations (i.e., it allows loops).

See also (ceyhan:Phd-thesis,ceyhan:arc-density-CS,ceyhan:test2014;textualpcds).

Usage

IndCSTe(pt1, pt2, t, M = c(1, 1, 1), re = NULL)

Arguments

pt1

A 2D point whose CS proximity region is constructed.

pt2

A 2D point. The function determines whether pt2 is inside the CS proximity region of pt1 or not.

t

A positive real number which serves as the expansion parameter in CS proximity region.

M

A 2D point in Cartesian coordinates or a 3D point in barycentric coordinates which serves as a center in the interior of the standard equilateral triangle \(T_e\); default is \(M=(1,1,1)\) i.e. the center of mass of \(T_e\).

re

The index of the edge region in \(T_e\) containing the point, either 1,2,3 or NULL (default is NULL).

Value

\(I(\)pt2 is in \(N_{CS}(pt1,t))\) for pt1, that is, returns 1 if pt2 is in \(N_{CS}(pt1,t)\), returns 0 otherwise

References

See Also

IndNCStri and IndNPETe

Examples

# NOT RUN {
A<-c(0,0); B<-c(1,0); C<-c(1/2,sqrt(3)/2);
Te<-rbind(A,B,C);
n<-10

set.seed(1)
dat<-runifTe(n)$gen.points

M<-as.numeric(runifTe(1)$g)  #try also M<-c(.6,.2)

t<-1

IndCSTe(dat[1,],dat[2,],t,M)
IndCSTe(dat[1,],dat[1,],t,M)

IndCSTe(dat[1,],dat[2,],t=4,M)
IndCSTe(dat[2,],dat[5,],t,M)
IndCSTe(c(.2,.5),dat[2,],t,M)

IndCSTe(c(.2,.5),c(.2,.5),t,M)

#or try
re<-reTeCM(dat[1,])$re
IndCSTe(dat[1,],dat[2,],t,M,re=re)

IndCSTe(dat[1,],dat[2,],t,M)

# }