# NOT RUN {
A<-c(1,1); B<-c(2,0); C<-c(1.5,2);
Tr<-rbind(A,B,C);
n<-10 #try also n<-20
set.seed(1)
dat<-runif.tri(n,Tr)$g
M<-as.numeric(runif.tri(1,Tr)$g) #try also M<-c(1.6,1.0)
NumArcsPEtri(dat,Tr,r=1.5,M)
PEarcdens.tri(dat,Tr,r=1.5,M)
PEarcdens.tri(dat,Tr,r=1.5,M,tri.cor = FALSE)
NumArcsPEtri(dat,Tr,r=1,M)
PEarcdens.tri(dat,Tr,r=1,M)
NumArcsPEtri(dat,Tr,r=1.5,M)
PEarcdens.tri(dat,Tr,r=1.5,M)
r<-2
PEarcdens.tri(dat,Tr,r,M)
dat.fr<-data.frame(a=dat)
PEarcdens.tri(dat.fr,Tr,r,M)
dat.fr<-data.frame(a=Tr)
PEarcdens.tri(dat,dat.fr,r,M)
# }
Run the code above in your browser using DataLab