1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 |
#============ png("#147_R_Lights.png" , units="in", width=15, height=7, res=500) set.seed(1) par(mar=c(0,0,0,0) , bg="black" ) plot( x=c(4,4) , y=c(0,20) , xlim=c(6,40) , ylim=c(2,18) , xaxt="n",yaxt="n", xlab="", ylab="", bty="n" , col=rgb(0.5,0.5,0.5, 0.8) , main="" , type="l" , lwd=0.1 ) # 1 for(i in 1:1000){ a=rnorm(1.5, 3 , n=1) points( x=c(4+a, 4-a) , y=c(0+rnorm(0.3, n=1), 60+rnorm(0.3, n=1)) , col=rgb(0.9,0.4,0.2, 0.5) , type="l", lwd=0.3) } # 2 for(i in 1:1000){ a=rnorm(1.5, 3 , n=1) points( x=c(13+a, 13-a) , y=c(-45+rnorm(0.3, n=1), 25+rnorm(0.3, n=1)) , col=rgb(1,0.5,0.7, 0.4) , type="l", lwd=0.3) } # 3 for(i in 1:1000){ a=rnorm(1.5, 4 , n=1) points( x=c(20+a, 18-a) , y=c(-15+rnorm(0.3, n=1), 30+rnorm(0.3, n=1)) , col=rgb(0.1,0.9,0.2, 0.4) , type="l", lwd=0.3) } # 4 for(i in 1:1000){ a=rnorm(1.5, 4 , n=1) points( x=c(25+a, 26-a) , y=c(-10+rnorm(0.3, n=1), 70+rnorm(0.3, n=1)) , col=rgb(1,0.6,1, 0.4) , type="l", lwd=0.3) } # 5 for(i in 1:1000){ a=rnorm(1.5, 4 , n=1) points( x=c(30+a, 29-a) , y=c(-25+rnorm(0.3, n=1), 75+rnorm(0.3, n=1)) , col= rgb(0.6,0.7,0.2, 0.3) , type="l", lwd=0.7) } # 6 for(i in 1:1000){ a=rnorm(1.5, 4 , n=1) points( x=c(26+a, 43-a) , y=c(-42+rnorm(0.3, n=1), 34+rnorm(0.3, n=1)) , col=rgb(1,0.2,0.5, 0.4) , type="l", lwd=0.7) } legend("bottomright", legend = "© Yan Holtz", bty = "n", text.col = "white", cex=0.9, angle=90) dev.off() |
Not what you are looking for ? Make a new search ![mediatagger]
Leave a Reply
Be the First to Comment!