N<-c(32,32) eg<-sim.data(type="mulmod",N=N) lf<-leafsfirst(eg) ngrid<-4 lf.redu<-treedisc(lf,eg,ngrid=ngrid) lf.plot<-lf.redu lf.plot$level<-c(0,1,1,1,2,2,3) # frame 1 plottree(lf.plot,modelabel=FALSE,axes=FALSE) # frame 2 stepsi<-lf$maxdis/(ngrid+1) rad<-seq(stepsi,lf$maxdis-stepsi,stepsi) roundrad<-round(rad,digits=3) dm<-draw.pcf(eg) contour(dm$x,dm$y,dm$z,levels=roundrad,drawlabels=FALSE,axes=FALSE)
N<-c(32,32) eg<-sim.data(type="mulmod",N=N) lf<-leafsfirst(eg) ngrid<-4 lf.redu<-treedisc(lf,eg,ngrid=ngrid) lf.plot<-lf.redu lf.plot$level<-c(0,1,1,1,2,2,3) info<-matrix("",length(lf.plot$level),1) for (i in 1:length(lf.plot$level)){ info[i]<-paste(as.character(round(lf.plot$center[1,i],digits=1)), as.character(round(lf.plot$center[2,i],digits=1))) } plottree(lf.plot,modelabel=TRUE,ptext=0.3,axes=FALSE,info=info,infopos=-0.2, xmarginright=0.1)
N<-c(32,32) eg<-sim.data(type="mulmod",N=N) lf<-leafsfirst(eg) ngrid<-4 lf.redu<-treedisc(lf,eg,ngrid=ngrid) lf.plot<-lf.redu lf.plot$level<-c(0,1,1,1,2,2,3) # frame 1 plotbary(lf.plot,ptext=0.2,modelabel=TRUE) # frame 2 plotbary(lf.plot,coordi=2,ptext=0.2,modelabel=TRUE)
# frame 1 a<-1 low<-0.05 end<-2 n<-1000 step<-(end-low)/n curve<-matrix(0,n,2) level<-matrix(0,n,1) for (i in 1:n){ r<-low+i*step phi<-a/r curve[i,1]<-r*cos(phi) curve[i,2]<-r*sin(phi) level[i]<-r } col<-"red" matplot(curve[,1],curve[,2],type="l",xlim=c(-1.5,1.5),ylim=c(-1,1),col=col, xlab="coordinate 1", ylab="coordinate 2") i0<-floor(n/34) low<-0.4 end<-3 n2<-1000 step<-(end-low)/n2 curve2<-matrix(0,n2,2) level2<-matrix(0,n2,1) for (i in 1:n2){ r<-low+i*step phi<-1/r curve2[i,1]<-curve[i0,1]+r*cos(phi)-low*cos(1/low)-i/n2 curve2[i,2]<-curve[i0,2]+r*sin(phi)-low*sin(1/low)-i/n2 level2[i]<-r } col2<-"blue" matplot(curve2[,1],curve2[,2],type="l",add=TRUE,col=col2) i0<-floor(n2/30) low<-0.4 end<-3 n3<-1000 step<-(end-low)/n3 curve3<-matrix(0,n3,2) level3<-matrix(0,n3,1) for (i in 1:n3){ r<-low+i*step phi<-1/r curve3[i,1]<-curve2[i0,1]+r*cos(phi)-low*cos(1/low)-i/n3/2 curve3[i,2]<-curve2[i0,2]+r*sin(phi)-low*sin(1/low)-i/n3/2 level3[i]<-r } col3<-"green" matplot(curve3[,1],curve3[,2],type="l",add=TRUE,col=col3) i0<-floor(n/16) low<-0.4 end<-3 n4<-1000 step<-(end-low)/n4 curve4<-matrix(0,n4,2) level4<-matrix(0,n4,1) for (i in 1:n4){ r<-low+i*step phi<-1/r curve4[i,1]<-curve[i0,1]-r*cos(phi)+low*cos(1/low) curve4[i,2]<-curve[i0,2]-r*sin(phi)+low*sin(1/low) level4[i]<-r } col4<-"orange" matplot(curve4[,1],curve4[,2],type="l",add=TRUE,col=col4) text(0.5,0.9,"A") text(1,0.73,"B") text(1.1,0.3,"C") text(-1,-0.65,"D") # frame 2 # level plot, coordinate 1 coordi<-1 dendat<-curve lenni<-dim(dendat)[1] level<-matrix(0,lenni,1) for (i in 1:lenni) level[i]<-sqrt(sum(dendat[i,]^2)) plot(dendat[,1],level,ylab="level",xlab="coordinate 1",type="l", col=col,xlim=c(-1.5,1.5),log="y") dendat<-curve2 lenni<-dim(dendat)[1] level<-matrix(0,lenni,1) for (i in 1:lenni) level[i]<-sqrt(sum(dendat[i,]^2)) matplot(dendat[,coordi],level,ylab="level",xlab="coordinate 1",type="l", add=TRUE,col=col2,xlim=c(-1.5,1.5))#,ylim=c(0,25)) dendat<-curve3 lenni<-dim(dendat)[1] level<-matrix(0,lenni,1) for (i in 1:lenni) level[i]<-sqrt(sum(dendat[i,]^2)) matplot(dendat[,coordi],level,ylab="level",xlab="coordinate 1",type="l", add=TRUE,col=col3,xlim=c(-1.5,1.5))#,ylim=c(0,25)) dendat<-curve4 lenni<-dim(dendat)[1] level<-matrix(0,lenni,1) for (i in 1:lenni) level[i]<-sqrt(sum(dendat[i,]^2)) matplot(dendat[,coordi],level,ylab="level",xlab="coordinate 1",type="l", add=TRUE,col=col4,xlim=c(-1.5,1.5))#,ylim=c(0,25)) text(0.0,0.72,"A") text(0.1,0.55,"B") text(0.5,0.55,"C") text(-1.2,1.85,"D") # frame 3 # level plot, coordinate 2 coordi<-2 dendat<-curve lenni<-dim(dendat)[1] level<-matrix(0,lenni,1) for (i in 1:lenni) level[i]<-sqrt(sum(dendat[i,]^2)) plot(dendat[,coordi],level,ylab="level",xlab="coordinate 2",type="l", col=col,xlim=c(-1.5,1.5),log="y") dendat<-curve2 lenni<-dim(dendat)[1] level<-matrix(0,lenni,1) for (i in 1:lenni) level[i]<-sqrt(sum(dendat[i,]^2)) matplot(dendat[,coordi],level,ylab="level",xlab="coordinate 2",type="l", add=TRUE,col=col2,xlim=c(-1.5,1.5))#,ylim=c(0,25)) dendat<-curve3 lenni<-dim(dendat)[1] level<-matrix(0,lenni,1) for (i in 1:lenni) level[i]<-sqrt(sum(dendat[i,]^2)) matplot(dendat[,coordi],level,ylab="level",xlab="coordinate 2",type="l", add=TRUE,col=col3,xlim=c(-1.5,1.5))#,ylim=c(0,25)) dendat<-curve4 lenni<-dim(dendat)[1] level<-matrix(0,lenni,1) for (i in 1:lenni) level[i]<-sqrt(sum(dendat[i,]^2)) matplot(dendat[,coordi],level,ylab="level",xlab="coordinate 2",type="l", add=TRUE,col=col4,xlim=c(-1.5,1.5))#,ylim=c(0,25)) text(1.1,1.5,"A") text(0.5,1.5,"B") text(0.27,0.7,"C") text(-0.8,1,"D")
N<-c(32,32) eg<-sim.data(type="mulmod",N=N) lf<-leafsfirst(eg) ngrid<-4 lf.redu<-treedisc(lf,eg,ngrid=ngrid) lf.plot<-lf.redu lf.plot$level<-c(0,1,1,1,2,2,3) info<-matrix("",length(lf.plot$level),1) for (i in 1:length(lf.plot$level)){ info[i]<-as.character(round(lf.plot$volume[i],digits=1)) } plottree(lf.plot,modelabel=TRUE,ptext=0.3,axes=FALSE,info=info,infopos=-0.5, xmarginright=0.1)
N<-c(32,32) eg<-sim.data(type="mulmod",N=N) lf<-leafsfirst(eg) ngrid<-4 lf.redu<-treedisc(lf,eg,ngrid=ngrid) lf.plot<-lf.redu lf.plot$level<-c(0,1,1,1,2,2,3) # frame 1 plotprof(lf.plot,ptext=0.2) # frame 2 plotvolu(lf.plot,ptext=0.2)
n<-200 seed<-5 dendat<-sim.data(type="mulmod",n=n,seed=seed) N<-c(32,32) h<-0.28 pcf<-pcf.kern(dendat,h,N) lst<-leafsfirst(pcf) lsta<-treedisc(lst,pcf,ngrid=120) # frame 1 plottree(lsta,ptext=0.003) # frame 2 plotvolu(lsta,ptext=0.003,modelabel=TRUE)