Ivan, não sei se é bem isto que você se refere, mas vc pode usar o ffmpeg para fazer um "filme" com vários gráficos. Neste exemplo [1] ele mostra, na última linha, o comando utilizado para gerar o vídeo. [1] http://www.drewconway.com/zia/?p=2800 []s, --- Fernando Mayer Universidade Federal de Santa Catarina - UFSC Departamento de Ecologia e Zoologia - ECZ/CCB URL: http://fernandomayer.github.com e-mail: fernandomayer [@] gmail.com 2012/8/28 Ivan Bezerra Allaman <ivanalaman@yahoo.com.br>:
Bom dia senhores!
Estou fazendo uma ilustração animada para demonstrar que a média é influenciada por valores extremos com o seguinte código:
png(file='exemplos/med0%1d.png',width=500,height=250) x1 <- list() medias <- list() for(i in 1:10){ x <- c(7,10,13.4,15.3,17.6,19.5,20.6,21.1) y <- rep(0,8) x1[[i]] <- c(x,i+21.1) medias[[i]] <- mean(x1[[i]]) plot(x,y,ylim=c(0,5),xlim=c(0,35),ylab='',xlab='',axes=F) axis(1) par(xpd=TRUE) points(mean(x),-0.2,col='blue',lwd=10,pch=2) points(x1[[i]][9],0,col='red',lwd=7,cex=0.1) points(medias[[i]],-0.2,col='blue',lwd=10,pch=2) } dev.off()
png(file='exemplos/med1%1d.png',width=500,height=250) x2 <- list() medias2 <- list() for(j in seq(10,1,-1)){ x2[[j]] <- c(12.3-j,x) medias2[[j]] <- mean(x2[[j]]) plot(x,y,ylim=c(0,5),xlim=c(0,35),ylab='',xlab='',axes=F) axis(1) par(xpd=TRUE) points(mean(x),-0.2,col='blue',lwd=10,pch=2) points(x1[[j]][9],0,col='red',lwd=7,cex=0.1) points(medias[[j]],-0.2,col='blue',lwd=10,pch=2) } dev.off()
png(file='exemplos/med2%1d.png',width=500,height=250) for(j in 1:10){ x2[[j]] <- c(12.3-j,x) medias2[[j]] <- mean(x2[[j]]) plot(x,y,ylim=c(0,5),xlim=c(0,35),ylab='',xlab='',axes=F) axis(1) par(xpd=TRUE) points(mean(x),-0.2,col='blue',lwd=10,pch=2) points(x2[[j]][1],0,col='red',lwd=7,cex=0.1) points(medias2[[j]],-0.2,col='blue',lwd=10,pch=2) } dev.off()
É possível juntar isso em um único código?
Desde já grato por qualquer ajuda.
(S,f,P) Allaman
\begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.