# NOT RUN {
A<-c(1,1); B<-c(2,0); C<-c(1.5,2);
Tr<-rbind(A,B,C);
n<-10
set.seed(1)
dat<-runif.tri(n,Tr)$g
dat<-matrix(dat,ncol=2)
M<-as.numeric(runif.tri(1,Tr)$g) #try also M<-c(1.6,1.0)
t<-.5 #try also t<-2
dat<-matrix(dat,ncol=2)
Xlim<-range(Tr[,1],dat[,1])
Ylim<-range(Tr[,2],dat[,2])
xd<-Xlim[2]-Xlim[1]
yd<-Ylim[2]-Ylim[1]
plotCSregsTri(dat[4,],Tr,t,M)
plotCSregsTri(dat,Tr,t,M)
if (dimension(M)==3) {M<-bary2cart(M,Tr)}
#need to run this when M is given in barycentric coordinates
plotCSregsTri(dat,Tr,t,M,main="CS Proximity Regions with t=.5",
xlab="",ylab="",xlim=Xlim+xd*c(-.05,.05),ylim=Ylim+yd*c(-.05,.05))
L<-rbind(M,M,M); R<-Tr
segments(L[,1], L[,2], R[,1], R[,2], lty=2)
txt<-rbind(Tr,M)
xc<-txt[,1]+c(-.02,.03,.03,.03)
yc<-txt[,2]+c(.02,.02,.02,.07)
txt.str<-c("A","B","C","M")
text(xc,yc,txt.str)
# }
Run the code above in your browser using DataLab