Ajuda em ilustração animada!
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}
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.
Boa noite Fernando! Acho que me expressei mal. Eu sei como faz uma figura animada, eu gostaria apenas de comprimir os tres loops em apenas um. De qualquer modo obrigado pela 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} ________________________________ De: Fernando Mayer <fernandomayer@gmail.com> Para: r-br@listas.c3sl.ufpr.br; Ivan Bezerra Allaman <ivanalaman@yahoo.com.br> Enviadas: Terça-feira, 28 de Agosto de 2012 22:52 Assunto: Re: [R-br] Ajuda em ilustração animada! 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.
participantes (2)
-
Fernando Mayer -
Ivan Bezerra Allaman