rm(list=ls())


pdf("DeepNets1.pdf")

plot(c(-1,12),c(1,12),type="n",xlab="",ylab="",xaxt="n",yaxt="n") ##,bty="n")

points(1,4,pch=21,cex=2.5,bg="lightpink")
points(1,6,pch=21,cex=2.5,bg="lightpink")
points(1,8,pch=21,cex=2.5,bg="lightpink")

for(i in 0:4) points(5,2*i+2,pch=21,cex=2.5,bg="red") ## in 1:4

points(9,6,pch=21,cex=2.5,bg="black")

for(i in 2:4){
      for(j in 0:4){
          segments(1,2*i,5,2*j+2,lwd=2,col="purple")
          }
     }	  

for(i in 0:4)segments(5,2*i+2,9,6,lwd=2,col="purple")
#for(i in 1:5)segments(5,2*i+1,10,6.5,lwd=2)

dev.off()




pdf("DeepNets2.pdf")

plot(0,0,xlim=c(0,6),c(0,8),type="n",xaxt="n",yaxt="n",xlab="",ylab="",bty="n")

points(1,3,pch=21,cex=2.5,bg="lightpink")
points(1,4,pch=21,cex=2.5,bg="lightpink")
points(1,5,pch=21,cex=2.5,bg="lightpink")

points(2,2,pch=21,cex=2.5,bg="red")
points(2,3,pch=21,cex=2.5,bg="red")
points(2,4,pch=21,cex=2.5,bg="red")
points(2,5,pch=21,cex=2.5,bg="red")
points(2,6,pch=21,cex=2.5,bg="red")

points(3,1,pch=21,cex=2.5,bg="red")
points(3,2,pch=21,cex=2.5,bg="red")
points(3,3,pch=21,cex=2.5,bg="red")
points(3,4,pch=21,cex=2.5,bg="red")
points(3,5,pch=21,cex=2.5,bg="red")
points(3,6,pch=21,cex=2.5,bg="red")
points(3,7,pch=21,cex=2.5,bg="red")

points(4,1,pch=21,cex=2.5,bg="red")
points(4,2,pch=21,cex=2.5,bg="red")
points(4,3,pch=21,cex=2.5,bg="red")
points(4,4,pch=21,cex=2.5,bg="red")
points(4,5,pch=21,cex=2.5,bg="red")
points(4,6,pch=21,cex=2.5,bg="red")
points(4,7,pch=21,cex=2.5,bg="red")

points(5,4,pch=21,cex=2.5,bg="black")

M1 = matrix(c(1,1,1,3,4,5),ncol=2)
M2 = matrix(c(2,2,2,2,2,2:6),ncol=2)
M3 = matrix(c(3,3,3,3,3,3,3,1:7),ncol=2)
M4 = matrix(c(rep(4,7),1:7),ncol=2)


for(i in 1:3){
     for(j in 1:5){
          segments(M1[i,1],M1[i,2],M2[j,1],M2[j,2],col="purple")
          }
     }


for(i in 1:5){
     for(j in 1:7){
          segments(M2[i,1],M2[i,2],M3[j,1],M3[j,2],col="purple")
          }
     }


for(i in 1:7){
     for(j in 1:7){
          segments(M3[i,1],M3[i,2],M4[j,1],M4[j,2],col="purple")
          }
     }


for(i in 1:7){
     segments(M4[i,1],M4[i,2],5,4,col="purple")
     }

dev.off()