c<-.4
a<-0; b<-10; int<-c(a,b)
Mc<-centMc(int,c)
nx<-10
xr<-range(a,b,Mc)
xf<-(xr[2]-xr[1])*.5
dat<-runif(nx,a,b)
Ext<-cl2Mc.int(dat,int,c)
Ext
summary(Ext)
plot(Ext)
cl2Mc.int(dat[1],int,c)
cMc<-cl2Mc.int(dat,int,c)
Xlim<-range(a,b,dat)
xd<-Xlim[2]-Xlim[1]
plot(cbind(a,0),xlab="",pch=".",xlim=Xlim+xd*c(-.05,.05))
abline(h=0)
abline(v=c(a,b,Mc),col=c(1,1,2),lty=2)
points(cbind(dat,0))
points(cbind(c(cMc$ext),0),pch=4,col=2)
text(cbind(c(a,b,Mc),-0.1),c("a","b","Mc"))
n<-10 #try also n<-20
dat<-runif(n,a-5,b+5)
cl2Mc.int(dat,c(a,b),c)
dat<-runif(n,a+b,b+10)
cl2Mc.int(dat,int,c)
c<-.4
a<-0; b<-10; int<-c(a,b)
n<-10
dat<-runif(n,a,b)
cl2Mc.int(dat,int,c)
Run the code above in your browser using DataLab