interval.indices.set | R Documentation |
Returns the indices of intervals for all the points in 1D data set,
Xp
, as a vector.
Intervals are based on Yp
and left end interval is labeled as 1,
the next interval as 2, and so on.
If there are duplicates of Yp
points,
only one point is retained for each duplicate value,
and a warning message is printed.
interval.indices.set(Xp, Yp)
Xp |
A set of 1D points for which the indices of intervals are to be determined. |
Yp |
A set of 1D points from which intervals are constructed. |
The vector
of indices of the intervals in which points in the 1D data set, Xp
, reside
Elvan Ceyhan
a<-0; b<-10; int<-c(a,b)
#nx is number of X points (target) and ny is number of Y points (nontarget)
nx<-15; 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) #try also Yp<-runif(ny,a+1,b-1)
ind<-interval.indices.set(Xp,Yp)
ind
jit<-.1
yjit<-runif(nx,-jit,jit)
Xlim<-range(a,b,Xp,Yp)
xd<-Xlim[2]-Xlim[1]
plot(cbind(a,0), xlab=" ", ylab=" ",xlim=Xlim+xd*c(-.05,.05),ylim=3*c(-jit,jit),pch=".")
abline(h=0)
points(Xp, yjit,pch=".",cex=3)
abline(v=Yp,lty=2)
text(Xp,yjit,labels=factor(ind))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.