plotPEarcs1D {pcds} | R Documentation |
The plot of the arcs of Proportional Edge Proximity Catch Digraphs (PE-PCDs)
for 1D data (vertices jittered along y
-coordinate) - multiple interval case
Description
Plots the arcs of PE-PCD whose vertices are the 1D points, Xp
.
PE proximity regions are constructed with
expansion parameter r \ge 1
and centrality parameter c \in (0,1)
and the intervals are based on Yp
points (i.e.
the intervalization is based on Yp
points).
That is, data set Xp
constitutes the vertices of the digraph and Yp
determines the end points of the intervals.
If there are duplicates of Yp
or Xp
points,
only one point is retained for each duplicate value,
and a warning message is printed.
For better visualization,
a uniform jitter from U(-Jit,Jit)
(default for Jit=.1
) is added to
the y
-direction where Jit
equals to the range of Xp
and Yp
multiplied by Jit
with default for Jit=.1
).
centers
is a logical argument, if TRUE
,
plot includes the centers of the intervals
as vertical lines in the plot,
else centers of the intervals are not plotted.
See also (Ceyhan (2012)).
Usage
plotPEarcs1D(
Xp,
Yp,
r,
c = 0.5,
Jit = 0.1,
main = NULL,
xlab = NULL,
ylab = NULL,
xlim = NULL,
ylim = NULL,
centers = FALSE,
...
)
Arguments
Xp |
A |
Yp |
A |
r |
A positive real number which serves as the expansion parameter in PE proximity region;
must be |
c |
A positive real number in |
Jit |
A positive real number that determines the amount of jitter along the |
main |
An overall title for the plot (default= |
xlab , ylab |
Titles of the |
xlim , ylim |
Two |
centers |
A logical argument, if |
... |
Additional |
Value
A plot of the arcs of PE-PCD whose vertices are the 1D data set Xp
in which vertices are jittered
along y
-axis for better visualization.
Author(s)
Elvan Ceyhan
References
Ceyhan E (2012). “The Distribution of the Relative Arc Density of a Family of Interval Catch Digraph Based on Uniform Data.” Metrika, 75(6), 761-793.
See Also
plotPEarcs.int
and plotCSarcs1D
Examples
r<-2
c<-.4
a<-0; b<-10; int<-c(a,b)
#nx is number of X points (target) and ny is number of Y points (nontarget)
nx<-20; ny<-4; #try also nx<-40; ny<-10 or nx<-1000; ny<-10;
set.seed(1)
xf<-(int[2]-int[1])*.1
Xp<-runif(nx,a-xf,b+xf)
Yp<-runif(ny,a,b)
Xlim=range(Xp,Yp)
Ylim=.1*c(-1,1)
jit<-.1
set.seed(1)
plotPEarcs1D(Xp,Yp,r=1.5,c=.3,jit,xlab="",ylab="",centers=TRUE)
set.seed(1)
plotPEarcs1D(Xp,Yp,r=2,c=.3,jit,xlab="",ylab="",centers=TRUE)