t<-2
c<-.4
a<-0; b<-10;
nx<-10; ny<-4
set.seed(1)
Xp<-runif(nx,a,b)
Yp<-runif(ny,a,b)
IM<-IncMatCS1D(Xp,Yp,t,c)
IM
dom.greedy(IM)
# \donttest{
dom.exact(IM) #might take a long time depending on nx
# }
IndUBdom(IM,5)
Arcs<-ArcsCSMI(Xp,Yp,t,c)
Arcs
summary(Arcs)
plot(Arcs)
IncMatCS1D(Xp,Yp+10,t,c)
t<-2
c<-.4
a<-0; b<-10;
#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;
Xp<-runif(nx,a,b)
Yp<-runif(ny,a,b)
IncMatCS1D(Xp,Yp,t,c)
Run the code above in your browser using DataLab