FELIZ NAVIDAD Y PRÓSPERO AÑO NUEVO

Este año, en Todo Econometría ha sido una verdadera aventura con la ayuda de los lectores, amigos y seguidores de nuestras redes sociales.

Nuestra participación en varias actividades de consultoría y clases en diferentes programas de Maestrías han hecho de este sitio y del que se dirige a ustedes muy feliz de poder ayudar a conseguir objetivos que nuestros clientes y nuestros alumnos se planteaban  como retos inalcanzables antes de conocernos.

A todos ustedes familia, amig@s, seguid@res, fans, les agradecemos por sus likes, sus comentarios que leemos y contestamos a todos lo mas antes posible.

Este cierre de año les dejo un pequeño programa modificado por mi para estas navidades para el paquete R.

FELIZ NAVIDAD Y UN 2017 PRÓSPERO

rplot_treerxmas

part
branch1=NULL,branch2=NULL,extend=NULL,
lwd=1,depth=0,col=’springgreen’)

par(mfrow=c(1,1),mar=c(5, 4, 4, 2) + 0.1)
segplot
if (is.null(tree)) return()
segments(tree$x0,tree$y0,tree$x1,tree$y1,
col=tree$col,
lwd=tree$lwd)
segplot(tree$branch1)
segplot(tree$branch2)
segplot(tree$extend)
}
#segplot(part)

grow
if (is.null(tree) ) return(NULL)

tree$lwd=tree$lwd*1.2
if (tree$lwd>2.5) tree$col <- ‘brown’
if (is.null(tree$extend)) {
tree$extend
x0=tree$x1,
y0=tree$y1,
x1=rnorm(1,1,.03)*(2*tree$x1-tree$x0),
y1=(rnorm(1,.98,.02)+.02*(tree$x1==tree$x0))*(2*tree$y1-tree$y0),
branch1=NULL,
branch2=NULL,
extend=NULL,
lwd=1,
depth=tree$depth,
col=tree$col
)
length=sqrt((tree$x1-tree$x0)^2 + (tree$y1-tree$y0)^2)
angle
branch
x0=(tree$x1+tree$x0)/2,
y0=(tree$y1+tree$y0)/2,
branch1=NULL,
branch2=NULL,
extend=NULL,
lwd=1,
depth=tree$depth,
col=tree$col
)
shift
branch$x0
branch$y0
length=length*rnorm(1,.5,.05)
co
branch$x1

branch$y1
tree$branch1
branch$x0
branch$y0
co
branch$x1
branch$y1
tree$branch2
} else {
tree$branch1
tree$branch2
tree$extend
}
tree$depth
if (tree$depth>2) tree$col <- ‘green’
if (tree$depth>4) tree$col <- ‘darkgreen’
if (tree$depth>6) tree$col <- ‘brown’

tree
}
tree
for (i in 1:9) tree
par(mar=c(0,0,0,0))
plot(x=c(-3,3),y=c(0,9),type=’n’,axes=FALSE,xlab=”,ylab=”)
segplot(tree)
# Add some presents
xPres = runif(10,-4.5,4.5)
xWidth = runif(10,0.1,0.5)
xHeight=runif(10,0,1)
for(i in 1:10){
rect(xPres[i]-xWidth[i],0,xPres[i]+xWidth[i],xHeight[i],col=sample(c(“blue”,”red”),size=1))
rect(xPres[i]-0.2*xWidth[i],0,xPres[i]+0.2*xWidth[i],xHeight[i],col=sample(c(“gold”,”grey87″),size=1))
}

points(x=runif(4,-2,2),y=rep(2,4),col=sample(c(“blue”,”red”),size=4,replace=T),cex=3,pch=19)
points(x=runif(4,-1,1),y=rep(2.5,4),col=sample(c(“blue”,”red”),size=4,replace=T),cex=3,pch=19)
points(x=runif(3,-2,2),y=rep(4,3),col=sample(c(“blue”,”red”),size=4,replace=T),cex=3,pch=19)
points(x=runif(4,-1,1),y=rep(6.25,4),col=sample(c(“blue”,”red”),size=4,replace=T),cex=3,pch=19)
points(0,8.5,pch=8,cex=5,col=”gold”,lwd=3)

install.packages(“audio”)
library(“audio”)
xmas_file
download.file(“https://sites.google.com/site/pocketecoworld/merrychristmas1.wav&#8221;, xmas_file, mode = “wb”)
xmas
play(xmas)

Gracias a la inspiración de este programa a stackoverflow, y a jtleek/gist:4369771.

 

P.D. Si te gustó no olvides darle al like y a compartir, gracias!

Anuncios

Responder

Introduce tus datos o haz clic en un icono para iniciar sesión:

Logo de WordPress.com

Estás comentando usando tu cuenta de WordPress.com. Cerrar sesión / Cambiar )

Imagen de Twitter

Estás comentando usando tu cuenta de Twitter. Cerrar sesión / Cambiar )

Foto de Facebook

Estás comentando usando tu cuenta de Facebook. Cerrar sesión / Cambiar )

Google+ photo

Estás comentando usando tu cuenta de Google+. Cerrar sesión / Cambiar )

Conectando a %s

Crea un sitio web o blog en WordPress.com

Subir ↑

A %d blogueros les gusta esto: