c<-.4
a<-0; b<-10; int<-c(a,b)
Mc<-centMc(int,c)
rv.mid.int(6,c,int)
rv.mid.int(3,c,int)
rv.mid.int(13,c,int)
rv.mid.int(4,c,int)
rv.mid.int(0,c,int)
rv.mid.int(-3,c,int)
n<-10 #try also n<-20
xr<-range(a,b,Mc)
xf<-(xr[2]-xr[1])*.5
dat<-runif(n,a,b)
Rv<-vector()
for (i in 1:n)
Rv<-c(Rv,rv.mid.int(dat[i],c,int)$rv)
Rv
jit<-.1
yjit<-runif(n,-jit,jit)
Xlim<-range(a,b,dat)
xd<-Xlim[2]-Xlim[1]
plot(cbind(Mc,0),main="vertex region indices for the points", xlab=" ", ylab=" ",
xlim=Xlim+xd*c(-.05,.05),ylim=3*range(yjit),pch=".",cex=3)
abline(h=0)
points(dat,yjit)
abline(v=c(a,b,Mc),lty=2,col=c(1,1,2))
text(dat,yjit,labels=factor(Rv))
text(cbind(c(a,b,Mc),.02),c("rv=1","rv=2","Mc"))
rv.mid.int(6,c,int)
Run the code above in your browser using DataLab