Pular para o conteúdo principal

Cartão de Natal no "R"

Achei bem legal, segue abaixo todo código, basta copiar e colar no seu R

Feliza natal :)

#instando "rgl" package:
install.packages("rgl")

#lendo library
library(rgl)

##Dados

r=1.0
h=5.0
m=h/r
x=seq(0:99)
y=x
X=x*cos(y)
Y=x*sin(y)
Z=-1*(m*x)
Z1=Z+rnorm(length(Z),Z,5)
sub<-matrix(subset(c(X,Y,Z),Z==max(Z)),ncol=3)
trunk<-matrix(c(0.54,0.84,-200,0.54,0.84,-650),byrow=TRUE,ncol=3)
sno<-c(0,rep(15,9))

##Plot card
open3d(FOV=1)
par3d(windowRect=c(100,100,600,600))
bg3d("aquamarine2")
plot3d(X,Y,Z,col="darkgreen",type="l",box=FALSE,axes=FALSE,lwd=10,zlim=c(min(Z),0),xlab="",ylab="",zlab="",top=TRUE,border=TRUE)
plot3d(X,Y,Z,col=c("yellow","red"),add=TRUE,type="s",radius=3)
points3d(sub,col="orange",size=10)
lines3d(trunk,col="brown",lwd=20)
text3d(matrix(c(0.54,-100,300),ncol=3),text="Happy Holidays",font=5,cex=1.9,color="darkred",adj=c(0.5,1))
text3d(matrix(c(0.54,0.84,-650),ncol=3),text="and",font=5,cex=1.5,color="darkred",adj=c(0.5,1))
text3d(matrix(c(0.54,75,-780),ncol=3),text="a Happy New Year",font=5,cex=1.9,color="darkred",adj=c(0.5,1))
mat<-par3d("scale")

##Play card
for(i in 1:10){
par3d(scale=mat,ignoreExtent=TRUE)
points3d(matrix(c(sample(-160:160,sno[i]),sample(-160:160,sno[i]),sample(-500:150,sno[i])),ncol=3),col="white",size=3)
play3d(spin3d(axis=c(0,0,1)),duration=1)
if(i==10){points3d(sub,col="darkorange",size=17)
plot3d(X,Y,Z,col=c("yellow","red"),type="s",radius=4,add=TRUE) } }
text3d(matrix(c(0,140,-900),ncol=3),text="From the EEB & Flow",font=1,cex=1.0,color="darkgreen",adj=c(0.5,1))

Comentários

  1. Cara, achei bem legal tb!!!
    Fiz a pouco no meu note aqui!
    Um abração!

    ResponderExcluir
  2. Abraço Menezes

    Pelo que vi tem os créditos no fim da programação, mas retirei do grupo de estatística do Facebook e da R-list

    ResponderExcluir

Postar um comentário

Postagens mais visitadas deste blog

LISTA FONTES DE PESQUISA ÚTEIS E CONFIÁVEIS

http://novo.periodicos.capes.gov.br/?option=com_pnews&component=NewsShow&view=pnewsnewsshow&cid=128&mn=0 (Novo portal periódicos CAPES) http://www.enap.gov.br/index.php?option=com_content&task=view&id=252&Itemid=65 (ENAP - Escola Nacional de Administração Pública) http://www.scielo.br/scielo.php?lng=pt (BASE DE DADOS SCIELO) http://www.periodicos.capes.gov.br/portugues/index.jsp (PERIÓDICOS CAPES) http://acessolivre.capes.gov.br/ http://www.dominiopublico.gov.br/pesquisa/PesquisaObraForm.jsp (DOMÍNIO PÚBLICO) www.ufrgs.br http://www.ea.ufrgs.br/ www.culturaacademica.com.br (UNESP) http://www.armazemmemoria.com. br/AjudaArmazem. aspx (história) http://www.4shared. com/dir/18704839 /b1b64358/ anthropology. html www.bvce.org (esse dá acesso a várias bibliotecas virtuais) http://bibliotecadearqueologia.blogspot.com/ http://www.ufgd.edu.br/historiaemreflexao/julho_dez_2008/edicoes-anteriores

Leitura da madrugada - Métodos para modelos aditivos não-paramétricos

O título é "Finite sample performance of kernel-based regression methods for non-parametric additive models under common bandwidth selection criterion" de CARLOS MARTINS-FILHO e KEYANG, preciso ler o artigo,para além do conhecimento, responder uma questão da lista de exercícios de econometria não-paramétrica.

Exemplo de estimação de VAR no "R"

A idéia nesse post é "ajudar" na estimativa de um Vetor auto-regressivo usando o programa "R"....OBS: os comandos estão escritos em vermelho. Os dados utilizados foram coletados no sistema de gerenciador de séries temporais do Bacen (link barra lateral) e no Ipeadata (link barra lateral). Para facilitar o exercício salvei o documento em minha conta no google docs, é só clicar a baixar a tabela pro teste, o nome do arquivo é "macro". A períodicidade é mensal, de jan de 2000 a dez 2010....as variáveis já estão estacionárias (para PIB, SELIC (dessazonalizado) e CÂMBIO aplicamos a primeira diferença), IBOVESPA (Ibovespa - variação percentual mensal) e IGP-M são estacionãrias em nível. O primeiro passo é importar os dados para o programa.....eu sempre tenho costume de, dentro do R, ir em "Arquivo" - "Mudar dir.." e escolher a pasta para importar os dados salvos.....no meu caso, salvo sempre dentro da própria pasta do R....em "