Merry Christmas tree & Happy New Year in R
Nishant Upadhyay
Eco dept.Univ of Pune/Infosysnishantup@gmail.com
2017-01-02
Merry-christmas-2016-with-r
pacman::p_load(plot3D)
trunk <- list(
x0=0, y0=0, z0=0,
x1=0, y1=0, z1=6,
extension=NULL,
branch1=NULL, branch2=NULL,
branch3=NULL, branch4=NULL,
depth=1,
lwd=1,
stub_color="brown"
)
draw_tree <- function(tree) {
#function to draw the tree recursively
if (is.null(tree)) return()
with(tree,{
segments3D(x0,y0,z0,x1,y1,z1,col=stub_color,
lwd=lwd,add=TRUE)
draw_tree(branch1)
draw_tree(branch2)
draw_tree(branch3)
draw_tree(extension)
})
}
extend <- function(tree) {
#creates an extension of the main branch
#in the same direction vector u=p1-p0
with(tree,{
growth_factor <- runif(1,min=0.9,max=1)
return(list(x0=x1,y0=y1,z0=z1,
x1=x1+growth_factor*(x1-x0),
y1=y1+growth_factor*(y1-y0),
z1=z1+growth_factor*(z1-z0),
extension=NULL,
branch1=NULL, branch2= NULL,
branch3=NULL, branch4= NULL,
depth=depth*0.9,
lwd=1,
stub_color="darkgreen"))
})
}
create_branch <- function(tree,c1,c2) {
#this function returns a branch for tree
#c1 and c2 are the components for the vector
#base {v1,v2}. v1 and v2 are two orthonrmal
#vectors, both perpedicular to the direction
#vector u=p1-p0 (p1 and p0 are the endpoints
#of the tree trunk).
#where_at is a scalar value in (0,1), indicating
#where along p0=(x0,y0,z0) and p1=(x1,y1,z1)
#the branch will grow.
with(tree,{
#vector u=(x1-x0,y1-y0,z1-z0) gives the
#direction of the tree trunk.
#stub_length is the length of that trunk.
u <- c(x1-x0,y1-y0,z1-z0)
stub_length <- sqrt(sum(u*u))
#growth_factor is how long the branch
#will be with respect to stub_length
growth_factor <- runif(1,min=0.7,max=0.8)
where_at <- runif(1,min=0.4,max=0.8)
#create two perpendicular vectors to u
#and set their norms to 1
i <- which(u!=0)[1]
oi <- setdiff(1:3,i)
v1 <- rep(0,3); v1[i] <- -u[oi[1]]; v1[oi[1]] <- u[i]
v2 <- rep(0,3); v2[i] <- -u[oi[2]]; v2[oi[2]] <- u[i]
v1 <- v1/sqrt(sum(v1*v1)); v2 <- v2/sqrt(sum(v2*v2))
#vector v is a linear combination of v1 and v2,
#hence it lies on the plane perpendicular to u
v <- v1*c1+v2*c2
#beta is the angle of separation between the
#branch and the trunk. The angle is larger
#for the "deeper" parts of the tree.
beta <- runif(1,min=0.7*depth,max=0.9*depth)
#calculate the scalar k for multiplying with v
#so that u+kv, the branch, will form a beta
#angle with the trunk.
k <- tan(beta)*sqrt(sum(u*u))/sqrt(sum(v*v))
#new_u is the direction vector of the branch,
#whose length will be the length of the branch
new_u <- u+k*v
new_u <- new_u / sqrt(sum(new_u*new_u))*
stub_length*growth_factor
#the new (x0,y0,z0) point of the branch
#is where along the main trunk the branch begins.
new_x0 <- x0+(x1-x0)*where_at
new_y0 <- y0+(y1-y0)*where_at
new_z0 <- z0+(z1-z0)*where_at
#the new (x1,y1,z1) point of the branch is
#its ending point.
new_x1 <- new_x0+new_u[1]
new_y1 <- new_y0+new_u[2]
new_z1 <- new_z0+new_u[3]
list(
x0=new_x0,y0=new_y0,z0=new_z0,
x1=new_x1,y1=new_y1,z1=new_z1,
extension=NULL,
branch1=NULL, branch2= NULL,
branch3=NULL, branch4= NULL,
depth=depth*0.8,
lwd=1,
stub_color="darkgreen")
})
}
grow_tree <- function(tree) {
if (is.null(tree) ) return(NULL)
tree$lwd <- tree$lwd*1.4
if (tree$lwd>2.0) tree$stub_color <- 'brown4'
if (is.null(tree$extension)) {
tree$extension <- extend(tree)
if (tree$depth<0.5) {
#make two branches at an angle of more than 160 degrees
rot1 <- runif(1,min=0,max=2*pi)
rot2 <- (runif(1,min=8*pi/9,max=10*pi/9)+rot1)%%(2*pi)
tree$branch1 <- create_branch(tree,cos(rot1),sin(rot1))
tree$branch2 <- create_branch(tree,cos(rot2),sin(rot2))
} else {
#make 3 branches with 120 degree angles between them and
#rotate the entire branch set.
rot <- runif(1,min=0,max=2*pi)
tree$branch1 <- create_branch(tree,-sin(rot),cos(rot))
tree$branch2 <- create_branch(tree,
sqrt(3)/2*cos(rot)+sin(rot)/2,
sqrt(3)/2*sin(rot)-cos(rot)/2)
tree$branch3 <- create_branch(tree,
-sqrt(3)/2*cos(rot)+sin(rot)/2,
-sqrt(3)/2*sin(rot)-cos(rot)/2)
}
} else {
tree$extension <- grow_tree(tree$extension)
tree$branch1 <- grow_tree(tree$branch1)
tree$branch2 <- grow_tree(tree$branch2)
if (tree$depth>=0.5)
tree$branch3 <- grow_tree(tree$branch3)
}
return(tree)
}
create_ornaments <- function(tree) {
if (is.null(tree)) return()
po <- (1-tree$depth)^6
ornament <- sample(c(T,F),size=1,prob=c(po,1-po))
co <- sample(c("red","darkgoldenrod4"),size=1,
prob=c(0.6,0.4))
if (ornament)
ornaments <<- rbind(ornaments,
data.frame(x=tree$x1,y=tree$y1,z=tree$z1,color=co))
create_ornaments(tree$branch1)
create_ornaments(tree$branch2)
create_ornaments(tree$branch3)
create_ornaments(tree$extension)
}
create_lights1 <- function(tree) {
if (is.null(tree)) return()
po <- (1-tree$depth)^4
light <- sample(c(T,F),size=1,prob=c(po,1-po))
if (light)
lights1 <<- rbind(lights1,
data.frame(x=tree$x1,y=tree$y1,z=tree$z1))
create_lights1(tree$branch1)
create_lights1(tree$branch2)
create_lights1(tree$branch3)
create_lights1(tree$extension)
}
create_lights2 <- function(tree) {
if (is.null(tree)) return()
po <- (1-tree$depth)^4
light <- sample(c(T,F),size=1,prob=c(po,1-po))
if (light)
lights2 <<- rbind(lights2,
data.frame(x=(tree$x1+tree$x0)/2,
y=(tree$y1+tree$y0)/2,
z=(tree$z1+tree$z0)/2))
create_lights2(tree$branch1)
create_lights2(tree$branch2)
create_lights2(tree$branch3)
create_lights2(tree$extension)
}
draw_ornaments <- function() {
with(ornaments,
{points3D(x=x,y=y,z=z,
pch=19,cex=1.2,col=as.character(color),
colkey=FALSE,add=TRUE)})
}
draw_lights1 <- function() {
with(lights1,
{points3D(x=x,y=y,z=z,pch="+",cex=0.8,col="white",
colkey=FALSE,add=TRUE)})
}
draw_lights2 <- function() {
with(lights2,
{points3D(x=x,y=y,z=z,pch="+",cex=0.8,col="yellow",
colkey=FALSE,add=TRUE)})
}
set.seed(20161224)
pine_tree <- trunk
for (i in 1:5) pine_tree <- grow_tree(pine_tree)
ornaments <- NULL;
lights1 <- NULL;
lights2 <- NULL;
create_ornaments(pine_tree)## NULLcreate_lights1(pine_tree)## NULLcreate_lights2(pine_tree)## NULL#png("tree%02d.png")
for (i in 9:10) {
perspbox(x=c(-25,25),y=c(-25,25),z=c(0,40),bty="n",
phi=8,theta=i*10,col="white",alpha=0)
draw_tree(pine_tree)
draw_ornaments()
switch((i%%4)+1,{},
{draw_lights1()},
{draw_lights2()},
{draw_lights1(); draw_lights2()})
}
#graphics.off()CHRISTMAS TREE WITH GGPLOT
rm(list = ls())
library(ggplot2)
# create data
x <- c(8,7,6,7,6,5,6,5,4,5,4,3,4,3,2,3,2,1,0.5,0.1)
dat1 <- data.frame(x1 = 1:length(x), x2 = x)
dat2 <- data.frame(x1 = 1:length(x), x2 = -x)
dat1$xvar <- dat2$xvar <- NA
dat1$yvar <- dat2$yvar <- NA
dat1$siz <- dat2$siz <- NA
dat1$col <- dat2$col <- NA
# set threshold for christmas balls
dec_threshold = -0.5
# create random places, sizes and colors for christmas balls
set.seed(2512)
for (row in 1:nrow(dat1)){
if (rnorm(1) > dec_threshold){
dat1$xvar[row] <- row
dat1$yvar[row] <- sample(1:dat1$x2[row]-1,1)
dat1$siz[row] <- runif(1,0.5,1.5)
dat1$col[row] <- sample(1:5, 1)
}
if (rnorm(1) > dec_threshold){
dat2$xvar[row] <- row
dat2$yvar[row] <- sample(1:dat2$x2[row],1)
dat2$siz[row] <- runif(1,0.5,1.5)
dat2$col[row] <- sample(1:5, 1)
}
}
# plot the christmas tree
ggplot() +
geom_bar(data = dat1, aes(x=x1, y=x2),stat = "identity", fill = '#31a354') +
geom_bar(data = dat2, aes(x=x1, y=x2),stat = "identity", fill = '#31a354') +
geom_point(data = dat1,aes(x = xvar, y = yvar, size = siz, colour = as.factor(col)) ) +
geom_point(data = dat2,aes(x = xvar, y = yvar, size = siz, colour = as.factor(col)) ) +
coord_flip() + theme_minimal()+ theme(legend.position="none",
axis.title.x=element_blank(),
axis.text.x=element_blank(),
axis.ticks.x=element_blank(),
axis.title.y=element_blank(),
axis.text.y=element_blank(),
axis.ticks.y=element_blank()) +
ggtitle('The Analytics Lab wishes you a Merry Christmas')## Warning: Removed 7 rows containing missing values (geom_point).## Warning: Removed 4 rows containing missing values (geom_point).
Quick R Christmas tree
library(ggplot2)
dataset <- c(9,8,7,6,5,4,4.5,4,3,2,1,0.8,0.2)
df <- data.frame(group = rep(c("l","d"),
each=length(dataset)),
x = 1:length(dataset),
y = c(dataset, dataset*-1))
ggplot() +
geom_bar(data = df, aes(x=x, y=y),stat = "identity", fill = '#00A650',width=.8) + coord_flip() + ggtitle("Happy New Year 2017") + theme_void()
FIREWORKS (IN R)
library(ggplot2)
rm(list = ls())
# First rocket
t1 = 1:30 # number of points
size1 = 2 # size of rocket
r1_center = c(-4,10) # end center
rocket1 = data.frame(x1 = rep(r1_center[1], length(t1))
, x2 = size1*cos(t1) + r1_center[1]
, y1 = rep(r1_center[2], length(t1))
, y2 = size1*sin(t1) + r1_center[2]
)
# Second rocket
t2 = 1:44
size2 = 3
r2_center = c(3,12)
rocket2 = data.frame(x1 = rep(r2_center[1], length(t2))
, x2 = size2*cos(t2) + r2_center[1]
, y1 = rep(r2_center[2], length(t2))
, y2 = size2*sin(t2) + r2_center[2]
)
# Third rocket
t3 = 1:44
size3 = 4
r3_center = c(-2,17)
rocket3 = data.frame(x1 = rep(r3_center[1], length(t3))
, x2 = size3*cos(t3) + r3_center[1]
, y1 = rep(r3_center[2], length(t3))
, y2 = size3*sin(t3) + r3_center[2]
)
# Fourth rocket
t4 = 1:44
size4 = 4
r4_center = c(4,20)
rocket4 = data.frame(x1 = rep(r4_center[1], length(t4))
, x2 = size4*cos(t4) + r4_center[1]
, y1 = rep(r4_center[2], length(t4))
, y2 = size4*sin(t4) + r4_center[2]
)
# Plot fireworks
ggplot() + theme(panel.background = element_rect(fill = '#252525', colour = '#252525'), panel.grid.major = element_blank(), panel.grid.minor = element_blank(),axis.line=element_blank(),axis.text.x=element_blank(),axis.text.y=element_blank(),axis.ticks=element_blank(),axis.title.x=element_blank(),axis.title.y=element_blank(),legend.position="none") +
# first rocket
geom_point(aes(x = x2, y = y2), data = rocket1, shape = 4, colour = '#74add1')+ geom_curve(aes(x = x1, y = y1, xend = x2, yend = y2 ), curvature = 0.1, data = rocket1, colour = '#74add1') + geom_curve(aes(x = 0, y = -7, xend = r1_center[1], yend = r1_center[2] ), curvature = 0.1, colour = '#878787') +
# second rocket
geom_point(aes(x = x2, y = y2), data = rocket2, shape = 4, colour = '#fed976')+ geom_curve(aes(x = x1, y = y1, xend = x2, yend = y2), curvature = 0.1, data = rocket2, colour = '#fed976') + geom_curve(aes(x = 0, y = -7, xend = r2_center[1], yend = r2_center[2] ), curvature = -0.1, colour = '#878787') +
# third rocket
geom_point(aes(x = x2, y = y2), data = rocket3, shape = 4,
colour ='#fa9fb5')+ geom_curve(aes(x = x1, y = y1, xend = x2, yend = y2 ), curvature = 0.1, data = rocket3, colour = '#fa9fb5') + geom_curve(aes(x = 0, y = -7, xend = r3_center[1], yend = r3_center[2] ), curvature = 0.1, colour = '#878787') +
# fouth rocket
geom_point(aes(x = x2, y = y2), data = rocket4, shape = 4, colour = '#addd8e')+ geom_curve(aes(x = x1, y = y1, xend = x2, yend = y2 ), curvature = 0.1, data = rocket4, colour = '#addd8e') + geom_curve(aes(x = 0, y = -7, xend = r4_center[1], yend = r4_center[2] ), curvature = -0.2, colour = '#878787') +
# title
ggtitle('Nishant wishes you a Happy and a Prosperous New Year') 
# save
#ggsave('Happy New Year.png', units = 'cm', width = 15, height = 20)Happy New Year animations-in-r-using-plotly
library(plotly)
#devtools::install_github("ropensci/plotly")
rm(list = ls())
gc()## used (Mb) gc trigger (Mb) max used (Mb)
## Ncells 612107 32.7 2110301 112.8 2637877 140.9
## Vcells 782323 6.0 30953105 236.2 40236841 307.0# Options for plotting ----
x <- 0.2
y <- 0.72
speed <- 250
nbkdrops <- 100
# Colorset for plot
# See http://colorhunt.co/
cols <- c("#FFC85B", "#379956","#234C63")
ncolors <- length(cols)
# Function to create random points by adding jitter to ----
# a starting set of points
n <- 1000 # Number of points
# Starting template
bkdrop.x <- runif(n, min = 0, max = 1)
bkdrop.y <- runif(n, min = 0, max = 1)
# Function Definition
bkdrop <- function(n = 1000, amount = 0.005){
x <- jitter(bkdrop.x, amount = amount)
y <- jitter(bkdrop.y, amount = amount)
df <- data.frame(x, y)
return(df)
}
# Make backdrops ----
# Each call to the backdrop function is a separate frame
# Number of frames is controlled by nbkdrops
bkdrop.df <- data.frame()
for(i in 1:nbkdrops){
temp <- bkdrop()
temp <- data.frame(temp, frame = i, color = sample(1:ncolors, size = nrow(temp), replace = T))
bkdrop.df <- rbind(bkdrop.df, temp)
}
# Make back lights ----
# Coordinates for backlight rectangles
# Will be plotted as line segments
bklight.x <- c(0.28, 0.18, 0.48)
bklight.y <- c(0.42, 0.62, 0.65)
bklight.xend <- c(0.63, 0.50, 0.75)
bklight.yend <- c(0.42, 0.62, 0.65)
# Function to create a dataframe containing coordinates, frame and
# color of each backlight segment
makebklight <- function(id){
bklight <- data.frame()
for(i in 1:nbkdrops){
temp <- data.frame(x = bklight.x[id],
y = bklight.y[id],
xend = bklight.xend[id],
yend = bklight.yend[id],
frame = i,
color = sample(1:ncolors, size = 1))
bklight <- rbind(bklight, temp)
}
return(bklight)
}
# Create backlight segments
bklight1 <- makebklight(1)
bklight2 <- makebklight(2)
bklight3 <- makebklight(3)
# Initialize colors for first frame
bklight1$color[1] <- 1
bklight2$color[1] <- 2
bklight3$color[1] <- 3
# Plot !! ----
p <- plot_ly(height = 800, width = 1024,
colors = cols,
frame = ~frame,
x = ~x,
y = ~y,
color = ~factor(color)) %>%
# Backdrop
add_markers(data = bkdrop.df,
opacity = 0.8,
marker = list(symbol = "star", size = 8),
hoverinfo = "none") %>%
# Add segments (for back lighting)
add_segments(data = bklight1,
xend = ~xend, yend = ~yend,
line = list(width = 150)) %>%
add_segments(data = bklight2,
xend = ~xend, yend = ~yend,
line = list(width = 150)) %>%
add_segments(data = bklight3,
xend = ~xend, yend = ~yend,
line = list(width = 150)) %>%
# Animation options
# See https://cpsievert.github.io/plotly_book/key-frame-animations.html
animation_opts(speed, easing = "linear", transition = 0) %>%
animation_button(x = 1, xanchor = "right", y = 1, yanchor = "bottom") %>%
animation_slider(hide = T) %>%
# Layout, annotations and shapes
layout(
showlegend = F,
xaxis = list(title = "", showgrid = F, zeroline = F, showticklabels = F, range = c(0, 1)),
yaxis = list(title = "", showgrid = F, zeroline = F, showticklabels = F, range = c(0, 1)),
annotations = list(
# For shadow
list(xref = "paper", yref = "paper",
xanchor = "left", yanchor = "top",
x = x + 0.002, y = y + 0.002,
showarrow = F,
text = "Happy New<br>Year !",
font = list(size = 100, family = "Times New Roman",
color = "black")),
list(xref = "paper", yref = "paper",
xanchor = "left", yanchor = "top",
x = x + 0.003, y = y + 0.003,
showarrow = F,
text = "Happy New<br>Year !",
font = list(size = 100, family = "Times New Roman",
color = "black")),
list(xref = "paper", yref = "paper",
xanchor = "left", yanchor = "top",
x = x + 0.004, y = y + 0.004,
showarrow = F,
text = "Happy New<br>Year !",
font = list(size = 100, family = "Times New Roman",
color = "black")),
# Actual
list(xref = "paper", yref = "paper",
xanchor = "left", yanchor = "top",
x = x, y = y,
showarrow = F,
text = "Happy New<br>Year !",
font = list(size = 100, family = "Times New Roman",
color = "#ff6666"))
),
shapes = list(
# Border
list(xref = "paper", yref = "paper",
x0 = 0, y0 = 0,
x1 = 1, y1 = 1,
type = "rect",
line = list(width = 10, color = cols[1])),
list(xref = "paper", yref = "paper",
x0 = 0.01, y0 = 0.01,
x1 = 0.99, y1 = 0.99,
type = "rect",
line = list(width = 10, color = cols[2])),
list(xref = "paper", yref = "paper",
x0 = 0.02, y0 = 0.02,
x1 = 0.98, y1 = 0.98,
type = "rect",
line = list(width = 10, color = cols[3])),
# Black outline
list(xref = "plot", yref = "plot",
path = "
M 0.50 0.53
L 0.50 0.50
L 0.18 0.50
L 0.18 0.73
L 0.48, 0.73",
type = "path",
line = list(width = 7, color = "black")),
list(xref = "plot", yref = "plot",
path = "
M 0.50 0.535
L 0.48 0.535
L 0.48 0.77
L 0.75 0.77
L 0.75 0.535
Z",
type = "path",
line = list(width = 7, color = "black")),
list(xref = "plot", yref = "plot",
path = "
M 0.28 0.5
L 0.28 0.31
L 0.63 0.31
L 0.63 0.535",
type = "path",
line = list(width = 7, color = "black"))
)
)
pdevtools::session_info()## Session info --------------------------------------------------------------## setting value
## version R version 3.3.2 (2016-10-31)
## system x86_64, mingw32
## ui RTerm
## language (EN)
## collate English_India.1252
## tz Asia/Calcutta
## date 2017-01-02## Packages ------------------------------------------------------------------## package * version date source
## assertthat 0.1 2013-12-06 CRAN (R 3.3.0)
## backports 1.0.4 2016-10-24 CRAN (R 3.3.1)
## base64enc 0.1-3 2015-07-28 CRAN (R 3.3.0)
## Cairo 1.5-9 2015-09-26 CRAN (R 3.3.2)
## colorspace 1.3-2 2016-12-14 CRAN (R 3.3.2)
## colourpicker 0.3 2016-12-05 CRAN (R 3.3.2)
## crosstalk 1.0.0 2016-12-21 CRAN (R 3.3.2)
## DBI 0.5-1 2016-09-10 CRAN (R 3.3.1)
## devtools 1.12.0 2016-06-24 CRAN (R 3.3.1)
## digest 0.6.10 2016-08-02 CRAN (R 3.3.1)
## dplyr 0.5.0 2016-06-24 CRAN (R 3.3.1)
## evaluate 0.10 2016-10-11 CRAN (R 3.3.1)
## ggplot2 * 2.2.1 2016-12-30 CRAN (R 3.3.2)
## gtable 0.2.0 2016-02-26 CRAN (R 3.2.3)
## htmltools 0.3.5 2016-03-21 CRAN (R 3.2.4)
## htmlwidgets 0.8 2016-11-09 CRAN (R 3.3.2)
## httpuv 1.3.3 2015-08-04 CRAN (R 3.2.1)
## httr 1.2.1 2016-07-03 CRAN (R 3.3.1)
## jsonlite 1.2 2016-12-31 CRAN (R 3.3.2)
## knitr 1.15.1 2016-11-22 CRAN (R 3.3.2)
## labeling 0.3 2014-08-23 CRAN (R 3.2.0)
## lazyeval 0.2.0 2016-06-12 CRAN (R 3.3.0)
## magrittr 1.5 2014-11-22 CRAN (R 3.2.1)
## memoise 1.0.0 2016-01-29 CRAN (R 3.2.3)
## mime 0.5 2016-07-07 CRAN (R 3.3.1)
## miniUI 0.1.1 2016-01-15 CRAN (R 3.2.4)
## misc3d 0.8-4 2013-01-25 CRAN (R 3.3.0)
## munsell 0.4.3 2016-02-13 CRAN (R 3.2.3)
## pacman 0.4.1 2016-03-30 CRAN (R 3.2.4)
## plot3D * 1.1 2016-01-13 CRAN (R 3.3.0)
## plotly * 4.5.6.9000 2017-01-02 Github (ropensci/plotly@9b8e28c)
## plyr 1.8.4 2016-06-08 CRAN (R 3.3.0)
## purrr 0.2.2 2016-06-18 CRAN (R 3.3.0)
## R6 2.2.0 2016-10-05 CRAN (R 3.3.1)
## Rcpp 0.12.8 2016-11-17 CRAN (R 3.3.2)
## rmarkdown 1.3 2016-12-21 CRAN (R 3.3.2)
## rprojroot 1.1 2016-10-29 CRAN (R 3.3.2)
## scales 0.4.1 2016-11-09 CRAN (R 3.3.2)
## shiny 0.14.2 2016-11-01 CRAN (R 3.3.2)
## stringi 1.1.2 2016-10-01 CRAN (R 3.3.1)
## stringr 1.1.0 2016-08-19 CRAN (R 3.3.1)
## tibble 1.2 2016-08-26 CRAN (R 3.3.1)
## tidyr 0.6.0 2016-08-12 CRAN (R 3.3.1)
## viridisLite 0.1.3 2016-03-12 CRAN (R 3.2.4)
## withr 1.0.2 2016-06-20 CRAN (R 3.3.0)
## xtable 1.8-2 2016-02-05 CRAN (R 3.2.3)
## yaml 2.1.14 2016-11-12 CRAN (R 3.3.2)Reference:
https://unamatematicaseltigre.blogspot.com/2016/12/merry-christmas-2016-with-r.html
http://www.theanalyticslab.nl/2016/12/25/christmas-tree-with-ggplot/
https://tomaztsql.wordpress.com/2016/12/30/quick-r-christmas-tree/